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VELOCITY-SPACE PLASMA INSTABILITIES OBSERVED IN A MIRROR MACHINE* 
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Lawrence Radiation Laboratory, University of California, Livermore, California 
(Received January 9, 1961) 


We have previously reported’? our first observa- 
tions of a deliberately induced velocity-space in- 
stability of magnetically compressed plasma, 
adiabatically confined in a mirror machine. Search 
for this instability was prompted by theoretical 
calculations, to our knowledge first performed by 
Rosenbluth,*® predicting that hydromagnetic in- 
stabilities could arise in magnetically confined 
plasmas, the particle distributions of which are 
anisotropic in velocity space. The observations 
reported in this Letter show that plasma instabil- 
ity can be induced in an otherwise stable plasma, 
under the conditions predicted by theory. The 
instability is detected by observing the greatly 
enhanced rate of transport of particles (across 
the flux lines of the confining magnetic field) 
which it causes, as compared to a stable case. 

Aside from its relevance to plasma research, 
the instability we have observed may be of more 
general interest because it can be made to appear 
at very low plasma densities. Thus it might arise 
in situations (such as in some high-current ac- 
celerators) where cooperative effects are nor- 
mally considered unimportant. 

The instability in question is the so-called 
“mirror instability,” predicted to occur in plas- 
mas where p, >p) (the pressure tensor compo- 
nent perpendicular to the lines is greater than 
that parallel to the lines), whenever the local 
value of 8 =(p,)/[(B’/81) +p, ] exceeds a local 
critical value 8,. In these regions where 8 >B,, 
theory predicts that random perturbations in the 
magnetic field should rapidly develop exponential- 
ly (until some nonlinear limit is reached). As 
suggested by its name, one can understand this 


instability by noting that within the plasma any 
local perturbation, depending on whether it in- 
creases or decreases the field, must create 
either a local “magnetic mirror,” or else a trap- 
ping region between two mirrors. If the plasma 
velocity distribution is anisotropic, this in turn 
must lead either to the turning back or to the 
trapping of some particles by the perturbation. 
The diamagnetic effects associated with these 
excluded or trapped particles will in either case 
then act to amplify the initial perturbation, the 
effects becoming regenerative if B>B,. One there- 
fore expects that in unstable regions turbulent 
“rippling” of the magnetic field would develop, 
destroying the adiabaticity of the particle motions, 
and thus leading to enhanced transport effects. 
Clearly, the more anisotropic the distribution, 
the smaller will be the value of 8. Conversely, 
in the mirror machine, confined plasmas which 
have reached diffusion equilibrium with the end 
losses are normally only slightly anisotropic and 
will therefore be characterized by large £,. values, 
typically 0.5 or greater.* To observe the mirror 
instability at smaller values of 8, the anisotropy 
must therefore be made anomalously large. 
Theoretical expressions from which B c can be 
calculated for an “infinite” plasma in a uniform 
magnetic field have been reported.*»*~” To cal- 
culate 8, for a finite plasma in pressure equilib- 
rium in a nonuniform magnetic field (as in the 
mirror machine), a condition derivable from the 
work of Newcomb*® will be used. The condition is 


8<B .=[31(¥p, /p, J(VB/B)|+1]~ stable, (1) 
the gradients being taken in the direction of the 
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magnetic field. This condition is to be applied at 
any point in the plasma. In deriving this expres- 
sion, individual particle motions have been as- 
sumed to be governed by the usual adiabatic in- 
variants, i.e., the magnetic field varies slowly 
in time, and particle orbit diameters are small.® 
In these experiments the confined anisotropic 
plasma is produced by the adiabatic magnetic 
compression of an injected “blob” of low-energy 
plasma. As previously discussed’®™ this proc- 
ess can be used in a mirror machine to generate 
a stably confined plasma of very high electron 
temperature (~25 kev); we induce the instability 
in these experiments by manipulation of the elec- 
tron component of the plasma, the ions playing 
an unimportant role. Because magnetic com- 
pression acts primarily to increase the “temper- 
ature” of only one component of particle motion, 
that perpendicular to the lines of force, anisotropy 
results. It is appropriate therefore to describe 
the final particle distribution by two “tempera- 
ture” parameters, T, and 7), and to take the 
midplane distribution function to be a “two-tem- 
perature Maxwellian,” 


f =exp{-4m[(v,?/T, )+(v )?/T I}, (2) 
within the allowed values of v;/vy. Calculating 
px from (2), (1) can be evaluated for various 
values of the mirror ratio, R =Bmax/B,, as a 
function of position between the midplane (B,) and 
the mirrors (Bmax). For values of t=7/T, <1, 
it is found that the minimum £, occurs at the 
midplane. Figure 1 shows §,(min) as a function 
of ¢, for various values of R. Since in our ex- 
periment R=2, Fig. 1 shows 8,(min)=?, to close 
approximation. Parenthetically, it can be shown 
that the curves of Fig. 1 also represent B-(min) 
for a midplane distribution function which is any 
differentiable function of the variable (v,?/T, 
+v\"/T), not merely the special exponential 
form assumed above. 

To observe the instability at small values of B, 
the anisotropy must be large; this requires a 
large magnetic compression ratio. But since col- 
lisions tend to destroy anisotropy, it is also nec- 
essary that the initial plasma density be low enough 
that collision effects are unimportant during the 
compression, but not so low that B lies below its 
critical value. Therefore, in order to produce 
the instability, it was necessary to reduce the 
initial density about a factor 10 below that nor- 
mally used." Magnetic compression ratios 
By /Bj= 10° or more were used, so that ¢t =6,.=107°. 
The peak magnetic field, Br, was 10 kilogauss. 
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eral values of the mirror ratio R. 


Figure 2 shows the experimental arrangement 
used. Scintillation probes were used to detect the 
fluxes of radially and axially escaping electrons. 
The plasma conditions, which determined whether 


or not instability occurred, were varied by chang- 


ing (a) the initial magnetic field, (b) the intensity 
and the timing (relative to the pulsed magnetic 
field) of the injected plasma burst, and (c) the 
residual gas pressure. 

Most of the observations consisted of correla- 
ting the intensity of signals received at the radial 
scintillation probe (see Fig. 2), with changes in 
the plasma or compression parameters. The 
most striking feature of the signals received by 
this probe was the extreme variation in intensity 
between stable and unstable cases. Whereas 
strong signals, several orders of magnitude above 
multiplier noise, were detectable in the unstable 
cases, signal levels dropped below detectability 
when the plasma was rendered stable. Figure 3 
illustrates the effect of changing plasma density 
(by varying the plasma-source-pulse-capacitor 
charging voltage). Although the signal received 
by the end scintillation probe increases mono- 
tonically until detector saturation occurs, the 
time-integrated radial probe signal increases 
only up to a critical density and then falls rapidly 
to zero with further increase in density. At this 
density we believe that collisional randomization 


becomes important and suppresses the instability. | 


This interpretation is apparently borne out by 
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7 experiments in which a neutral “scattering gas” 
er (helium) was admitted to the chamber. As the 
helium density was raised above about 3 x10” 


- | 
“4 cm~’, the ratio of radial-to-axial signal dropped 
g rapidly to zero. From scattering theory one can 
: calculate that, for electrons of mean energy of 
i 20 kev, this density of helium scattering centers 
__ | would smooth out all anisotropies in about 300 
al microseconds, i.e., in a time comparable to the 
| compression time of 500 microseconds. 

The velocity of plasma transport (across the 
magnetic field lines) caused by the instability 
yf was estimated by three methods: (1) Calculation 

i from the rate of collapse of flux tubes during 
ove ; compression; the plasma collected during com- 
le pression must possess at least enough velocity 


y : to “paddle upstream” to the detectors. (2) Time 





3 _ correlation between signals received by scintil- 
y | lation probes at two different radial positions. 
' ie (3) Observation of the effect of a movable obstacle 
i || Which projects into the chamber at a point dia- 
: metrically opposite the radial probe. By inserting 
this probe into the chamber to a radius less than 
_ _ that of the radial scintillator probe we decreased 
dly | its signals; from these data and from the known 
is mean precession time of electrons around the 
on chamber (about 8 microseconds), the radial trans- 
lity. | port velocity can be calculated. The mean veloci- 
F ties deduced by the above three methods ranged 
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between 2 x10* and 2x10° cm/sec. These veloci- 
ties are more than 10’ times the “classical” dif- 
fusion velocities appropriate to the plasma density 
and temperature, and are at least 10* times the 
upper limits experimentally deduced" for the 
stable, higher density, plasmas studied earlier 

in the same apparatus. 

Associated with the presence or absence of the 
instability one expects substantial changes in the 
radial density distribution of the compressed 
plasma. By using a movable scintillator probe 
located just outside one of the mirrors, the radial 
distribution of the escaping electrons was meas- 
ured; from such measurements radial density 
distributions can be inferred. Figure 4 shows 
these data for two cases: (1) plasma density in 
the range where instabilities are observed; 

(2) “low” plasma density, i.e., B<B,- The dif- 
ferences are striking. Similar constricting ef- 
fects were found upon increasing the initial field 
(lowering By /B;, therefore reducing anisotropy), 
even though the theoretical radial compression 
itself was thereby reduced. 

Approximate quantitative agreement between 
the theoretical 8, and experimental 8 values was 
obtained from the following considerations: The 
theoretically predicted axial variation of plasma 
density [calculated by using Eq. (2)] was compared 
with the experimental data (obtained by axial move- 
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pulse-capacitor charging voltage. 


ment of a scintillator probe about the midplane) 

to determine the anisotropy. In the cases analyzed 
t=107°, so that, from Fig. 1, B-=10™°. From 
calorimetric measurements’ of the energy flux 
carried radially by the escaping plasma, and from 
the measured plasma volume, experimental values 
of 8 could be deduced. For a plasma that had dif- 
fused out to the scintillator, these ranged between 
1x10 and 4x107°, i.e., between about 1 and 4 
times the “theoretical” value of B,. (The cor- 
responding plasma density was therefore of order 
10° cm~*.) However, the value of 8 prior to ex- 
pansion of the plasma out to the detector should 

be larger than its value after expansion. This 
means that we believe that whenever the pressure 
of the compressed plasma was initially only slight- 
ly above the critical value, only the most aniso- 
tropic components were able to diffuse out to the 
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FIG. 4. End escaping flux as a function of radius 
for an initial magnetic field of 0.0 gauss for medium 
and low densities. 


radial probe, the pressure of the rest of the 
plasma dropping below £8, (so that the plasma 
became “frozen” to the field lines) before it could 
reach the detector. We have obtained experimen- 
tal evidence for this interesting selective trans- 
port effect by method (3), above, used to measure 
radial transport velocity. These measurements 
showed a clear separation of the radial probe 
signals into two groups: (1) a “fast” group, ap- 
pearing within about 500 usec of the initiation 
of compression; (2) a “slow” group, appearing 
several milliseconds after peak field, i.e., not 
until inner field lines of the collapsing field had 
expanded out to the detector. For the “slow” 
group a relative motion of the obstacle and de- 
tector of only 0.005 inch was sufficient to extin- 
guish most of the detected signal. The “late” 
plasma particles appeared therefore to be frozen 
to the expanding field lines. In support of this 
interpretation, it was found from the axial varia- 
tion of the radial flux about the midplane that only 
the “fast” group was characterized by a very 
large ratio of p,/py. These measurements also 
apparently resolve an old mystery; we had earlier 
discovered’® that under certain conditions the 
outermost energetic electrons of an otherwise 
stable compressed plasma core could apparently 
diffuse rapidly to the chamber walls, producing 
x rays. Their rate of radial transport was much 
greater than that predicted by classical theory. 
We now believe that this effect was early evidence 
of the “mirror” instability here discussed, but 
occurring only in the outer, low-density, portion 
of the plasma core. 

In summary, we believe that we have established 
the existence and some of the manifestations of a 
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hydromagnetic velocity-space instability of the 
“mirror” type. In our experiments, using high- 
energy electrons, it has been observed to cause 
an enhanced rate of transport of particles across 
field lines. It may also cause other effects, such 
as more rapid thermalization of the plasma, or 
the emission of plasma radiation. These effects, 
or those which might arise when the instability is 
associated with anisotropies in the ionic compo- 
nent of the plasma rather than the electronic 
component, we have not investigated. Although 

it is possible to avoid the instability by control- 
ling the plasma parameters, we have shown that 
it may occur in plasma compression experiments. 
It also seems likely that any method of creating 

a hot plasma which results in anomalously large 
anisotropies (such as the method of injecting 
highly directed high-energy particles), may stim- 
ulate this instability during buildup, even if the 
final plasma state aimed for is stable. Finally, 
since this instability feeds on perturbation in 
energy density, and is not of electrostatic origin, 
it should persist in the limit of very low particle 
densities, provided the anisotropy becomes cor- 
respondingly large. For this reason, it (or re- 
lated instabilities) might occur in some astro- 
physical situations or in particle accelerators, 
even though the conditions are such that coopera- 
tive effects would normally be considered unim- 
portant. 
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SUPERCONDUCTIVITY IN Nb,Sn AT HIGH CURRENT DENSITY IN A 
MAGNETIC FIELD OF 88 kgauss 


J. E. Kunzler, E. Buehler, F. S. L. Hsu, and J. H. Wernick 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received January 9, 1961) 


We have observed superconductivity in Nb,Sn 
at average current densities exceeding 100000 
amperes/cm? in magnetic fields as large as 88 
kgauss. The nature of the variation of the criti- 
cal current (the maximum current at a given 
field for which there is no energy dissipation) 
with magnetic field shows that superconductivity 
extends to still higher fields. Existing theory 
does not account for these observations. In addi- 





tion to some remarkable implications concerning 
superconductivity, these observations suggest 
the feasibility of constructing superconducting 
solenoid magnets capable of fields approaching 
100 kgauss, such as are desired as laboratory 
facilities and for containing plasmas for nuclear 
fusion reactions.*»? 

The highest values of critical magnetic fields 
previously reported for high current densities 
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have been less than 18 kgauss.*»* However, 
values as large as 30 kgauss have been reported® 
for low current densities. 

The discovery of the existence of Nb,Sn with 
its high critical temperature of 18.0°K by 
Matthias, Geballe, Geller, and Corenzwit® im- 
mediately suggested that the critical field might 
be high. The existence of superconductivity in 
Nb,Sn at high fields is also indicated by the re- 
cent susceptibility observations of Bozorth, 
Williams, and Davis.’ 

Nb,Sn is very brittle and is thus not readily 
fabricated into specified configurations. How- 
ever, a method of minimizing this difficulty was 
developed during the current investigation. Two 


distinctly different types of samples were studied. 


The first type consisted of rectangular rods cut 
from an ingot which was prepared by first sinter- 
ing a stoichiometric mixture of Sn and Nb pow- 
ders at 1800°C and then melting the compact in a 
zirconia crucible in an argon atmosphere at about 
2400°C. X-ray diffraction data confirmed the 
existence of essentially a single homogeneous 
phase of the 8-wolfram structure. The rods 
were about 2 cm long and had thicknesses varying 
from 0.025 cm to 0.063 cm. The second type 
consisted of small niobium tubes which contained 
thin cores of Nb,Sn. These were prepared by 
packing 0.6 cm o.d.—0.3 cm i.d. niobium tubes 
with either a mechanical mixture of powdered 
Nb,Sn + a 10 weight percent excess of powdered 
Sn or similar mixtures of unreacted Sn and Nb 
powders. The ends of the tubes were closed with 
Nb plugs and the sealed tubes were mechanically 
reduced in size to 0.038 cm o.d. Each tubular 
sample was heated to a prescribed temperature 
between 970°C and 1400°C for a predetermined 
period of time, as long as 24 hours. 

Current and potential leads were attached to 
each of the samples using an ultrasonic soldering 
iron and indium solder. The potential leads were 
about 1 cm apart. The samples were mounted 
transverse to the axis of a solenoid magnet cap- 
able of a maximum steady-state field of 88 
kgauss. Observations were made by selecting a 
fixed value for the magnetic field and slowly in- 
creasing the current through the sample until a 
potential drop across the sample of a few hun- 
dredths of one microvolt was observed. The re- 
sults for three sizes of parallelepiped samples 
cut from the Nb,Sn ingot are shown in Fig. 1. 

In addition to the high magnetic fields and cur- 
rent densities for which superconductivity exists 
in Nb,Sn, there are two other interesting fea- 
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for three sample sizes at 1.5°K and 4.2°K. The mag- 
netic field was perpendicular to the current direction. 
Each experimental point represents the maximum cur- 
rent, at the value of magnetic field indicated, for 
which no voltage drop along the sample was observed, 
the smallest detectable voltage being a few hundreths 
of one microvolt. 


tures of the data of Fig. 1. First, the average 
current density does not scale as the perimeter 
of the cross-sectional area as is expected for a 
“soft” superconductor, nor as the cross-sec- 
tional area as is observed for “hard” supercon- 
ductors, such as mechanically deformed Mo,Re,* 
but is about midway between the two types of be- 
havior. This variation with sample size suggests 
that appreciable current is being carried both by 
the surface and by “filaments.” Secondly, the 
variation of the maximum permissible applied 
magnetic field with temperature, at constant 
current, is far from parabolic. Although 7, is 
about 18°K for Nb,Sn,° the magnetic field in- 
creases more than 50% between 4.2°K and 1.5°K. 
We have not been able to account for these ob- 
servations in terms of existing theory. 

Critical current vs applied magnetic field 
curves for several characteristic samples of 
“Nb,Sn” cores clad with Nb are shown in Fig. 2. 
It is of considerable interest that all of the curves 
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FIG. 2. Critical current vs applied magnetic field 
for Nb-clad cores of “ Nb;Sn.” The o.d. of the cores 
was about 0.015 cm and the o.d. of the Nb jackets was 
about 0.038 cm. “+ 10% Sn’ in the table legend means 
10 wt. % more Sn than is required to form Nb;Sn as- 
suming no reaction with the Nb tube. The magnetic 
field was perpendicular to the current direction. Each 
experimental point represents the maximum current, 
at the value of magnetic field indicated, for which no 
voltage drop along the sample was observed, the 
smallest detectable voltage being a few hundreths of 
one microvolt. 


represent significantly higher average current 
densities at the highest fields than the samples 
cut from the Nb,Sn ingot. The single point shown 
on Curve 1 in Fig. 2 represents an average cur- 
rent density nearly 50 times higher than the data 
of Fig. 1. The average current density of this 
“Nb,Sn” core exceeds 100000 amperes/cm?. The 
Nb tubes do not contribute since it is found that 
annealed niobium becomes normal in fields well 
below 20 kgauss. A comparison of Curve 3 of 

Fig. 2 with Curve 5 suggests that a stoichiometric 
mixture of Nb and Sn powders is preferable to an 





excess of Sn. A comparison of these curves with 
Curves 1 and 2 indicates that it is preferable to 
start with the elemental powders, with or with- 
out an excess of Sn, rather than with powdered 
Nb,Sn + excess Sn powder. Also, the curves of 
Fig. 2 and similar data for heat treatment at 
temperatures as high as 1400°C show that a low- 
er temperature of heat treatment consistently 
yields material capable of higher critical current 
density. This rapid trend with temperature sug- 
gests that still higher current densities might be 
obtained at still lower reaction temperatures. 

The data of Fig 2 and data from related sam- 
ples suggest that the difference in properties 
between pure dense Nb,Sn and the Nb-clad sam- 
ples is not principally one of chemistry but ra- 
ther is in some manner associated with the physi- 
cal state of the “Nb,Sn” in the core. However, 
the critical temperature of the fused dense Nb,Sn 
was found to be 17.8°K (when extrapolated to the 
limit of zero current through the sample) while 
that of the clad sample having the highest critical 
current density (Curve 1 of Fig. 2) was found to 
be 17.9°K. Although the pertinent physics of the 
situation is not yet clear, it is tentatively con- 
cluded that the conditions of preparation of the 
clad samples are such as to lead to a structure 
containing large numbers of “filaments.” It is 
fortunate that these conditions are compatible 
with those which are expected to permit fabrica- 
tion of the material into high-field superconduct- 
ing electromagnets. 

We are indebted to B. T. Matthias for encour- 
aging the investigation of Nb,Sn for high-field 
studies. Encouragement and stimulation from 
R. Kompfner and M. Tanenbaum contributed 
significantly to this work as did assistance from 
H. Dail and C. V. Wahl. 
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SUPERCONDUCTORS AS QUANTUM DETECTORS FOR MICROWAVE AND 
SUB- MILLIMETER-WAVE RADIATION* 


E. Burstein, D. N. Langenberg, and B. N. Taylor 
University of Pennsylvania, Philadelphia, Pennsylvania 
(Received January 6, 1961) 


Giaever' and Nicol, Shapiro, and Smith? have 
reported experiments on electron tunneling 
through metal oxide barriers between a metal 
and a superconductor, M|B|S, and between two 
superconductors, S,|B|S,. The tunneling pro- 
cess which involves “unpaired” charge carriers 
can be simply described in terms of a quasi- 
particle energy band picture in which the super- 
conductor is represented, at temperatures be- 
low its transition temperature, as an intrinsic 
semiconductor with a temperature- dependent 
energy gap, Eq= 2€,, and equal effective masses 
for electrons and holes [ Fig. 1(a) |. 

For small bias voltages, V,, such that |eVp| 
is less than E¢/2 in the case of M|B|S and less 
than (Eg +Eq@ )/2 in the case of S,|B|S,, the 
tunneling current depends on the density of ther- 
mally excited carriers. For larger bias voltages, 
the tunneling current is essentially independent 
of the density of thermally excited carriers and 
results in the “injection” of electrons and holes 
into the conduction and filled bands of the super- 
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conductor. The concept of “photo-injection” of 
carriers by optical excitation across the energy 
gap of the superconductor follows quite naturally. 
These considerations suggest the use of “low- 
voltage” tunneling in M|B|S and S|B|S structures 
for the quantum detection of microwave and sub- 
millimeter-wave radiation in a manner which is 
analogous to the detection of visible and near in- 
frared radiation by p-” junctions in semicon- 
ductors [Fig. 1(b)]. Such superconductor photo- 
detectors would have long-wavelength limits, 
Amax’ determined by the energy gap of the su- 
perconductor used. Thus, for aluminum \max 
= 3.9 mm corresponding to E¢=3.2x10~* ev, and 
for lead max = 0.46 mm corresponding to E,=2.7 
x 107% ev.» 
The magnitude of the tunneling current due to 
optically excited carriers will depend, among 
other things, on the lifetime of the carriers. The 
processes which limit the lifetime include (a) the 
recombination of electrons and holes across the 
gap, (b) the pairing of electrons in the conduc- 
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FIG. 1. (a) Quasi-particle 
energy band diagrams for an 





M|B|S structure and an S,|B|S, 
structure in which the two 
superconductors are the same. 











(b) Tunnelling of optically ex- 
cited carriers in M|B|S and 
S,|B|S, structures under “low 
bias voltage’ conditions. 


a 




















92 




















Se Tee eA ee 












VoLtuME 6, NUMBER 3 





PHYSICAL REVIEW LETTERS 


FEesBruARY 1, 1961 





tion band and the corresponding pairing of holes 
in the filled band, and (c) in the case of M|B|S 
structures, the simultaneous tunneling of elec- 
trons and holes into the metal from the super- 
conductor and their subsequent transition to the 
Fermi surface of the metal. At the present time, 
estimates for the lifetimes for these processes 
which may be radiative or nonradiative are not 
available. We have carried out an estimate of 
the electron-hole radiative recombination life- 
time, Tope for lead at 2°K using the theory of van 
Roosbroeck and Shockley*® (Fig. 2). The depend- 
ence on frequency of 7x, where n+ix is the com- 
plex refractive index, was derived from theoret- 
ical expressions for 0,+i0,, the complex con- 
ductivity, obtained by Mattis and Bardeen,‘ while 
the density of thermally generated carriers, n;, 
was calculated from the quasi-particle density 

of states of superconductors given by the Bar- 
deen, Cooper, and Schrieffer theory® and a value 
of N(0) = 2.2 x10”8/cm* which Gold® proposed for 
the density of states at the Fermi surface of lead 
in the normal state. The theoretical values of 

0, and o, are in reasonable agreement with the 
experimentally determined values obtained by 
Glover and Tinkham’ and Ginsberg and Tinkham® 
but do not show any of the structure near the 
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FIG. 2. The dependence of n *« and of Un*x, which 
is proportional to the radiative recombination rate per 
unit frequency interval, on u=liw/kT for lead at 2°K. 








band gap which is present in the experimental 
spectrum. We obtain n;=5.5x10"*/cm® and 
Topt = 0.4 sec. It is, therefore, very unlikely 
that the radiative recombination of electrons and 
holes is the dominant process which determines 
the lifetime. Measurements of the tunneling cur- 
rent due to optically and electrically injected 
carriers should provide an effective means for 
obtaining information about the lifetime, and, in 
appropriate configurations, about the diffusion 
length of the unpaired carriers in superconduc- 
tors. 

B. Rosenblum of RCA Laboratories has pointed 
out that photoexcited carriers in the metal of the 
M|B\S structure may also contribute to the tun- 
neling current. This would extend the photore- 
sponse to longer wavelengths with A»,a, con- 
trolled by the bias voltage. The magnitude of the 
photoresponse of the metal would be expected to 
be much smaller than that of the superconductor 
since the lifetime of the excited carriers in the 
metal is very much shorter than the correspond- 
ing lifetime in the superconductor. 

The metal and superconductor films used in the 
superconductor photodetectors must be much 
smaller than the wavelength of the radiation in 
the films in order to decrease the reflection and, 
thereby, to increase the absorption of the inci- 
dent radiation by the superconductor. An S,|B|S, 
structure is an ideal configuration for such de- 
tectors. On the one hand, it avoids the annihila- 
tion of electrons and holes by simultaneous tun- 
neling of electrons and holes from the supercon- 
ductor into the metal which can occur in M|B|S 
structures, and on the other hand it allows an 
optimum absorption of the incident radiation by 
the superconductor layers. The absorption of 
radiation by an S,|B|S, structure can be calcula- 
ted by considering the composite structure as a 
single superconductor film. Curves of the ab- 
sorptivity, A, and of A/(hw/E ~) which is pro- 
portional to the number of photons absorbed per 
unit intensity of incident radiation, are plotted 
as a function of hw/E in Fig. 3 for lead films 
having thicknesses of 5, 15, and 150 angstroms. 
We see that the absorptivity of the 15-angstrom 
film approaches the maximum value of 0.5 at 
hw/EG=2.5. However, the 150-angstrom film 
still exhibits an appreciable absorptivity, having 
a value of 0.15 at he/E ¢~3. The photoresponse 
of the S,|B|S, structure may be expected to differ 
somewhat from the photon absorption spectrum 
as a result of two effects: (a) a dependence of 
the lifetime of the excited carriers on energy, 
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FIG. 3. Absorptivity and photon absorption spectra 
for superconducting lead films of 5, 15, and 150-ang- 
strom thicknesses. 


and (b) an increase in quantum efficiency for 
high-energy photons due to secondary impact 
ionization processes, as has been observed in 
the case of photoconductors.® The latter effect 
would increase the photoresponse at the higher 
frequencies and, thereby, would compensate for 
the decrease in photoresponse resulting from the 
decrease in the number of photons absorbed per 
unit intensity of radiation with increase in fre- 
quency. However, the photoresponse may also 
be expected to exhibit a sharp drop-off at fre- 
quencies greater than the “collision” frequency 
of the carriers (not shown in Fig. 3) where o, 
exhibits a 1/w? dependence. 

Operating conditions for the superconductor 
photodetector are straightforward. The super- 
conductor should be maintained at a temperature 
well below its transition temperature in order to 
diminish the thermal generation of carriers. 
Under these conditions (kT <<E) the thermally 
excited carriers obey Maxwell- Boltzmann sta- 
tistics and the equilibrium density of excited 
carriers is, to a good approximation, given by 
n;=N(0)(nE (kT) exp(-E G/2kT)/2 which has the 
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same temperature dependence as a one-dimen- 
sional intrinsic semiconductor. For optimum 
response the temperature is lowered to a point 
where the rate of optical generation of carriers 
by background radiation is greater than the rate 
of thermal generation.’ Under these conditions 
the fractional change in carrier density produced 
by a small radiation signal is given by An/n; 
=G,/Gp, where Gg is the rate of generation by 
the radiation signal and Gp is the rate of genera- 
tion by the background radiation accompanying 
the radiation signal. The frequency response of 
the photodetector should be limited only by the 
lifetime of the excited carriers. It should also 
be noted that the photodetector can be operated 
with bias voltages of either polarity, as con- 
trasted to the semiconductor photodiode. 

It is a pleasure to acknowledge the valuable 
discussions we have had with R. Amado, M. 
Cohen, R. Parmenter, A. Rose, B. Rosenblum, 
A. Rothwarf, P. Stiles, and R. deZafra. This 
research was supported in part by the U. S. 
Atomic Energy Commission. 
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STIMULATED OPTICAL EMISSION FROM EXCHANGE-COUPLED IONS OF Cr*** IN ALO, 


if Irwin Wieder and Lynn R. Sarles 
re Varian Associates, Palo Alto, California 
Z (Received December 19, 1960) 





















| The advantages of chromium ions in aluminum little more complicated. The fact that it exhibits 
oxide crystals (ruby) for double resonance ex- enhanced emission while the R, line in this sam- 
periments have been pointed out'»? and utilized ple does not is a result of the relatively unpopu- 
- in several microwave-optical experiments.* ° lated terminal state for the 7010A transition 
Recently, stimulated emission® and coherence’ (negligible optical absorption for this wavelength 
of the characteristic “R,” fluorescence at at 77°K). Since the number of atoms radiating 
6943 A has been observed in pink ruby (~0.07% 7010A light saturates at intermediate energy 
Crt** by weight). We wish to report the obser- (Fig. 2), it appears that the total number of equi- 
| vation of stimulated emission at wavelengths valent ion pairs which have this line in their 
y 7010 A and 7040 A from transitions in red ruby spectrum is excited at this energy. This number 


(~0.7% Cr*** by weight) which arise® from ex- 
change coupling between neighboring chromium 


ions. 
The observations were made by monitoring the — t— 1 Millisecond 
' 


| | 


fluorescence of a 4-in. x13-in. red crystal rod 
following an intense excitation flash. By virtue 
of the broad optical absorption bands and the ef- 
i} ficient radiationless transfer to sharp excited 
states, these levels accumulate an excess popu- 
lation which may satisfy the requirements® for 
stimulated optical emission. At high excitation 
energy the fluorescence at 7010 A and 7040 A 
exhibits the initial decrease from the natural 
radiative lifetime which is characteristic of 
stimulated emission. However, the decay curves 
show significant differences from the previous 
2 results with the R, line. Figure 1 is a compari- 
son of the R, line in pink ruby, the 7010A line in 
if red ruby, and the R, line in red ruby, all under 
i similar conditions of excitation at 77°K. The de- 
cay curves at low excitation are all similar in 
appearance to Fig. 1(c), and in two cases the 
8, : shapes change progressively with increasing ex- 
citation until they reach the curves shown in 
59). Fig. 1(a) and Fig. 1(b). A plot of the areas under 
the fluorescence curves of R, and 7010 Aasa 
function of flash energy is shown in Fig. 2. 
Our interpretation of the data is as follows: It 
is apparent that the R, line in red ruby does not 
exhibit stimulated emission. This follows from 
the fact that the absorbed flash energy is not 
sufficient to overcome the ground-state popula- 
tion, and is corroborated by the approximately 
linear relationship between R, and flash energy 
in Fig. 2. (The increased linewidth in red ruby 
would not change conditions enough to prevent FIG. 1. Fluorescence decay curves of pink and red 
stimulated emission from occurring.) The be- ruby taken at 77°K under approximately the same high 
havior of the 7010A line (and 7040A line) is a excitation conditions. 





oe 


(a) Pink 6934 A 











(b) Red 7010 A 





(c) Red 6934 A 
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FIG. 2. A plot of area under fluorescence decay 
curves for R, and 7010 A in the same red ruby at 77°K 
as a function of excitation flash energy. 


is not sufficiently greater than the threshold 
number needed for stimulated emission to exhibit 
the degree of lifetime shortening observed in 
pink ruby [compare Fig. 1(a) and Fig. 1(b) |. 


These experimental results also yield informa- 
tion on the excitation processes of the coupled 
ions. Since more than half of the single ions are 
in the ground state even at high excitation, direct 
excitation via the broad absorption bands is not 
likely for those ion pairs near the center of the 
rod. On the other hand, thermal transfer of ex- 
citation between single ions and paired ions is 
ruled out by the saturation behavior of the 7010A 
line. One possibility which seems promising is 
radiative excitation by trapped 6934A light. As 
a preliminary check we have observed absorp- 
tion lines in the red sample which overlap the 
6934A line and have also observed the 7010A 
fluorescence when the red rod is excited by a 
pulse of 6934A light from a pink sample. Fur- 
ther experiments to clarify the nature of the ex- 
citation processes of exchange-coupled ions are 
in progress. 
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SIMULTANEOUS OPTICAL MASER ACTION IN TWO RUBY SATELLITE LINES 


A. L. Schawlow and G. E. Devlin 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received December 19, 1960) 


We have obtained pulsed maser oscillations 
simultaneously at wavelengths of 7009 A and 
7041 A in concentrated ruby. The lines involved 
are two of the strongest of the satellite or neigh- 
bor lines,’ which have been known and studied 
for many years but are not yet fully analyzed.?"° 
It is known’ that these lines disappear at low con- 
centrations, but appear strongly in emission at 
chromium concentrations of the order of 0.5%. 
In absorption, they are weak and disappear en- 
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tirely at low temperatures. 

Since these lines terminate on states which are 
far enough above the ground states to be empty 
at low temperatures and since their upper states 
can be populated easily by white light illumina- 
tion, it was proposed that the strongest of them 
be used for an optical maser.® A very low esti- 
mate of the fluorescence efficiency (0.1%) was 
used, but with the advantage of an empty lower 
state, pulsed operation seemed possible at low 
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temperatures. Now that the fluorescence effi- 
ciency for the main R lines is known to be nearly 
unity,” it is apparent that much less pumping 

light is needed. Moreover, recent investigations 
of the 7009A line in absorption have shown that 

its lower state is about 100 cm™ above the ground 
state. Thus liquid hydrogen or neon temperatures 
would be sufficiently low to depopulate this lower 
state. 

However, operation should also be possible at 
higher temperatures. This is expected because 
even at 77°K the line is so strong in emission, 
but so weak in absorption. It seems possible that 
the very numerous individual chromium ions ab- 
sorb the pumping radiation and transmit it to the 
much rarer paired ions which emit the N lines.® 
Thus even if we have to invert the population of 
the N-line levels by pumping to empty the lower 
level, this can be accomplished with only enough 
pumping light to excite a small fraction of the 
isolated ions. 

Maser operation was sought first in a ruby rod 
containing about 0.5% Cr,O,, having a diameter 
of 2 mm anda length of 4cm. The ends of the 
rod were silvered, one of them completely, and 
the other to a transmission of about 1%. The rod 
was immersed in liquid nitrogen in a small clear 
Dewar flask, with its upper end supported in a 
cone which allowed emergence only of light 
coming through the silver film. Pumping radia- 
tion was supplied by a General Electric FT524 
xenon flash tube. When a 400-microfarad capa- 
citor bank was charged to 3500 volts and dis- 
charged through the flash tube, the relative in- 
tensity of the 7041A line was greatly enhanced 
over the other lines in the spectrum. At about 
3700 volts it was joined by the 7009A line. Fig- 
ure 1 shows these two lines photographed in a 
Single flash at about 3900 volts. For compari- 
son, part of the spectrograph slit was covered 
by a concentrated ruby at room temperature. 
When the latter was illuminated by a microscope 
lamp, it gave a comparison spectrum similar to 
that of the ruby below the maser threshold except 
for a small wavelength shift of all lines. 

The observation of enhancement thresholds is 
a strong indication of maser action.* Further 
confirmation was obtained by observation of 
parallelism in the output radiation, as shown by 
a sharply defined spot on a distant surface. The 
beam width, well above the pumping threshold 
energy, was about 1 degree. This approximate 
parallelism implies considerable spatial coher- 
ence in the output.® 





(a) 


(b) 





FIG. 1. (a) Spectrum of a concentrated ruby at room 
temperature with a 10-second exposure. (b) Spectrum 
of a ruby rod at 77°K showing N-line maser oscilla- 
tions as described in the text, with a single flash last- 
ing about 1/2000 second. All lines are shifted slightly 
because of the temperature difference between (a) and 
(b). The wavelength scale is in angstroms. 


It may seem surprising that optical maser ac- 
tion occurs at two frequencies simultaneously. 

If the upper levels of the two lines are identical, 
or relax rapidly to each other, as with the E and 
2A levels of the isolated Cr** ions, this will not 
happen. In such a case, when maser action starts 
to empty one of the levels rapidly, it prevents 
the additional buildup of population needed for 
maser action from the other level. However, it 
had already been suspected that the 7041A line 
has an origin not closely connected with the 
7009A and other neighbor lines. The evidence 
for this different origin is that the relative in- 
tensity of the two lines varies erratically from 
sample to sample. The intensity of the 7009A 
line varies with concentration in much the same 
way as most of its neighbors, whereas the 7041A 
line is usually weaker, but occasionally brighter 
than the 7009A line.’° Occurrence of simulta- 
neous maser action is confirmation of their se- 
parate origin. 

The pumping power needed for this sample is 
no less than that required for R,-line maser ac- 
tion in dilute ruby.**® There are many ways in 
which the power required can be reduced. Of 
these, the simplest is to cool it further, and 
this is being studied. Preliminary attempts pro- 
duced too much heat in the sample, so that its 
temperature rose far above that of the bath. In- 
deed, heating may have increased the threshold 
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power even when the sample was in liquid nitro- 
gen. It is believed that more careful filtering of 
the light source to remove all but the useful 
pumping radiation will permit low-temperature 
operation at greatly reduced pump power. 

It has thus been shown that optical maser oper- 
ation at the wavelength of the 70094 line (N,) 
does occur in concentrated ruby as predicted.® 
Moreover, simultaneous maser oscillations at 
7041 A (N,) show that these two lines are not 
closely related, and probably come from differ- 
ent kinds of chromium ion pairs. 

We wish to thank W. L. Bond for polishing the 
ruby rod, D. F. Nelson for lending us the flash 
tube mount and power supply, and E. M. Kelly 
for silvering the ends of the rod. We are also 
indebted to W. S. Boyle for pointing out the im- 
portance of sample heating by the flash, and to 
A. M. Clogston and D. L. Wood for helpful dis- 
cussions on the origin and properties of the 
satellite lines. 
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INTERPRETATION OF THE Fe*” ISOMER SHIFT 


L. R. Walker, G. K. Wertheim, and V. Jaccarino 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received January 3, 1961) 


We have made a systematic study of the Méss- 
bauer effect’ of Fe*”’ in di- and tri-valent iron 
compounds and in d-group metals. The observed 
“isomer shift’” measures the total s-electron 
density at the nucleus. The shift in compounds 
is shown to depend mainly upon the 3d configura- 
tion of iron involved and to a lesser extent upon 
the “chemical” bond. The Hartree-Fock calcula- 
tions of Watson* on the various 3d configurations 
of iron are combined with the data on the shifts 
in the most ionic compounds to obtain a calibra- 
tion of the shift in terms of s-electron density. 
This enables us to estimate the difference in 
charge radius of the ground state and isomeric 
state of Fe*’. We associate the shift of an Fe 
solute ion in a d-group metal with the addition of 
some fraction of a 4s electron to an Fe 3d’ con- 
figuration; an estimate of the 4s-electron wave 
function density at the nucleus from the Fermi- 
Segré-Goudsmit (FSG) formula‘ enables this 
fraction to be determined. 

Kistner and Sunyar® first observed in the recoil- 
free emission and resonant absorption of the 
14.4-kev nuclear gamma ray of Fe*”’, that the 
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energies of the emitted and absorbed gamma rays 
were noticeably different if the emitter and ab- 
sorber were two dissimilar lattices containing 
iron. The origin of this effect is as follows. 

The ground and isomeric levels of the nucleus 
have different effective charge radii; the elec- 
trostatic interaction with the electronic charge 
is then different in the two states and the gamma- 
ray energy is consequently changed (relative to 
its value for a point nucleus) by an amount pro- 
portional to the s-electron density at the nucleus. 
If the s-electron density is different for the ab- 
sorber and emitter, the difference in gamma- 
ray energies, E,-E,, is defined as the isomer 
shift. 


EB -E,= Ze |R, “a, 1y(0) 1-1 y(0), 17], (1) 





where Rj, and Ro, are the radii of the isomeric 
and ground states and |4(0), 1? and |y(0),1? are 
the total s-electron densities at the nucleus for 
absorber and emitter, respectively.® It is to be 
noted that the sign of the shift has its origin in 
the fact that the level energy is lower the more 
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compact the charge distribution. 
An examination of the restricted Hartree-Fock 
calculations of Watson® shows that there are 
significant differences in the value of |¥,,(0)|? 
for different 3d configurations of iron; the change 
in |¥,s5(0)|? and |y,5(0)|? is substantially smaller. 
It is perhaps worthwhile to point out that un- 
restricted Hartree-Fock calculations, which ex- 
hibit the imbalance of up and down spin density 
(exchange polarization effect), appear to give the 
same total density at the nucleus as the restricted 
procedure.’ The variation of |¥,,(0)|? is such as 
to correspond to different degrees of shielding 
of 3s by 3d electrons. To calibrate the observed 
shifts in terms of total s-electron density, we 
associate the difference in the shifts for the most 
ionic Fe?* and Fe** compounds with the differ- 
ence in Watson’s values of 2) ,=3|%,,(0)|? for 3d° 
and 3d° configurations. When an Fe atom is in- 
troduced into a d-group host metal, 2) |¥%)<(0)!? 
will also contain a contribution from 4s conduc- 
tion electrons. Hartree-Fock calculations do 
not exist for single 4s electrons outside 3d” con- 
figurations, but the FSG* formula, combined with 
the known term value of 3d”4s configurations,® 


provides an entirely adequate estimate of |y,,(0)|*. 


In Fig. 1 is shown a possible interpretation of 
the observed shifts making use of the above ideas. 
The total s-electron density in atomic units is 
plotted as ordinate. Watson’s values for 
22) n=1!%ps(0)!? for the Fe 3d” configurations 
from n =4 to n =8 are indicated on the left. The 
scale on the right for Méssbauer center- of-gravi- 
ty shifts relative to stainless steel is established 
by identifying the shifts in the most ionic Fe** 
and Fe** with Watson’s densities for 3d° and 3d°, 
respectively. The solid straight lines represent 
s-electron densities for hypothetical 3d"4s* con- 
figurations. They are drawn on the assumption 
that the density for such a configuration is of the 
form |y(3d") |? +x1y,,(0)1?, where |¥,.(0)|? is 
calculated from the FSG formula for a single 4s 
electron outside the 34” configuration. This 
assumes no screening of inner s electrons by 
the single 4s electron.® Curves for the configur- 
ations 3d®°-*4s* and 3d’-*4s* are indicated by 
dashed lines and were obtained by extrapolation 
and interpolation of the FSG formula. The ex- 
perimental data are given in Table I. 

The shifts for Fe in various metals have been 
represented as horizontal lines of a length suf- 
ficient to cover what appear to be the most plaus- 
ible configurations. Since for Fe in Fe metal 
there are certainly 8 electrons to be accounted 
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A= 4s ELECTRON CONTRIBUTION IN % 


FIG. 1. A possible interpretation of the Fe’ Méss- 


bauer isomer shifts in various solids. The total s-elec- 


tron density is plotted as a function of the percentage 
of 4s character for various d-electron configurations. 
The reasons for placing the experimental data on given 
theoretical curves are discussed in the text. The con- 
stant C =11873 a)’. 


for, the evidence seems to point clearly to the 
configuration 3d"4s for this case. The small- 
ness of the spread in shifts between the different 
metals indicates strongly that the configuration 
of the solute Fe is substantially the same in the 
metals investigated. For the Fe** and Fe** 
compounds, the experimental data are entered 
upon the 3d°4s* and 3d°4s* curves, respectively. 
This is consistent with the idea that in these 
covalent compounds the 4s atomic orbitals are 
partially occupied by electrons from the ligand 
ions (bonding orbitals).*°** The data for the 
“ferro” and “ferri” cyanides of potassium have 
not been plotted since it is not clear how to fit 
them into the above scheme. Indeed, since they 
have ground states which do not follow Hund’s 
rule, it is unlikely that Watson’s calculations are 
applicable to them. The fact that the isomer 
shifts in both cyanides are very small relative 
to each other and to stainless steel appears to be 
fortuitous. 

From the observed isomer shift and the cal- 
culated difference in |y(0)|? for the d® and d® 
free-ion configurations we may compute the dif- 
ference of the excited and ground-state charge 
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Table I. Observed shifts in gamma-ray energy in 
various iron compounds and d-group metals measured 
relative to type 310 stainless steel. (Source and ab- 
sorber at room temperature.) The measured shifts 
contain contributions from the second-order Doppler 
shift and the zero-point energy shift in addition to the 
isomer shift discussed here. Measurements at 77°K® 
indicate that the second-order Doppler shift contribu- 
tions are small compared to the total shift. More- 
over, since the Doppler shifts are likely to be of sim- 
ilar magnitudes in the compounds considered, they 
will cancel, since we are ultimately concerned only 
with differences in the shifts between similar mater- 
ials. The same is true of the zero-point energy shift. 
The uncertainty in comparing metals with salts is 
considerably greater. To convert the shifts from 
cm/sec to Mc/sec, one should multiply the values by 
116. The errors indicated reflect the uncertainty in 
the last significant figure. 








Shift 
(cm/sec) 
3d® FeF, (single crystal) 0.140 +5 
3d® KFeF; 0.139 +5 
3d® FeSO,: 7H,O 0.140 +5 
3d® FeCl, 4H,O 0.130 +5 
3d® FeS ~0.11 +1 
3d$ Fe,(SO,)3° 6H,O 0.052 +5 
sd Fe,0;” 0.047 £5 
3a® Yttrium-iron garnet, 
octahedral 0.057 +5 
3d$ Yttrium-iron garnet, 
tetrahedral 0.026 +5 
FeS, (pyrites)° 0.048 +5 
FeS, (marcasite) 0.048 +5 
Metals Fe? 0.015 +5 
Co 0.012 +5 
Ni 0.015 +5 
Mn -0. 008 +2 
Cr -0. 005 +2 
Mo -0. 001 +2 
Cyanides K,Fe(CN),*3H,O 0.0083 +10 
K3Fe(CN)<¢ 0.0000 +10 








@R, L. Ingalls, G. Lang, and S. DeBenedetti, Bull. 
Am. Phys. Soc. 5, 429 (1960). 
e reference 5. 
“See reference 6. 
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radii using Eq. (1), assuming the usual charge 
radius dependence on mass number, i.e., 
R=1.20A"*x107"3 cm. We obtain 6R/R =1.8x10~° 
as the fractional change in the charge radius, 
with the effective radius of the ground state 
larger than that of the excited state. This result 
is not unexpected” in sign and magnitude. The 
shell model predicts that the ground state for 3 
odd neutrons outside a closed shell (28 neutrons), 
corresponding to a hole in the 2p,, shell, is one 
for which J=3/2, as is the case for the 14-kev 
excited state. The large spin-orbit coupling pre- 
cludes the possibility of a 2p,,. configuration 
being the ground state. 

Using the radial moments for an isotropic 
square well’® and even assuming a proton ex- 
citation corresponding to A/=3 or 4, one still 
obtains a value of only 2x107° for 6R/R. (The 
more likely case of a neutron excitation would 
leave the charge radius unchanged in a first ap- 
proximation.) It is most unlikely therefore that 
the ground state is a simple shell state. It is 
interesting to note that the charge radius change 
expected for the addition of one particle (isotope 
shift at A =57) is 6R/R = 4(64/A) =5.9x10. 

In the isomeric transition of Hg’” gr =1/2— 
Igy =13/2), it was estimated? that the charge re- 
distribution corresponded to an increase of one- 
fourth to one-fifth of that observed experimentally 
for the addition of one neutron at A=157. (The 
excited state has the larger charge distribution 
in this case.) Since both states are identifiable 
as single-particle neutron states, it is clear that 
a general charge redistribution accompanies the 
isomeric transition. 

We would like to thank W. E. Blumberg, A. M. 
Clogston, and M. Goldhaber for several critical 
discussions. 
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EVIDENCE FOR THE EXISTENCE OF SURFACE EXCITONS IN TETRACENE MICROCRYSTALS* 


J. Katul and A. B. Zahlan 
Physics Department, American University of Beirut, Beirut, Lebanon 
(Received January 9, 1961) 


During the past few years there have been a 
number of predictions of new types of excitons in 
molecular crystals.’ 

Selivanenko’ predicted the existence of surface 
excitons in molecular crystals with a high heat 
of fusion. Ferguson and Schneider® observed a 
b-polarized shoulder in a 0.2y anthracene film 
which was not observed in thicker films. These 
authors believed that it was unlikely that this 
band be due to a surface exciton. However, 
McClure? interpreted this 25100-cm™! band as 
the b-polarized factor group component of the 
pure electronic transition occurring in these sur- 
face molecules. 

We have observed two absorption bands in thin 
tetracene (microcrystalline) films prepared by 
rapid sublimation under conditions which produce 
a high surface-to-volume ratio. Rapid sublima- 
tion is essential” to increase the density of sur- 
face states. 

No absorption above 5200 A was observed in 
thick crystals grown from the vapor under vac- 
uum. However, when these same crystals were 
sublimed rapidly, under vacuum, onto a quartz 
disk forming a highly imperfect surface, a band 
at about 5400 A was observed. The two bands in 
the range 530-600 my could not be observed in 
the solutions made by dissolving the films in 
suitable solvents. This proves that these bands 
were not due to an impurity produced during the 
preparation. 

Films were prepared under a pressure of 30 , 
but with different rates of sublimation. No ab- 
sorption in the range 520-600 my occurred in 


slowly sublimed films. However, films prepared 
by rapid sublimation (at 130°C) exhibited one 
strong band when the film was about 0.5 » and 
two strong bands at a thickness of 1- 1.5 yu. The 
1 » need not be the thickness of the microcrys- 
tals involved. This observation of a critical 
thickness is an indication that the 5330A band (I) 
involves either a minimum crystalline length 

or a crystalline structure that can exist in stable 
form only above a certain thickness. Below such 
a thickness either the strains or the large num- 
ber of imperfections, due to the high (surface/ 
volume) ratio, make such a crystalline struc- 
ture unstable. Table I summarizes sample data 
on the positions and extinction coefficients. Band 
II at 5350 A shifts continuously to about 6000 A 
with increasing thickness. Band I is very sharp 


Table I. Dependence of the position and intensity of 
the surface excitons on thickness. 








Wavelength Thickness €x1073 Bandwidth 
(A) (4) (cm-*) (em=) 
5344 (I) 1-2 25- 75 ( 80°K) 
5650 (II) 0.4-2 100 - 300 (300°K) 
5340 (II) 0.4 4.2 
5410 (II) 1 4.4 
5540 (II) 1.3 4.9 
5680 (II) 1.5 4.4 
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FIG. 1. In (a) and (b), bands I and II are clearly seen 
at 80°K. The thickness of the film in (a) is about 1.5 4 
and in (b) it is about 0.5 4. Note the shifts of I and II 
with thickness. In (c) the continuous shift of II with 
thickness is seen clearly at 300°K. Note that I appears 
suddenly at a critical thickness (see text). The 0-0’ 
of the molecular 'A,,.~ ~'B,_* transition is shown for 
comparison. A large Hilger spectrograph with glass 
optics was used throughout. 


and shifts with thickness (see Fig. 1) about 100 A 
over the range 1-2 uw. The critical thickness 
depends on temperature. At 80°K, I is observed 
in a thinner film that at 300°K. 

Besides the shift in the position of band II with 
thickness, its intensity may increase or decrease 
as it shifts to longer wavelength. We do not re- 
port oscillator strength, and due to variations in 
bandwidth with thickness, the intensity of the 
band is not linearly proportional to the extinction 
coefficient. 

Heating in air up to 130°C for 3 hours did not 
induce any absorption in films where the bands 
are absent. Furthermore, a similar heat treat- 
ment did not change the bands when these were 
present. Tetracene is extremely sensitive to 
oxidation in solution but stable in the solid state. 
The lack of sensitivity of these bands to oxida- 
tion means that these states involve several 
“crystalline” layers of tetracene, and due to the 
packing in a crystal, the oxygen molecules can- 
not penetrate deep enough to produce observable 
damage. 


102 


Using a 10° magnification and polarized light, 
a difference between films that do not display 
“surface states” and those that do was observed: 
The nonabsorbing films consisted of microcrys- 
tals whose linear dimensions were of the order 
of 5 microns, about 1 micron thick, and appeared 
to be of the flaky type. As for the “absorbing” 
films, their dimensions were just at the limit of 
visibility and must be less than a micron. It was 
not possible to determine the shape, or dimen- 
sions, of these microcrystals. Furthermore, a 
large number of needle-like crystals were ob- 
served along with these microcrystals. Thus two 
different crystalline forms exist in these films. 
I and II might then arise in the same crystal or 
in two different types of crystalline structure. 

Two arguments in favor of surface states 
rather than surface imperfections are the shift 
of the energy with thickness and the sharpness 
of the bands. The shift in the energy is of the 
order of 10%, i.e., 0.2 ev, from the thin to the 
thick edge of a film. It has been pointed out® that 
surface imperfections could lead to new states 
but theoretical and experimental proofs do not 
exist to support this hypothesis. Surface imper- 
fections are, however, important perturbations 
of surface states. Since such imperfections act 
to decrease the degree of surface order they 
also act to raise the energy levels. This might 
account for the observed shift with thickness—it 
being assumed that film thickness determines 
approximately the size of the microcrystals. 

Since there should be a large variety of surface 
imperfections, there should be an equally broad 
distribution in the energy levels. This would 
lead to broad absorption bands, whereas we ob- 
serve a very sharp band I and band II is still 
sharper than the molecular ‘A, , -'Bo,,* bands 
(these bands have half-widths of 480, 340, 480, 
780 cm™~*) in the solid state at 90°K. 





*This work was supported by the Research Corpora- 
tion, American Chemical Society Petroleum Research 
Fund, and the American University Arts and Sciences 
Rockefeller Research Grant. 

'D. S. McClure, Solid State Physics, edited by F. 
Seitz and D. Turnbull (Academic Press, Inc., New 
York, 1959), Vol. 8, p. 1. 

2A. S. Selivanenko, J. Exptl. Theoret. Phys. 
(U.S.S.R.) 32, 75 (1957) (translation: Soviet Phys. - 
JETP 5, 79 (1957)]. 

3J. Ferguson and W. A. Schneider, J. Chem. Phys. 
28, 761 (1958). 
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CYCLOTRON AND PARAMAGNETIC RESONANCE IN STRAINED CRYSTALS 


G. E. Pikus and G. L. Bir 
Semiconductor Institute, Academy of Sciences, Leningrad, U.S.S.R. 
(Received January 11, 1961) 


Recent studies on cyclotron and paramagnetic 
resonance in strained p-Si’»? have confirmed the 
earlier developed theory,** and have shown that 
investigations of this kind give important informa- 
tion concerning band structures and the nature of 
impurity centers. In connection with this, some 
new possibilities offered by measurements of the 
resonance effect in strained crystals are pointed 
out. 

I. Measurements of the cyclotron resonance 
in strained Ge and Si make it possible to find not 
only the values of the constants A, B, and D de- 
termining the spectrum E(k) in these crystals® 
but also the ratio of the deformation potential 
constants, b/d, which is very important for the 
theory of scattering.* Thus, according to refer- 
ence 3, when the crystal is strained along [110], 
the inverse effective masses on the main axes, 
[110], [110], and [001], are equal to 


AzD 1+yB/D ,,D 1-yB/D ,. By 
2 +/)”’ 2 (1+/)”’ (1+77)”’ 
(1) 
where 
_ 2(b/d) cq, 
Cy Cy. 


Here we assume that the upper choice of signs 
corresponds to de <0 while the lower one cor- 
responds to de>0O. The signs of the constants B 
and D cannot be determined by cyclotron reso- 
nance. Such experiments as well as measure- 
ments of the piezoresistance at high temperatures 
allow us to determine only the signs of the pro- 
ducts Bb and Dd. From the data’ as well as from 
Smith’s experiments® it is evident that for Si, 
Bb>0O and Dd>0. 

Il. The use of strained crystals with an InSb- 
type lattice makes it possible to obtain semicon- 
ductors with some special kind of spectrum when 
the surfaces of constant energy have toroidal 
shape and the extremum is in the ring. Taking 
into consideration terms linear® in k at low k for 
the lower band corresponding to the minimal en- 
ergy of the holes, we obtain from the calculations 


7 


E(k) = E(k) + ain! Da, bey. (2) 
€ i,j 


Here Ek) has significance for the strained p-Ge 


corresponding to that determined by Eq. (18) of 
reference 3, K is Kane’s constant,’ 6, is de- 
termined by Eq. (16) of reference 3, and 


a =3b%e -e€ )*+4d%e *+e 7), 
%X yy 22 xy  =—«KB 
@ =-2V3(& )@de +2d*e € , etc. 

xy € xy NZ yz 


This formula is just for the strains when the band 
splitting, which is equal to 2(§,<)”, exceeds kT 
and the difference E(0)-E(k) yin, i-e., under 
conditions similar to those attained during the 
above-mentioned experiments.’ 

Toroidal surfaces occur, in particular, in the 
strains on the axes [001] and [111]: When de >0, 


E(k) =(A +$B)(k +k, +(A-B)k?; (3) 
when be <0, 

E(k)=(A -3B)(k, #k, )+(A+B)k,*; (4) 
where 


v3 IKI 
2 


2-4 2 2 - 
R k.. +, ’ k AziB° 


10 


Correspondingly, when de >0, 
E(k) -(4 2 )ee,s ko? “(4 . - k,?, (5) 
2v3. v3" 
where 
___!Kl 
a8 A +D/2NS 


On the contrary, when de <0, 


E(k) = ~+-)a: (a2 )a, +k, ” 
v3. 


where 


_ v2 |KI 
A+DN3 


The properties of semiconductors with toroidal 
energy surfaces have been studied in detail by 
Rashba et al.’° Of great interest is the investi- 
gation of the cyclotron resonance and the com- 
bined resonance predicted by Rashba in these 
semiconductors. 
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Ill. In strained crystals p-Ge and p-Si it seems 
possible to observe a spin resonance of free car- 
riers. In these nonstrained crystals such a reso- 
nance cannot be seen, because as a result of the 
large spin-orbit interaction an effective magnetic 
moment of the hole depends on its impulse direc- 
tion. Therefore, instead of a resonance line, a 
broad resonance spectrum would result. 

In strained crystals all the holes with low k 
possess the same effective magnetic moment de- 
pending on the strain direction, and their reso- 
nance frequency is determined by 


(hw P= (ug*K*/8 XH, H), (6) 


where 


(A,B) =28,,A B,, (7) 


i." [(s +b(A - 3c) + sale," + € 2» 


a os y2 _ a 
B,, =VSdfVSde¢, € [28 ) b(A se Ih etc., 


where A=€,,+€), +€22 and X is Luttinger’s 
constant." In the case of strain along [001] and 
[111], 


thw P= Hy Kg)" +8,"H |”), (8) 
where gy" =(142), g,?=(141), H,’=H,’+H,’, 
and p., is the Bohr magneton. Here the upper 
sign corresponds to be >0 and de >0 while the 
lower one corresponds to be <0 and de <0. 

The matrix element determining the probability 
of the transition under the influence of an alter- 
nating magnetic field h is equal to 


- => 


IWis |? s Ko (H, H) (h, h) me (H, h)? 


= 9 
(H, H) ™ 





€ 


where (A, B) may be determined according to (7). 
For the indicated directions [001] and [111], when 
be >0 or de>O, 


[hit]? +3[hHt] .? 
H#+3H?” 
& 





|W,217 =x" yp," (10) 


when be <0 or de<0, W,,=0. 

Significant changes in the resonance frequency 
of the hole of the acceptor center occur when the 
splitting of the ground fourfold-degenerate state” 
due to the strain exceeds the term splitting in the 
magnetic field. However, in the case of strains ob- 
tained in practice,” a main contribution to the y 
function of the impurity center comes from states 


104 





with large k for which Eq. (6) does not apply. 
That is why the resonance frequency of the bound 
hole depends not only on the constants b and d but | 
also on the 7 function.? Therefore it is not to be 
believed, as is accepted in Kleiner and Roth, 
that the term splitting is equal to the splitting of 
the valence bands at k=0, for the term splitting 
also depends on the change of the corresponding 
effective masses. 

Only in cases when the strain is so large that 
the main contribution to the y function comes 
from states of low k, i.e., when the band splitting 
2(& «)* exceeds the activation energy of the im- 
purity center £;, will the resonance frequency be 
the same as for free holes. Therefore, in the 
case of lower strains this frequency must depend 
notably on the strain value according to 

= v2 
4/%» 2( 6) /E;. (11) 
In the above-mentioned experiments* Ag, /eg, =30 
which, according to (11), gives for the constant 
b a value about 0.5 ev. 

We have no grounds to believe that in this in- 
stance the ratio g\ /g, undergoes considerable 
changes for the bound states when compared with 
the free ones. Therefore, comparing the results 
of the experiments? with Eq. (8), the conclusion 
may probably be drawn that in Si both constants 
b and d are negative; hence B and D are also 
negative, which corresponds to the theoretical 
estimate.°® 
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SELECTION RULE FOR THE INTERACTION OF MICROWAVE ULTRASONICS WITH SPINS 


N. S. Shiren and E. B. Tucker 
General Electric Research Laboratory, Schenectady, New York 
(Received December 27, 1960) 


Theory and results for iron-group paramag- 
netism in crystals’ show that the dominant con- 
tribution, of noncubic crystalline fields, to the 
ground-state spin Hamiltonian is quadratic in the 
spin. This indicates a quadrupolar form for the 
relevant spin operator in spin resonance transi- 
tions induced by lattice vibrations. Quadrupole 
selection rules for spin- phonon transitions have 
been explicitly derived by Mattuck and Strand- 
berg? for non S-state ions of the iron group. 
However, their low-frequency acoustic experi- 
ments® were unable to verify the selection rule. 
Using microwave ultrasonics in conjunction with 
standard spin resonance techniques,* we have 
experimentally confirmed the quadrupole selec- 
tion rule for a non S-state ion, Cr** in Al,O,, 
and also for two S-state ions, Mn** and Fe* in 
MgO. 

9-kMc/sec longitudinal phonons were generated 
in quartz rods bonded to the paramagnetic host 
crystals.°*® The wave propagation direction was 
parallel to the c axis in the ALO,, and parallel 
to a cubic axis in the MgO. CW resonance satu- 
ration techniques were used to observe the in- 
teraction of the ultrasonics with the spin reso- 
nances, and relative spin-phonon matrix ele- 
ments were extracted by comparison of the satu- 
ration by ultrasonics with saturation by a micro- 
wave source, at the resonance frequency, which 
was not coherent with the electron spin resonance 
spectrometer power.® Measurements of the ma- 
trix elements were made as a function of the 
angle @ between the dc magnetic field direction 
and the wave axis. 

The dominant quadrupole nature of the relevant 
spin operator is shown by comparison of the ob- 
served angular variations with those computed, 
using the quadrupole spin operator. Good agree- 
ment was obtained for all three ions. As an ex- 
ample, Fig. 1 shows the experimental and theo- 
retical angular variations of the square of the 





matrix element for the (-3/2, -1/2) transition 
of Fe** in MgO. In the absence of crystal field 
terms in the spin Hamiltonian, mgs is a good 
quantum number along H4y,, independent of the 
angle 6, and the probability for 4m, =+1 transi- 
tions, caused by a longitudinal wave, varies as 
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FIG. 1. Angular variation of the spin-phonon ma- 


trix element for the (-3/2, -1/2) transition of Fe** in 
MgO. The theoretical expression is normalized to the 
average of the experimental values. The small back- 
ground level is believed due to generation of other 
acoustic modes, either in the quartz-MgO bond, or in 
the MgO crystal. 
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cos*@ sin*6. However, for Fe** ions in MgO the 
cubic field causes a small admixing of states* 
and consequently, a slight deviation of the transi- 
tion probability from a simple cos*6 sin*6@ de- 
pendence. Therefore, exact wave functions have 
been used in computing the theoretical angular 
variation. 

A quadrupole spin operator cannot induce 
transitions between pure -1/2 and +1/2 states, 
and the near-zero ultrasonic effect observed on 
this transition for all three ions is further evi- 
dence for the quadrupole interaction. 

The authors would like to thank Dr. G. D. Wat- 


kins for helpful discussions, and Mr. T. G. 
Kazyaka for his aid in taking the data. 
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POPULATION INVERSION AND CONTINUOUS OPTICAL MASER OSCILLATION 
IN A GAS DISCHARGE CONTAINING A He-Ne MIXTURE 


A. Javan, W. R. Bennett, Jr., and D. R. Herriott 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received December 30, 1960) 


This Letter is intended to give a summary of 
the results of experimental determinations of 
several physical properties of a gaseous dis- 
charge consisting of a He-Ne mixture which have 
led to the successful operation of a continuous- 
wave maser at five different wavelengths in the 
near-infrared. Population inversions are achieved 
between several Ne levels by means of excitation 
transfer from the metastable He(2 °S) to the 2s 
levels of Ne.’»? The maser oscillation takes 
place in a narrow beam with a diameter of 0.45 
inch. The power in our strongest beam at the 
wavelength 11530 A is 15 milliwatts. The meas- 
ured linewidth is in the range of 10 to 80 kc/sec. 
The angular spread of the beam is less than one 
minute of arc. 

In the present system, the He(2 °S) metastables 
are used as carriers of energy to excite the levels 
of Ne. Due to the resonant nature of this process, 
the excitation cross section is large only to those 
levels of Ne which fall in energy within a few kT 
of that of the He(2 *S). Figure 1 gives the per- 
tinent energy levels of He and Ne. As may be 
seen from Fig. 1,° the energy separations be- 
tween the He(2 °S) and Ne(2s) states are small 
enough to insure an appreciable degree of trans- 
fer to the 2s states; however, the energy gap 
between the 2 °S and 2p levels is much too large 
to permit direct population of the 2 states in the 
collision process. Although the 2s, and 2s, levels 
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of Ne may radiate to the Ne ground state, in the 
limit of complete resonance trapping their life- 
times are determined primarily by radiative de- 
cay to the 2p levels. Under this condition, the 
lifetimes of all of the 2s levels are about one 
order of magnitude longer than those of the cor- 
responding 2p states. (The decay of the 2p levels 
is due to their radiative decay to the 1s levels.) 
Thus population inversions may be obtained on 
each of the thirty allowed 2s — 2p transitions (see 
below). 

In order to determine the transfer cross sec- 
tion and spurious sources of excitation of Ne 
levels, the afterglow of a pulsed rf discharge 
containing varying amounts of He and Ne was 
studied in detail. The time variation of the den- 
sities of the 2s and 2p levels was observed by the 
fluorescent decay of these levels using appro- 
priate isolated transitions. The time variation 
of the 2 °S density was determined through the 
absorption of the 2 *P +2 °S He transition. 

The conditions were found under which the 
dominant mode of excitation of Ne levels arose 
from the He(2 °S) transfer. Among other sources 
of excitation of the Ne levels it was found that 
inelastic collisions of energetic electrons with 
Ne metastables? (i.e., the 1s levels*) were by 
far the most important. It should be noted that 
certain amounts of 2p excitation, arising from 
sources other than 2s -2p cascade, can be tol- 
erated without upsetting the presence of a popula- 
tion inversion. For instance, let us consider 
the population inversion between the 2s, and 2p, 
levels. The 2p, is optically connected to only one 
of the 2s levels, namely the 2s,. The Einstein A 
coefficient for the 2s, ~- 2p, transition was meas- 
ured and found to be at least 25 times smaller 
than the rate of decay of the 2p, level to within 
an accuracy of 20%. (This measurement was 
done by means of two independent techniques. 
The first was a spectroscopic technique to be 
described in a later publication. The second is 
discussed below.) This means that even an order 
of magnitude added excitation of 2p, over that due 
to the cascade from the 2s, still yields a popula- 
tion of the 2s, levels larger than twice the popula- 
tion of the 2p,. In view of this effect, the ad- 
justment of the electron density was found to be 
most critical. 

In the afterglow, the energetic electrons ther- 
malize within 10 to 20 usec. The excitation 
transfer from He(2 *S), however, continues until 
the time when all of the 2 °S states are quenched. 
This effect gives rise to a decay of the 2s levels 


identical to the time variation of the He(2 °S). 
The afterglow of the 2p levels also shows an iden- 
tical decay due to the cascade transition from the 
2s levels. Our observations indicate that the 
time variations of the 2 °S, 2s, 2p levels in the 
afterglow are indeed identical to a high degree 

of accuracy. The possibility of the direct trans- 
fer from the 2 °S to the 2p levels was ruled out 
by a process of relative intensity measurement 
at a fixed time in the afterglow to be described 
elsewhere. 

Study of the fluorescent decay rate of the 2s 
and 2p levels as a function of several Ne pres- 
sures gave a total velocity-averaged inelastic 
cross section of (co) =(3.7+0.5)x107!" cm? for 
destroying He(2 °*S) metastables in two-body col- 
lisions with Ne ground-state atoms at 300°K. 
This cross section can safely be assigned to the 
excitation transfer to the 2s levels of Ne since 
we have been able to show that other sources of 
destruction of the He(2 °S) by Ne ground-state 
atoms such as the one leading® to a formation of 
HeNe* are negligible. Our cross-section meas- 
urement at the same time furnishes the value of 
the diffusion coefficient of 2 °*S He metastables. 
Our result is in excellent agreement with the 
value of (470+ 25 cm?/sec) determined by Phelps.® 

The radiative lifetimes of the 2s and 2p levels 
were determined using techniques described 
previously.” Our tentative values of the lifetimes 
of the 2s,, 2s,, 2s,, and 2s, fall in the range 
from about 100 to 200 mysec. Our results for 
the 2p levels are in general agreement with 
Ladenberg’s® results. These measurements are 
still in progress and will be reported more fully 
at a later date. Each of the 2s levels decays into 
several of the ten 2p levels. The measured and 
estimated branching ratios for the transitions 
originating from each of the 2s levels, together 
with their total rate of radiative decay obtained 
from the above lifetimes, give the values of the 
individual Einstein A coefficients connecting the 
various 2s and 2p levels. 

From the early stages of this experiment we 
were able to observe single-pass amplification 
of a transmitted signal at an infrared frequency 
where the inverted population was expected. 
Several electronic detection schemes were de- 
vised for this purpose. One of these schemes 
responded only during a small portion of the 
time and was in synchronism with the rf pulsed 
discharge. This technique enabled us to observe 
a larger value of amplification in the afterglow 
than during the glow; the latter is to be expected 
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since the added excitation of the 2p levels by 
electron impact is absent in the afterglow. How- 
ever, in addition to extreme difficulties encount- 
ered in the interpretation of the results, the gain 
measurement alone was not sufficient to shed 
insight into the nature of the physical processes 
involved in the discharge. 

In view of the above considerations, we were 
able to determine an optimum condition for the 
maser action prior to the observation of oscilla- 
tion as described below. 

The continuous-wave maser oscillation is ob- 
served in a discharge containing 0.1 mm Hg of 
Ne and 1 mm Hg of He. The optical feedback® 
is obtained through high-reflectance Fabry-Perot 
plates placed within the gas chamber. This cham- 
ber consists of a long quartz tube with an inside 
diameter of 1.5 cm and length of 80cm. The 
quartz tube is terminated at each end with a lar- 
ger metal chamber containing the high-reflectance 
plates. Flexible bellows are used in the end 
chambers to allow external mechanical tuning of 
the Fabry-Perot plates. Two optically flat win- 
dows at the extremes of this system are pro- 
vided to allow the maser beam to leave the in- 
strument undistorted. The separation of the 
plates is one meter. The discharge is excited 
by means of external electrodes using a 28- 
Mc/sec generator. The estimated power dis- 
sipation is around 50 watts. 

The high reflectance is achieved by means of 
13-layer evaporated dielectric films on the hun- 
dredth wavelength flat surface of the fused silica 
plates. The reflectance is 98.9% with 0.3% 
transmission in the wavelength range 11000 A 
to 12000 A. Outside of this range of wavelength 
the reflectance decreases rapidly. 

The initial adjustment of the plates prior to 
starting of the discharge is done using an auto- 
collimator. We have observed oscillations at the 
following wavelengths: 11180, 11530, 11600, 
11990, and 12070 A. The strongest oscillation 
occurs at 11530 A with an output power of 15 mw. 

An image converter tube is used visually and 
with a camera to observe and record structure 
of the beam. Two types of beam patterns have 
been studied. 

A cross section of the emergent beam can be 
intercepted by the photocathode of the image con- 
verter and the intensity distribution recorded. 
This near-field pattern is shown in the two upper 
photographs of Fig. 2 for two Fabry-Perot plate 
angles. 

The intensity pattern that the beam would form 
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PLATES PARALLEL 6 SECOND ANGLE 


PATTERN IN MASER CAVITY 





PLATES PARALLEL 6 SECOND ANGLE 


ANGULAR DISTRIBUTION IN BEAM 


FIG. 2. Far-field and near-field patterns from 
maser beam at two plate angles. 


at an infinite distance can be observed by imaging 
the maser beam to a point and observing the dis- 
tribution of light in this image with an image con- 
verter microscope. The intensity distribution 
in this image represents the angular distribution 
in the maser beam. The two lower pictures in 
Fig. 2 show this far-field pattern at the indicated 
plate positions. The structures seen in the near- 
and far-field patterns would indicate the existence 
of several modes at a finite plate angle. The 
patterns shown in Fig. 2 were made at constant 
exposure and are therefore greatly overexposed 
by the increased intensity with parallel plates. 
The half-intensity width of the parallel plate out- 
put is less than one minute of arc. 

The phase variation over the area of the beam 
was studied by placing a double or multiple slit 
in the beam and recording the diffraction pattern. 
A nine-slit mask of 0.050-in. spacing and 0.002- 
in. slit width covering the full 0.45-in. aperture 
of the beam results in a diffraction pattern shown 
in Fig. 3. The minima are very low, indicating 
small phase variation over the aperture. A 
double-slit pattern at 0.0175-in. is also shown. 

Since the linewidth of the maser oscillation is 
many orders of magnitude narrower than the res- 
olution of the best available spectrometers or 
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9 SLITS 0.050” SPACING 0.002” WIDTH 





0.0175” SPACING 


2 SLITS 0.0014” WIDTH 


FIG. 3. Far-field diffraction patterns obtained by 
inserting multiple slits in the maser beam. 


interferometers, standard optical techniques 
could not be used to measure the linewidth. Of 
the several possible methods examined, the most 
practical was found to be one based on analysis 
of the Fourier spectrum observed through a 
photomultiplier tube of the maser output. The 
photomultiplier is essentially a perfect square- 
law detector and should respond to beat frequen- 
cies. Hence, one can observe the spectral dis- 
tribution obtained by beating all of the components 
of the line. The exact width of the original line 
determined from this distribution is independent 
of detailed assumptions on the line shape to with- 
in a factor of about two, providing the original 
line is not multiply peaked. In practice, this 
spectral distribution will be superimposed on top 
of the flat spectrum obtained from shot noise in 
the photomultiplier tube, and the latter may be 
used to check the bandwidth of the system. 

The system was first checked by observing the 
spectrum from an intense line emitted in an 
ordinary discharge tube using a spectrum analy- 
zer. The latter was flat within the frequency 
range of our analyzer, covering from 100 cps to 
13 Mc/sec. Using a grating spectrometer to iso- 
late the strong line at 11530 A, the maser output 
was next examined. Our initial observations in- 
dicated a strong signal starting at zero frequency 
and extending to a half-width of about 100 kc/sec. 
Further examination showed the presence of an 
additional series of anharmonically related beats 
centered about frequencies other than zero. The 
location of these beats depended critically on the 
Fabry-Perot plate alignment. As the angular 
alignment of the plates was changed from 0 to 3 
seconds, the beats became more widely spaced 


and weaker in intensity. Several possible spurious 
sources of beats were carefully examined and 
found not to be responsible for the observed sig- 
nal. The possibility of “blocking oscillation” of 
the maser giving rise to an amplitude modulation 
of the output signal can be ruled out by the follow- 
ing observation: The peak frequency of a beat 
note is only a function of the plate adjustment 

and is independent of the factor by which the 
maser gain exceeds the oscillation threshold. The 
most plausible explanation found was that the 
beats arose from simultaneous oscillation in sev- 
eral different modes separated by frequencies in 
the Mc/sec range. Such a conclusion is not in 
disagreement with detailed calculations of the 
oscillatory modes by Fox and Li.’° A distinct 
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FIG. 4. Frequency spectra of the beat notes be- 


tween various maser beam components. 
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correlation was found between the visual struc- 
ture of the beam observed through the infrared 
image converter, and the presence of the beats 
as the angular alignment of the plates was varied. 
Figure 4(a) shows several beat notes ranging 
from 50 kc/sec to 2 Mc/sec. Figures 4(b) and 
4(c) show the structure of these beats on an ex- 
panded scale. It should be noted that in our pres- 
ent system the angular adjustment of the plates 

is also accompanied by a slight change in plate 
separation. The important conclusion to be drawn 
here regards the linewidth of the maser oscilla- 
tion. As is evident from Figs. 4(b) and 4(c), the 
beat notes imply linewidths in the range from 10 
to 80 kc/sec. 

The use of excitation transfer for production of 
inverted population was suggested’»” by one of 
the authors. The He-Ne mixture described above 
is the first gaseous system which has led to 
maser oscillations at optical frequencies. Re- 
cently some evidence for the presence of inverted 
populations in a Hg-Zn mixture has been reported” 
which supports the wide applicability of this prin- 
ciple.*»? 

In the course of this research, a number of 
co-workers have contributed in the design, con- 
struction, and measurement. We are indebted to 
F. Muller, P. S. Kubik, A. R. Strnad, D. MacNair, 


E. Koch, G. J. Wolff, and G. E. Reitter. Dr. A. F. 
Turner of Bausch& Lomb has been very helpful 

in supplying the special evaporated reflectance 
films. In particular we would like to acknowledge 
the expert and patient assistance given by Mr. 

E. A. Ballik in the later stages of this experi- 
ment. 
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MASS SPECTROGRAPHIC IDENTIFICATION OF THE ION OBSERVED IN 
HYDROGEN MOBILITY EXPERIMENTS* 


W. S. Barnes, D. W. Martin, and E. W. McDaniel 
Engineering Experiment Station and School of Physics, 
Georgia Institute of Technology, Atlanta, Georgia 
(Received January 9, 1961) 


In a Letter appearing in a recent issue of this 
journal, Varney’ pointed out that previous studies 
of the mobility of ions in pure hydrogen have pro- 
bably been misinterpreted. In the majority of 
experiments performed to date, only one ion has 
been observed, and it has commonly been thought 
to be H,*. Indications of a second ion, generally 
assumed to be H*, have been seen in a few ex- 
periments but the actual existence of this ion re- 
mains in doubt. Upon consideration of a study 
by Stevenson and Schissler? of the cross section 
for formation of H,* by the reaction H,* + H,- 
H,* +H, Varney has made the following asser- 
tions: (1) The H,* ion is never observed in hydro- 
gen mobility experiments; (2) the ion which in- 
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variably appears in these experiments is Hg,” 
(3) it is reasonable to assume that the second 
ion, if actually present, is H*. The purpose of 
this Letter is to report evidence which tends to 
corroborate Varney in each of these assertions. 
We have recently constructed and put into 

operation apparatus designed to permit identifi- 
cation of the ions appearing in mobility experi- 
ments and to yield information concerning the 
various types of reactions which can occur at 
low energy between ions and molecules.* Ions 
are produced by electron bombardment inside a 
long drift tube containing gas at a pressure up 
to about 700 microns. The ions diffuse down the 
drift tube under the influence of a weak electric 
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field. A sample of the ion population at the end 
of the drift tube is extracted through a two-stage 
differential pumping chamber and passed into a 
60° magnetic deflection mass spectrometer. The 
number of ion- molecule collisions in the drift 
tube may be varied over a wide range by changing 
the source position and/or the gas pressure. In- 
formation concerning the nature and probability 
of the reactions occurring is revealed by the re- 
sulting changes in the ionic mass spectrum. 

In studies of hydrogen, the three known species 
(H*, H,*, and H,") are all observed at very low 
drift-tube pressures, p, and for short drift dis- 
tances, d. The H,* ion is by far the most abun- 
dant throughout the range covered. As pd is in- 
creased, the abundances of H* and H,” relative 
to that of H,* diminish, the H,* abundance falling 
off more rapidly than that of H*. Taking into ac- 
count the increasing ion losses due to diffusion, 
recombination, etc., it appears that the H,* po- 
pulation increases mainly at the expense of H,*. 
This is consistent with Stevenson and Schissler’s 
finding that the reaction H,* + H,~H,* +H has a 
large cross section. The persistence of H* ions 
in our spectrum to higher pd values than that at 
which H,* disappears indicates that the cross 
section for the reaction H*+H,~-H,* is much 
smaller, consistent with the conclusions of other 


investigators. 

In mobility experiments, pd is always made 
rather large in order that a steady state, con- 
stant drift velocity be achieved. Values of pd of 
10 cm-mm Hg, and higher, are commonly used. 
At pd 21.5, only the ky ion is seen in our spec- 
trometer. This fact appears to corroborate 
Varney’s first and second assertions. Further- 
more, since the H* ion persists to much higher 
pd than does the H,* ion, it seems reasonable to 
assume that if a second ion is observed, it is H*. 
Of course, the possibility exists that the second 
ion sometimes reported is an impurity. 

A more complete discussion of the work re- 
ported here will be presented in another paper 
which will also contain a detailed description of 
the apparatus. 





*This research was supported by the U. S. Air 
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COLLISIONAL DETACHMENT IN MOLECULAR OXYGEN* 


A. V. Phelps and J. L. Pack 
Westinghouse Research Laboratories, Pittsburgh, Pennsylvania 
(Received December 30, 1960) 


The measurements reported in this Letter show 
that the electron affinity of the negative ion of 
molecular oxygen formed in low-energy swarm 
experiments is 0.46+0.02 ev. This value is 
three times that used in the analyses of recent 
photodetachment experiments.’ In addition, we 
show that the frequency of electron detachment 
collisions between oxygen negative ions anda 
thermal distribution of oxygen or nitrogen mole- 
cules in the lower ionosphere is at least two 
orders of magnitude smaller than the accepted 
value.” Therefore, a new mechanism for elec- 
tron detachment must be found in order to ex- 
plain the appreciable electron densities found in 
the lower ionosphere at night. 

The experimental apparatus used in these 


measurements is the same as that used by Chanin, 
Phelps, and Biondi® to measure electron attach- 
ment coefficients. In the present experiments, 
either (a) a pulse of voltage is applied to the grid 
to collect negative ions and electrons in the vi- 
cinity of the grid and thus reduce the current 
reaching the collector in proportion to the sum 

of electron and negative-ion currents reaching 

the grid, or (b) a pulse of high-frequency voltage 
(~1 Mc/sec) is used to collect only electrons* 

and reduce the collector current in proportion to 
the electron current reaching the grid. The wave- 
forms obtained are shown in Fig. 1. The upper 
waveform for T =423°K shows an initial peak due 
to electrons which have crossed the tube without 
attachment. This peak is followed by an exponen- 
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FIG. 1. Change in collector current as a function of 
time of application of voltage pulse to control grid at 
various temperatures for N =10'* molecules/cc, E/N 
=3.1x 107!" volt cm’, and a drift distance of 2.5 cm. 
The durations of the light pulse and the grid voltage 
pulses were 120 usec. 


tially increasing current due to negative ions. 
The lower waveform shows the initial electron 
peak and a delayed component due to electrons 
which are detached from the ions as they drift 
across the tube. At higher temperatures the ion- 
current waveform is distorted and the relative 
magnitude of the delayed electron current is in- 
creased. These waveforms agree well with theo- 
retical predictions. For a sufficiently narrow 
pulse width (At), the ratio of the electron current 
immediately after the initial pulse (Jz) to the 
peak value of the initial electron-current pulse 
(Ig) is given by Ig/Ig =vatevgAt, where vg and 

vq are the frequencies of electron attachment 
and detachment collisions and ¢g is the electron 
transit time. Under the conditions of these ex- 
periments, vt, is obtained from the ratio of the 
total current (I) to the current of electrons 
reaching the grid without attaching to the oxygen 
(les). Thus, vgte =In(I;/Ies). Igs is obtained by 
applying a large enough grid voltage pulse at the 
time of the initial peak to collect all of the elec- 
trons. The resultant values of vg are directly 
proportional to the oxygen density, N. At the 
value of E/p used, the negative ions are in 
thermal equilibrium with the gas. 

Figure 2 shows the measured values of vg/N 
as a function of the oxygen temperature. Also 
shown are the results of measurements of the 
three-body attachment coefficient, v,/N’*, for 
thermal electrons in oxygen.* Similar values 
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FIG. 2. Thermal attachment and detachment coeffi- 


cients in molecular oxygen as a function of gas tempera- 


ture. 


have been obtained by van Lint, Wikner, and 
Trueblood.® The results presented are consis- 
tent with the assumption that we are studying the 
reaction e +20,=-0,~+0O,. By equating rates of 
attachment and detachment and using the law of 
mass action, one finds that 


va = N(O, ) 
v ,N(O,) NN(Oz) 


_ g(O.7) ae ws eAE 
~ g8(02) (os — (27). 


where the N’s and g’s are, respectively, the 
densities and statistical weights of the negative 
ions, the electrons, and the oxygen molecules 
and AE is the electron affinity. Meyerott, Land- 
shoff, and Magee® state that 2(0,~)/geg(O,) =2/3. 
The data of Fig. 2 show that AE =0.46+ 0.02 ev. 
This value is within the limits set by Loeb* and 
intermediate between the previously proposed 
values.” Further experiments will be necessary 
to determine whether the O,~ ions are formed in 
the ground vibrational state or in a vibrationally 
excited state. There is no evidence of vibrational 
relaxation due to the approximately 5 x 10° ion- 
molecule collisions occurring during the exper- 
iment. 

If our electron affinity is combined with the 
thermal attachment coefficients in Fig. 2, the 
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calculated collisional detachment coefficient at 
230°K is 4x10~*°cm*/sec compared to values of 
107?” cm*/sec or greater used in the analysis of 
ionospheric data.” Preliminary measurements 

in dry air show that, as expected, nitrogen has 
very little effect on the collisional detachment 
rate. Therefore, if the negative ions formed in 
the ionosphere are O,~, then the detachment ob- 
served at night must be due to an excited species, 
such as vibrationally excited O, resulting from 
the recombination of oxygen molecules, elec- 
tronic excitation, etc. In addition, if the O,~ 

ions in the photodetachment experiments of Burch, 
Smith, and Branscomb’ are assumed to be in the 
same vibrational state as those described above, 
then the radiative attachment coefficient for 230°K 
electrons in molecular oxygen will be increased 
from about 3 x 107" cm*/sec to as much as 
2x107* cm*/sec. 


*This work was supported in part by the Air Research 
and Development Command. 
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NEW ISOTOPE OF CARBON: C’® 


S. Hinds and R. Middleton 
Atomic Weapons Research Establishment, Aldermaston, Berksuire, England 


and 


A. E. Litherland* and D. J. Pullen 
Clarendon Laboratory, Oxford, England 
(Received December 30, 1960) 


Estimates of the neutron binding energies of a 
number of neutron-rich light nuclei have been 
made recently by Zel’dovich.’ In particular, 
from these estimates the C"*(t,p)C'® reaction was 
expected to have a Q value lying between -3 and 
-4 Mev. This reaction has now been observed 
and the measured Q value is -3.014+ 0.016 Mev. 
In addition, the decay of C** has been studied 
and the measured half-life is 0.744 0.03 second. 

A 6-Mev triton beam from the Aldermaston 
electrostatic generator was used to bombard a 
thin C™ target and the emitted protons were ana- 
lyzed with a broad-range magnetic spectrograph.” 
The targets were prepared by polymerizing acet- 
ylene gas containing 33% C™ directly onto nickel 
backing foils. Several exposures were made 
with the spectrograph using different magnetic 
field strengths and at angles ranging from 15° to 
65°. On five of these an intense proton group 
was observed which could not be accounted for 
by the known target impurities. Further analysis 
established that this group varied in energy char- 






acteristically of a mass-14 target nucleus and 
probably corresponded to the formation of C’* in 
its ground state. The kinematics of the reaction, 
however, do not allow the C**(t,p)C** reaction to 
be distinguished from either the C'*(He*, p)N** or 
the N’*(t, p)N** reactions. The He* contamination 
of the ion beam was measured to be less than 
0.1% by observing the elastically scattered par- 
ticles and was therefore unlikely to give rise to 
significant reaction groups. It was also noted 
that, although several intense groups were ob- 
served arising from the C**(t,p)C™ reaction, 
none could be attributed to the C'*(He*, p)N™* re- 
action. To eliminate the possibility of the 
N'*(¢,p)N’® reaction, two exposures were made 
with a target containing nitrogen. No group was 
observed with the same energy as the group as- 
cribed to C**. Also no proton groups were ob- 
served corresponding to the known lower excited 
states of N’* in the spectra obtained with the C* 
target, and only a weak continuum of protons was 
observed from a typical nickel target backing. 
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The energy of the incident triton beam was de- 
termined in terms of the Q value of the C!*(t, p)C** 
reaction’ and from this the Q value of the C(t, 
p)C** reaction was calculated to be -3.014+ 0.016 
Mev. This is consistent with a C’* mass excess 
of 13.694+ 0.017 Mev and a mass of 16.014702 
+ 0.000017 mass units, both referred to the scale 
on which the mass excess of C™ is zero. 

An energy level diagram showing the possible 
modes of decay of C’* is shown in Fig. 1.*~® As 
a change of parity is required for the high-energy 
beta transitions, the beta decay of C** was ex- 
pected to leave N** in neutron-unstable 1+ states. 
This expected delayed-neutron emission has been 
observed experimentally and used to determine 
the half-life of C*®. 

A one-microampere triton beam from the elec- 
trostatic generator was interrupted by a me- 
chanical shutter some 20 feet away from the neu- 
tron detector which was a block of polystyrene, 
45 cm long by 45 cm in diameter, containing five 
BF, counters. A hole of 10-cm diameter through 
the polystyrene block allowed the C* nuclei to be 
formed at the center of the assembly. The effi- 
ciency of the detector was between 1 and 2%. 

The pulses from the neutron detector and asso- 
ciated electronic equipment were fed into a cir- 
cuit which converted their time of arrival into a 
pulse height. The time distribution of the pulses 
could then be displayed on a 100-channel pulse- 
amplitude analyzer immediately after the inter- 
ruption of the beam. The time scale of the ana- 
lyzer was calibrated using the 50-cps mains fre- 
quency. 
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FIG. 1. An energy level diagram showing the pos- 
sible modes of decay of C**. 
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FIG. 2. A logarithmic plot of the decay of the de- 
layed-neutron intensity from C’*. Also shown are the 
decay curves for the other two delayed-neutron emit- 
ters Li® and N!", 


Targets containing Li’, C, and O'* were 
bombarded by the 6- Mev triton beam for periods 
between 2 and 10 seconds. The delayed neutrons 
from, respectively, Li®, O'*, and N”’ were dis- 
played in the time-analyzer for periods ranging 
from 2 to 25 seconds. A fraction of a second de- 
lay between the end of each beam pulse and the 


start of the counting period was provided manually. 


No delayed neutrons were observed from the tar- 
get backings though there was a very small back- 
ground which was thought to arise from the triton 
beam striking the mechanical shutter. Figure 2 
shows the decay of the delayed neutrons following 
the Li’(t,p)Li®, the C*(t,p)C’*®, and the O'*(t, a)N” 
reactions. The half-lives of Li® and N’” were 
found to be 0.17+ 0.01 and 4.204 0.08 seconds, re- 
spectively, and these are in good agreement with 
the published values.* The half-life of C’* was 
found to be 0.74+ 0.03 second. 

The measured mass of C** implies that the 
first T=2 state in O** is at approximately 23.0 
Mev. In the future it would be interesting to de- 
termine the neutron-unstable states of N** to 
which C’* decays and to determine the beta-ray 
branching ratio to the low-lying negative-parity 
states of N**. 
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MATRIX ELEMENTS IN THE FORBIDDEN BETA DECAY OF Ce’ 


D. D. Hoppes 
National Bureau of Standards, Washington, D. C. and Catholic University of America, Washington, D. C. 


and 


E. Ambler, R. W. Hayward, and R. S. Kaeser 
National Bureau of Standards, Washington, D. C. 
(Received January 6, 1961) 


We have performed measurements on oriented 
Ce’* in which the energy and angular distribution 
of the beta particles are determined with respect 
to the nuclear spin orientation. These measure- 
ments give information on the relative contribu- 
tion of the various nuclear matrix elements con- 
tributing to the beta transitions. 

The 33-day Ce’ undergoes beta decay to Pr’ 
via two beta groups according to the schemes 


and 


The 70% abundant £, transition with a maximum 
energy of 435 kev and a AJ=0 will depend on all 
six first forbidden (real) matrix elements: 
m(io-r), IM (ys), mir), m(oxr), mM (a), and 31L(4Bj 5). 
The 30% abundant ground-state £, transition with 
a maximum energy of 580 kev anda AJ=1 will 
depend on only the last four matrix elements, 

e., those of tensor rank one and two. 

The energy and angular distribution function, 
which has been calculated by Morita and Morita’ 
and by Bincer,” has the form 


i 


W( p-d) =No(E)+N,(E) Ef,P\( dD) +N,(B) 2 f,P,( p-J) 


+N (6) fPA( 5-9), 


where N;(E) gives the energy dependence of the 






term of order k and contains products of the re- 
duced matrix elements. The orientation param- 
eter fp describes the nuclear orientation, given 
in terms of an average over the population of the 
magnetic sublevels. The electron momentum and 
energy are p and E, respectively, and Pp( pe -J) 

is a Legendre polynomial. N,(£) gives the usual 
energy distribution for an unoriented beta emitter 
apart from the statistical factors. 

The experimental apparatus and the procedure 
for the reduction of the data are similar to those 
used in earlier studies on the angular distribution 
of beta particles from oriented nuclei.* A single 
crystal of neodymium ethyl sulfate used for both 
cooling and orienting purposes is mounted with 
the c axis vertical in the demagnetization appara- 
tus. On the uppermost surface of this crystal 
is grown a thin surface layer containing the Ce’ 
activity. Located about one centimeter above this 
is a thin anthracene scintillator which provides 
the input to a 100-channel pulse-height analyzer 
for the analysis of the beta spectra. Three 2-in. 
by 2-in. NaI gamma counters are located equa- 
torially about 15 cm from the axis of the appara- 
tus. These counters monitor the gamma-ray 
anisotropy which, combined with nuclear hyper- 
fine splitting data for Ce’*? neodymium ethyl 
sulfate, allows one to determine the nuclear 
orientation parameters fp. The outputs of the 
gamma-ray counters are also compared for time 
coincidence with the output of the beta channel by 
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means of a fast-slow coincidence circuit in order 
to separate the distribution of the 435-kev inner 
beta group from the 580-kev ground-state transi- 
tion. 

The quantities N,(E)/N,(E) and N,(E)/N,(£) can 
be measured independently. The measurements 
of the ratio N,(E)/N,(Z) were made by determin- 
ing the beta distribution when the sample has 
been aligned by the Bleaney method after cooling 
the sample by adiabatic demagnetization from a 
horizontal magnetic field of 23 kilogauss at 1.0°K. 
After a counting period of about thirty minutes 
the temperature of the sample was raised to 1°K 
again, destroying the nuclear alignment. The 
normalization for the gamma-ray anisotropy de- 
pended on a reliable measurement of a uniform 
gamma-ray distribution. In the case of nuclear 
alignment /, and f, are always zero. 

The ratio N,(E)/N,(E) was determined under 
similar conditions except that here a magnetic 
polarizing field of 200 gauss was applied along 
the c axis of the source-crystal so that the Ce’ 
nuclei were polarized by the Rose-Gorter method. 
Since this polarizing field was applied adiabatic- 
ally, the sample temperature was not as low as 
in the case when no field was applied. By taking 
different combinations of the ratio of beta counts 
before and after orientation for two opposite 
field directions one can obtain values for both 
N,/N, and N,/N,. The accuracy of the latter is 
poorer than in the alignment experiments where 
a larger f, is achieved. Counting was continued 
for alternating polarizing field directions for 
periods up to one hour. At the end of this period 
the sample was warmed and the orientation 
destroyed. The contribution of the f, term is 
calculated to be negligible at the temperatures 
achieved. 

The observed distributions were corrected for 
backscattering effects in the source and for the 
finite energy and angular resolution of the appara- 
tus; the limits on our experimental values include 
estimates of errors in these corrections as well 
as those due to counting statistics. 

The experimentally determined distribution 
function for the ground-state 580-kev transition 
is 


W( p-J)/N,=1-(1.114 0.10) F fiPr(b-d) 


2 i 
+(0.364 0.10) 2 f,P,( ped), 


where we have evaluated the coefficients for an 
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average £ =1.85 in relativistic units. A rough 
measurement of the energy distribution N,(E)/ 
N,(E) is also obtained. 

The interpretation of these results requires 
that accurate radial Coulomb wave functions for 
the electron be employed. We have used exact 
radial wave functions taking into account nuclear 
finite-size effects. These functions have been 
calculated in a similar manner to those calcu- 
lated by Bhalla and Rose.* 

Each Np(E) contains products and squares of 
the matrix elements, each multiplied with energy- 
dependent factors. In principle one can get three 
independent nonlinear simultaneous equations for 
any individual N,(E£) at three energies or for all 
three N,(E) at one energy. However, in the first 
case the energy dependence of N;(£) on the matrix 
elements is so insensitive that in practice only a 
rough limitation on the matrix elements can be 
determined. There is no evidence for a deviation 
from a statistical shape of N,(Z) as measured by 
Freedman and Engelkemier® and Jones and Jen- 
sen,° and of N,(£) within the accuracy of our ex- 
periments. 

For the purposes of illustration only, we em- 
ploy an approximation where the radial functions 
are expanded in powers of §=a@Z/2R and retain 
only those terms involving the highest power of 
€ multiplying each matrix element product in 
each particle parameter (see reference 2). This 
approximation is useful only when the matrix 
element ratios are such that energy-independent 
factors dominate, giving the statistical shapes 
to the spectra. When we do this we obtain rather 
simple expressions for the N,(E), which are 
evaluated at the average energy £ =1.85 and com- 
pared with experiment in Figs. 1 and 2. We find 
that the allowable values of certain matrix ele- 
ment combinations are greatly restricted by ex- 
periment. In order to fix the relative values of 
the individual matrix elements, we must resort 
to a model to fix at least one ratio theoretically. 

Fortunately the ground states of both ,,.Ce,,’"* 
and ,.Pr,,'** have configurations that should be 
adequately described by shell-model wave func- 
tions. For Ce™ the measured spin and magnetic 
moment establish the ground state of the eighty- 
third neutron as an (f,,.)' configuration. For Pr’ 
the measured spin, the magnetic moment of +3.92 
nuclear magnetons, and the quadrupole moment 
of -0.054 barn are consistent with a mixture’ of 
configurations (g1.)°(d,.)' and (gy2)°(dya)° for the 
protons outside the closed shell at 50. With these 
configurations the ratios of certain nuclear matrix 
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elements are predicted uniquely and are, further- 
more, completely independent of any details of 
the nuclear radial wave functions. In addition, 
admixtures of other configurations such as 
(gx2)*(dy2)°, with appropriate initial state proton 
configurations, do not alter these ratios. Admix- 
tures of higher seniority states involving the 
rearrangement of the relevant nucleons in the g,, 
and d,, shells will not contribute at all to these 
matrix element ratios. 

We have 


am(ér)/aio xP) = -1; mB, V/s xP) = -3(2). 


The other ratio 31(@)/N(o xr) involves the rela- 
tivistic operator @ which normally connects 

large with small components of relativistic four- 
component wave functions. Since most nuclear 
models involve only nonrelativistic wave functions, 
we must use a nonrelativistic operator as ob- 
tained from a canonical transformation® yielding 

a ~iv/M. This ratio is calculated to be 


4 
mu (a) _ A he’ rar (+5) “sed 


mM(oXF) a f ru, (r)* u (rar 





where the value depends on the details of the 
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nuclear radial wave functions. 

The radial wave functions are not predicted by 
the shell model without the introduction of a 
specific potential, so that we treat s(a@)/m(o xr) 
as a parameter to be determined by the experi- 
ment after the theoretical values of .(ir)/(o xr) 
and M(4B 5) /m(o xr) have been introduced. In the 
evaluation of this parameter, we use the exact 
lepton radial wave functions referred to above. 

In Fig. 3 we show the theoretical predictions 

for N,/N, and N,/N, as functions of [Cysn(a) ]/ 
[Cam(oxr)]. In this figure we have assumed 
nuclear radii R =1.2A”x107"* cm and R=1.3A™“ 
x10™** cm, and that C4 =-1.21Cy. The value of 
om(a) /: m(o xr) allowed by experiment is also given 
in Table I for the two values of the nuclear radius. 

For comparison we have also included the con- 
ventional approximation’’* for the radial functions 


Table I. Experimental values of 9n (@)/sn(@*f) for 
several evaluations of the lepton radial wave functions. 











Lepton radial 
wave functions employed 
“Exact,” R=1.34x107" cm -29.441.4 
“Exact,” R=1.24x10" cm -29.5+2.3 
“Approx.,” R=1.3A“3 x10 cm -36.6+1.0 
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[CaM  By;)]/IC 4G x¥)] have been fixed by adopting 
the single-particle values and C4=-1.21Cy. The solid 
line corresponds to R = 1.3. A”*x107"* om, and the 
dashed line corresponds to R = 1.2. A”*x10-" cm. 


expanded in powers of é for a point nucleus but 
evaluated at R =1.3A”*x107* cm. It is instruc- 
tive to see how poor the approximation is com- 
pared to the more exact treatment. 

Figure 3 indicates that for certain choices of 
the radius 90(a)/s0(o xr) may be double valued. 
The lower value can be rejected as being incom- 
patible with the observed spectrum shape. It is 
interesting to note that the large coefficient for 


the P, term is consistent with the measured shape. 


The errors quoted in Table I do not include any 
uncertainty in the correct ratio for C4/Cy in 
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this particular transition. It is also interesting 
to note that a radius much smaller than R=1.2A“ 
x107** cm would be incompatible with the experi- 
mental results. 

Several theoretical estimates of om (a@) /mir) 
based on various models have been discussed by 
Rose and Osborn,® Ahrens and Feenberg,° and 
Pursey,’° and our value for 90(@)/s0(oxr) lies at 
about the negative mean of these estimates. 

The absolute magnitude of the nuclear matrix 
elements may be fixed by fitting the calculated 
total transition probability to the observed ft 
value, viz., logft=7.7. Information on the radial 
integrals and configurational mixing is thereby 
obtained. The observed ft is several orders of 
magnitude greater than the lower limit that is 
predicted using our matrix-element values when 
one assumes a pure single-particle configuration 
with complete overlap of the radial functions. 

Finally, we have also a preliminary determi- 
nation of the 8-correlation function calculated 
by Morita and Morita" for the 435-kev, AJ=0, 
transition to the excited state. However, the six 
contributing matrix elements cannot be uniquely 
defined by the limited number of parameters thus 
far determined. We are not presenting our re- 
sults on this transition until further work is 
completed, at which time we shall present a more 
complete exposition of the experimental results. 
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ISOTOPIC SPIN ADMIXTURES IN Ge®*~Ga®*+Zn® 8 TRANSITIONS* 


W. P. Alford and J. B. French 
Department of Physics and Astronomy, University of Rochester, Rochester, New York 
(Received January 9, 1961) 


There has been recent discussion,’ stimulated 
by the theory of the conserved vector interac- 
tion,? concerning the goodness of the isotopic 
spin quantum number T in low-lying states of in- 
termediate-weight nuclei. Isotopic spin admix- 
tures have been measured by determining the 
Fermi contribution to a AJ=0 (J #0), AT=1 
8 transition which proceeds mainly via a Gamow- 
Teller interaction. 

The fact that the odd-odd nucleus Ga® has spin 
zero*** implies that both of the transitions 
oe™+,,Ga~+,,Zn® involve a Ad =0 (J =0) iso- 
topic-spin-forbidden pure Fermi transition. It 
then follows that the ft values® (log ft >6.8 for the 
first, logft='7.88 for the second) supply rather 
precise determinations of isotopic-spin admix- 
tures. The three ground states must have dif- 
ferent T values (since otherwise one or both of 
the transitions would be superallowed) and we 
can assume, without affecting the argument, that 
T=1, 2, 3, respectively, which is most probable. 
We note first that the relativistic correction to 
the vector 8 decay interaction does not contribute 
in a 0*+0* decay.* Then if mesonic effects are 
negligible or conserve isotopic spin, we see that 
the Ge® 8 decay proceeds via an admixture in its 
own ground state of the T=2 level which is the 
analog of the Ga® ground state. Similarly Ga® 
decays via an admixture in its ground state of 
the analog of the Zn® ground state. We label 
these amplitudes as a,, and a,,, respectively. For 
a transition of this variety with T, T,=T, Tz+1 
the ft value is given by 


(ff)"*=C la I(T-7,)(T+T, +1), 


where the constant C is determined by the O* 
decay’ (T=1, T,=0, a=1, ft=3100). We then 
find immediately that 


lap, !?<3x1074; |ag.!?=1.41075. 


These small values do not of course give the 
total admixture since other 0* states, which do 
not contribute to the 8 decay, should be admixed. 
But the admixtures determined here are of ex- 
actly the same nature as those determined less 
accurately in the (8,7) experiments: and, be- 
sides that, there are no experiments which de- 
termine the total T impurity. It is likely too that 
the total impurity is quite comparable with the 





values given for a single component since the 
density of low-lying 0* states which may be ad- 
mixed should be quite small.* Small values 
would also be expected from the calculations of 
MacDonald,° of Bouchiat,'° and of Iwao." 

It is probable that the ,,Ga*~,,Zn™ decay sup- 
plies an analogous example since the Ga™ ground- 
state spin is probably” 0*. In this case logft=6.6 
which would give, for Ga“, |a,,|7=4x10"*. A 
more certain determination of the Ga“ spin would 
be worthwhile. It would also be of interest to 
measure more accurately the ft value for Ge® 
and to locate the contributing 0* state in Ge®, 
Ga®*, and Ga™. 

We acknowledge a useful conversation concern- 
ing the (Sy) experiments with Dr. S. D. Bloom. 





*Supported in part by the U. S. Atomic Energy Com- 
mission. 

'S. D. Bloom, L. G. Mann, and J. A. Miskel, Phys. 
Rev. Letters 5, 326 (1960), and references therein to 
earlier work. 

*R. P. Feynman and M. Gell-Mann, Phys. Rev. 109, 
193 (1958); R. P. Feynman, Proceedings of the Tenth 
Annual International Conference on High-Energy Phys- 
ics, 1960 (Interscience Publishers, New York, 1960), 
p. 501. 

3J. C. Hubbs, W. A. Nierenberg, H. A. Shugart, 
and J. L. Worcester, Phys. Rev. 107, 723 (1957); 

J. L. Worcester, W. A. Nierenberg, R. Marrus, and 
J. C. Hubbs, Bull. Am. Phys. Soc. 2, 383 (1957). 

4A. Schwarzschild and L. Grodzins, Phys. Rev. 
119, 276 (1960). 

K. Way, F. Everling, G. H. Fuller, N. B. Gove, 
C. L. McGinnis, and R. Nakasima, Nuclear Data 
Sheets (National Research Council, Washington, D. C.). 

SA. Altman and W. M. MacDonald, Phys. Rev. 
Letters 1, 456 (1958). 

™R. K. Bardin, C. A. Barnes, W. A. Fowler, and 
P. A. Seeger, Phys. Rev. Letters 5, 323 (1960). 

®The argument could fail if the other 0* states were 
much more strongly admixed. On the other hand, the 
admixing Coulomb matrix element corresponding to 
Qs, is 4x 10-AE, where AE is the excitation energy of 
the first 0+, T=3 state in Ga®, and, assuming that 
AE is a few Mev, this is probably not an anomalously 
small value. 

*w. M. MacDonald, Phys. Rev. 110, 1420 (1958). 

°C. C. Bouchiat, Phys. Rev. 118, 540 (1960). 

11g, Iwao, University of Rochester thesis, 1960 (un- 
published). 

27, H. Jacobi, H. A. Howe, and J. R. Richardson, 
Phys. Rev. 117, 1086 (1960). 
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DEUTERON ELECTRIC QUADRUPOLE MOMENT* 


J. P. Auffray 
Institute of Mathematical Sciences, New York University, New York, New York 
(Received December 19, 1960) 


The nuclear radio-frequency spectra of molecu- 
lar hydrogen and deuterium, first observed by 
Rabi et al.,’ were reinvestigated by Ramsey et al.’ 
From these experiments the value of the quad- 
rupole coupling constant eQq, in the v=0, J=1 
vibrational-rotational level of the electronic 
ground state of D,, was determined with an esti- 
mated accuracy of better than 1 part in 10°. 

An accurate value for the deuteron electric 
quadrupole moment Q can be derived from this 
result if g, the electric field gradient along the 
molecular axis at one deuteron, is evaluated - 
theoretically. This was done by Nordsieck® and 
later by Newell* using a trial molecular wave 
function which gave a molecular binding energy 
Dz, of 4.566 ev.® 

The present author has calculated gq from a 
trial function which gives Dg =4.728 ev. This is 
much closer to the experimental value® of 4.747 
ev. The trial function was constructed in the 
framework of the Born-Oppenheimer approxima- 
tion. Its electronic part consisted of an expan- 
sion in the electronic coordinates and the inter- 
electronic distance, of the type first investigated 
by James and Coolidge,’ and recently by Kolos 
and Roothaan.® 

The value of g obtained from this function yields 
Q =2.82x107?" cm”. This is about 3% larger than 
the currently accepted value based on Newell’s 
result. 


The discrepancy is larger than the estimate 
given by Newell.* Newell, however, obtained this 
estimate by an admittedly crude procedure. In 
addition it should be noted that g, in these cal- 
culations, is not stationary. Consequently no 
rigorous estimate of the error in qg can actually 
be made. 

The details of this calculation will be given ina 
forthcoming publication. 





*The work presented in this paper is supported by 
the AEC Computing and Applied Mathematics Center, 
Institute of Mathematical Sciences, New York Univer- 
sity, under contract with the U. S. Atomic Energy 
Commission. 

‘J. M. B. Kellogg, I. I. Rabi, N. F. Ramsey, Jr., 
and J. R. Zacharias, Phys. Rev. 57, 677 (1940). 

2H. G. Kolsky, T. E. Phipps, Jr., N. F. Ramsey, 
and H. B. Silsbee, Phys. Rev. 87, 395 (1952). 

3A. Nordsieck, Phys. Rev. 58, 310 (1940). 

‘G. F. Newell, Phys. Rev. 78, 711 (1950). 

5Converted in electron volts from Newell’s value in 
atomic units (see reference 4), using 1 a.u. = 27.210 
ev. 

®°G. Herzberg, Spectra of Diatomic Molecules 
(D. Van Nostrand and Company, Inc., Princeton, 
New Jersey, 1950). 

"H. M. James and A. S. Coolidge, J. Chem. Phys. 
1, 825 (1933). 

‘Ww. Kolos and C. C. J. Roothaan, Revs. Modern 
Phys. 32, 219 (1960). 








POSSIBLE EXISTENCE OF A NEW K’ MESON* 


J. Tiomno, A. L. L. Videira, and N. Zagury 
Centro Brasileiro de Pesquisas Fisicas and Faculdade Nacional de Filosofia, Rio de Janeiro, Brazil 
(Received December 15, 1960) 


It has been experimentally established that the 
A particles produced in the reaction 


p+n~-K°+A (1) 


are strongly peaked backward in the c.m. system. 
Pais’ has tried to explain this result by assum- 
ing the existence of a KKn interaction and that the 
main contribution to reaction (1) comes from the 
Feynman graph of Fig. 1, assuming x* to bea 

K* meson. He has obtained a reasonable qualita- 


120 





tive agreement with the available experimental 
data at 1.1 Bev. However, his scheme, supposing 
opposite parities for K* and K°, leads to many 
problems of difficult solution, arising from the 
failure of charge independence of the 7 interaction. 
One such difficulty results from the lack of a 
KKr° interaction, which gives rise to a mass of 
charged n’s smaller than that of 7°.” 

In the present paper we shall try, maintaining 
the convenient features of the Pais scheme, to 
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TT i 
FIG. 1. Feynman 7. 
diagram of the assumed | x* 
predominant contribu- 
tion to K°A production. ——P A. _ 
p » A 


overcome such difficulties by assuming the ex- 
istence of another K’ meson of the same hyper- 
charge and isospin (1/2) as those of the ordinary 
K meson but of opposite parity.* 

Reaction (1) proceeds now mainly via the same 
graph as in Fig. 1 but with the intermediate ,* 
being a K’+ meson instead of K*. 

In order to compute the differential cross sec- 
tions of reaction (1) we add to the usual d’Espagnat 
Prentki interaction Hamiltonian density a KK’7 
and a AK’N term: 


2m ,f (KTK’ +K’TK)-7+g(K’AIN +H.c.), (2) 


where my is the mass of the K meson and I is 1 
or y, according to the parity of K’. 

Using the Born approximation, we get for the 
differential cross section: 


do/dQ 
1 "9h cosé tm ,m/E .E p 
z z , 
[mi +m" = B= 2E El +0 Uy cos@) |* al 


, (3b) 
T 





A = 4m (f*/4a Vg /40)P EE ANE +E )’p 
where py is the mass of the K’ meson and @ is the 
angle between 7 and A in the c.m. system. The 
plus sign corresponds to ['=1 (K’ scalar), and 
the minus to I =y, (K’ pseudoscalar). 

Although the use of the Born approximation is 


Table I. 
fitting the experimental data to the theoretical curve. 








not well justified, we compared (3) with the avail- 
able experimental distributions at several values 
of the energy with the best statistics.*»® In this 
way we determined in each case the value of yu 
which gave the best fit. These values of yu (or pu”) 
are given in Table I, for both scalar and pseudo- 


scalar K’. The values of fg/4n are also indicated 
for the cases where total cross sections were 
available. The most probable value was found to 


be 
(4a) 
(4b) 


The last possibility can be immediately excluded 
since it would lead to a fast decay of K into K’ +7. 
Thus we conclude that the K’ (if it exists) is 
scalar, K being pseudoscalar in agreement with 
other experimental indications.°® 

In Fig. 2 all experimental data are plotted to- 
gether to give an idea of the goodness of fit with 
a unique mass in the case of scalar K’. In order 
to do this we introduce the notation (c.m. system): 


pu =605+ 49 Mev 
wu =176+46 Mev 


(K’ scalar), 


(K’ pseudoscalar). 


x(6,E)=2E E (1 +v_v,-C0s6) -(m_*+m,*), (5a) 
- do/dQ L(6, E)d2 
2°" (6,2)= Eyre, ance EN (5b) 
118, B)=1+9 9, e080 +m m pF AF es (5c) 
Thus Eqs. (3) imply that 
z(6,E) =u" +x(6,E), (6) 


for all energies (E =E,,!ab) and angles. 

The x and z coordinates for each experimental 
point plotted were obtained from (5) using p =605 
Mev. The values of (1/c)do/dQ and cosé@ used 
were average values in the intervals in which the 
angular distribution histograms, represented as 
functions of cos@, were divided. The straight 


Values of the mass of the hypothetical K’ meson and the product of the coupling constants obtained by 














tT,» a K’ scalar K’ pseudoscalar 
(Bev) u (Bev) fg/4n yu? (Bev?) fig/4n 
0.902? 0.90 +0. 20 see 0.048 +0. 034 vee 
0.907? 0.66 +0. 12 ia 0. 063 +0. 035 tee 
: 0.910° 0.45 +0.11 0.33 +0.07 -0.012 +0. 029 0.70 +0.10 
0.960° 0.49 40.11 0.43 +0.09 0.011 +0. 099 0.91 +0.40 
1. 200° 0.59 +0.11 0.31 +0. 06 0.034 +0. 039 0.55 +0.08 
1.300° 0.80 +0. 12 0.48 +0.09 0.103 +0.'085 0.69 +0. 08 











bgee reference 4. 
“See reference 5. 








«i lab ig the kinetic energy of the pion in the lab system. 
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lines indicated correspond to Eq. (6) for values 
of » =605, and 605+ 49 Mev. The fit of the ex- 
perimental points in Fig. 2 with these theoretical 
lines seems to be very good except for the points 
to the right which correspond to cos@~1, as 
should be expected. These points also have 
poorer statistics. The fit in the case of pseudo- 
scalar K’ was not so good. 

Actually the experimental points plotted are 
only appropriate for the comparison with the 
theoretical line for u = 605 Mev since z,,,, de- 
pends on up. However, for the two other values 
(605+ 49 Mev) the experimental points would be 
only slightly shifted. If we interpret these re- 
sults as indicating the existence of the K’ meson, 
this meson should have a mass 4 >my+m, for 
otherwise it would have been already observed. 
Then K’ should be very short-lived since the in- 
teraction leading to K’~K+z is strong. However, 








TT K' K K 
—---- ----- ae eee ae eee 
| ' 
|K TT 
| | 
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~ 
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FIG. 3. Diagrams of (a) 7 production of K’; (b) 
scattering of K producing a K’. 
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FIG. 2. Comparison of the 
experimental data on the angu- 
lar distribution of A at several 
energies with theoretical curves 
(dashed lines). The values of 
the mass of the hypothetical K’ 
meson are indicated. The co- 


4- 0,902 Bev 


0 - 0.907 Bev 





°o-09! Bev 
epoaen ordinates x and z are defined 
@-.2 Be by Eqs. (6). 
A-t3 Bev 
_ eee yaaa a 
10 L2 


it cannot be excluded that the K’ meson does not 
e..3t, representing only a pair of 7 and K, which 
are exchanged between the baryons and a KKnaa 
vertex. The existence of the K’ can be tested in- 
directly by the study of the reaction 


n+N—K+7+A (or £), (7) 


if the graph of Fig. 3(a) gives a significant con- 
tribution. Thus there should be a characteristic 
peaking in the energy distribution of A (or 2) in 
the c.m. system at a point depending on the mass 
of K’. The values (4a) should not, however, be 
taken very seriously due to the approximations 
made which are more satisfactory for higher 7 
energies. As seen from Table I, larger values 
of » are obtained at higher energies. Also in 
this peak the A (and 2) should have some back- 
ward preference and strong An angular correla- 
tions should be observed. The most spectacular 
effects of similar kind would occur, however, in 
the reaction 


K+N~K+m+N, (8) 


due to the graph of Fig. 3(b). 

A more detailed version of the present paper 
will appear in Anais da Academia Brasileira de 
Ciéncias. 





*The essential results of this paper were reported 
by J. Tiomno at the 1960 Rochester Conference during 
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the discussion of Gell-Mann’s paper. This work was 
supported in part by the Conselho Nacional de Pes- 
quisas, Rio de Janeiro. 

‘A. Pais, Phys. Rev. 112, 624 (1958). 

23. Tiomno (unpublished). 

3See also M. Gell-Mann, Proceedings of the 1960 
Annual International Conference on High-Energy Phys- 
ics at Rochester (Interscience Publishers, New York, 
1960). 











4M. Schwartz, Proceedings of the 1960 Annual In- 





ternational Conference on High-Energy Physics at 
Rochester (Interscience Publishers, New York, 1960). 
We are thankful to Dr. Schwartz who kindly supplied 
us with the data. 

°F. Eisler et al., Nuovo cimento 10, 468 (1958). 

®°M. Moravesik, Phys. Rev. Letters 2, 352 (1959); 
C. M. Meltzer, Ph. D. thesis, Duke University, 1960 
(unpublished). 








SOLAR-PRODUCED COSMIC RAYS NEAR THE NORTH AND SOUTH POLES* 


M. A. Pomerantz, S. P. Duggal, and K. Nagashima 
Bartol Research Foundation of The Franklin Institute, Swarthmore, Pennsylvania 
(Received January 16, 1961) 


The sudden arrival of solar-produced cosmic 
radiation of sufficient energy to be detectable at 
ground stations, a relatively rare phenomenon, 
has for the first time been observed simultane- 
ously by the two neutron monitors which are 
currently closest to the North and South geomag- 
netic poles. Two events, each of which appeared 
to be associated with a solar flare of Importance 
3+, occurred on November 12 and 15, 1960, re- 
spectively, the former displaying some remark- 
able features. 

The locations of the Arctic and Antarctic sta- 
tions operated by the Bartol Research Foundation 
are listed in Table I, and the relevant data are 
plotted in Fig. 1. The intensity is expressed in 
terms of percent deviation from the preflare 
mean (0000 to 1200 U. T. on November 12) of 
barometric pressure-corrected counts, recorded 
during 15-minute intervals. It is apparent that, 
on November 12, the increase at both stations 
manifested an unusual double-hump structure, 
whereas the November 15 event is similar to 


previous flare-associated intensity enhancements. 


It is interesting to note that, on November 12, 
greater increases in intensity were recorded at 
the Antarctic station during the period of the 


first hump. On the contrary, the increases at 
both stations were approximately equal there- 
after, although there is some uncertainty in the 
comparison of the amplitudes of the second maxi- 
mum on this date, owing to a 45-minute gap in 

the record during the weekly routine Ra-Be source 
runs. 

As a quantitative measure of the relative in- 
tensity enhancement at the two polar stations, 
the ratios of the increases, AJ(McMurdo)/A/(Thule), 
are also plotted in Fig. 1. The 12-hour preflare 
mean counting rate on each of the two respective 
dates is selected as reference level for the per- 
cent deviation, AJ. This ratio attains a maximum 
about 30 minutes after onset, and approaches 
unity after approximately 5 hours on November 
12. The first peak in the ratio AJ(McMurdo)/ 
Al(Thule) occurs during the rising portion of the 
first hump, while the second peak occurs between 
the two humps. A single peak of shorter duration 
also occurs during the rising portion of the event 
of November 15. 

The onset times, and times and amplitudes of 
maxima, are listed in Table II. Onset is defined 
as the first 2-minute interval during which, to- 
gether with five consecutive equal intervals, the 


Table I. Locations of the polar neutron monitor stations. 








Geographic coordinates Altitude Geomagnetic coordinates 

Station Latitude Longitude (meters) Latitude Longitude 
Thule N 76.6° W 68. 8° 260 +88. 0° ..¥ 
McMurdo S 77.9° E 166.6° 48 -79.0° 294. 3° 
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intensity exceeds by 20 that which prevailed 
during one hour immediately preceding the flare. 
Maximum is defined as the 2-minute interval 
during which the intensity exceeds by 2c that 
during the intervals on either side. Recorded 
counts uncorrected for pressure were selected 
for these time determinations. The indicated un- 
certainties arise solely from statistical consider- 
ations. Since, in the case of McMurdo, three 
equal maxima satisfying the above conditions oc- 
curred on November 12, the total period including 
all of these has been designated as the time of 
maximum. The amplitudes of the maxima are 
intensities expressed as percent deviations from 


NOV. 15, 1960 


the mean preflare intensity on November 12, 
pressure corrections having been applied. 

Table II reveals that, within the experimental 
uncertainties, the onset times of the November 
12 and 15 events are the same at both stations. 
However, in the former case, the initial inten- 
sity maximum is attained somewhat later, and 
is greater, at McMurdo than at Thule, in con- 
trast with the November 15 increase. 

The observation of these increases at polar 
stations is incompatible with simple impact 
zone theory. Magnetized clouds in the inter- 
planetary space, perhaps associated with the 
magnetic storm which commenced at 1350 U. T. 


Table Il. Comparison of the two events as observed at the polar stations. 








Date Hump Times and amplitudes Thule McMurdo 
November 12, 1960 Onset time, U. T. 1357 +1 1355 +1 
I Time of maximum, U. T. 1608-1610 1624-1640 
Amplitude of maximum 103 % 118% 
0 Time of maximum, U. T. 1958-2000 
Amplitude of maximum 135% see 
November 15, 1960 Onset time, U. T. 0239 +1 0239 +1 
Time of maximum, U. T. 0436-0438 0436-0438 
Amplitude of maximum 83% 85 % 
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on November 12,” may have been responsible 
for the deflection of the solar-produced particles 
into the polar regions. 

The double-hump structure of the increase in 
the intensity of the nucleonic component has not 
been observed previously. The most remarkable 


feature is the difference in the ratio, AJ(McMurdo)/ 


Al(Thule), in the first hump as compared with the 
second. For an isotropic distribution of cosmic- 
ray intensity, the ratio is approximately unity. 
Hence, a value of the ratio exceeding unity is 
indicative of an anisotropy in the intensity dis- 
tribution of the solar-produced cosmic rays in 
the vicinity of the earth. As is seen in Fig. 1, 
the anisotropy prevailed during the first hump 
and isotropy has been almost established near 
the beginning of the second hump. A possible 
interpretation may be that, during the first hump, 
the earth was located near the boundary of the 
modulating region, which later completely sur- 
rounded the earth. 

It is interesting to note that the anistropy was 
evident only during the exceptionally long dura- 
tion (about 6 hours) of the Ha flare. However, 
the significance of this coincidence is not clear, 


especially since a Forbush-type decrease com- 
menced® at approximately the same time as the 
termination of the optical flare. 

It is a pleasure to acknowledge the important 
contributions to the program by Hugo A. C. Neu- 
burg and Jamie Chapman who are in charge of the 
stations at McMurdo and Thule, respectively. 
Appreciation is expressed to the Geophysics 
Research Directorate for maintaining the Arctic 
Station at its Polar Research Facility, operated 
by the American Geographical Society under con- 
tract with the U. S. Air Force. The Antarctic 
Station is part of the U. S. Antarctic Research 
Program. 





*This work was supported in part by grants from the 
National Science Foundation, and was also assisted by 
the Office of Naval Research. 

‘A, Schliiter, Z. Naturforsch. 6a, 613 (1951). 

*We are grateful to Mr. Robert Clements for tracings 
of the magnetograms recorded at the New Zealand Sta- 
tion at Scott Base. 

3K. G. McCracken, D. C. Rose, G. Schwachheim, 
R. Palmeira, and T. Thambyahpillai (private com- 
munications). 





CLOUD-CHAMBER OBSERVATIONS OF PRIMARY COSMIC-RAY ELECTRONS* 


James A. Earl 
School of Physics, University of Minnesota, Minneapolis, Minnesota 
(Received December 16, 1960) 


This Letter reports new measurements of the 
flux of high-energy (0.3 Bev to 3.0 Bev) electrons 
and gamma rays at the top of the atmosphere. In 
1949, an experiment using essentially identical 
equipment (a multiplate cloud chamber carried 
to high altitudes by a balloon) was performed by 
Critchfield, Ney, and Oleksa.’ While no evidence 
for primary cosmic-ray electrons was found in 
that experiment, improvements in cloud-chamber 
technique and higher balloon altitude capability 
now make it possible to demonstrate the existence 
of a small flux of primary electrons. 

Table I summarizes the data pertaining to the 
balloon flight and to the characteristics of the 
cloud chamber. Electrons and gamma rays were 
identified by the characteristic electron- photon 
cascade showers which they produced in the lead 
plates of the cloud chamber. These showers, 
which appear as a sequence of cones of minimum- 
ionizing nonpenetrating tracks emerging from the 





lower surfaces of the lead plates, are readily dis- 
tinguished from events produced by nuclear inter- 
actions. The selection criterion that the axis of 
accepted showers had to pass through the illumi- 
nated areas of both the top and the bottom of the 
cloud chamber ensured that the axis passed 
through all five plates. In passing through this 
thickness of material a nuclear shower would 
almost certainly reveal its nature through the 
presence of penetrating particles and/or heavily 
ionizing evaporation tracks. Figure 1 shows a 
typical example of a shower produced by an elec- 
tron. The incident electron appears in the top sec- 
tion of the chamber as a single minimum -ionizing 
track which lies on the axis of the shower. Since 

it is very unlikely that a background track or a 
back-scattered track from the shower would coin- 
cide accurately with the shower axis, it appears 
certain that events such as the one shown in 
Fig. 1 are initiated by high-energy electrons. 
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Table I. Data on the balloon flight and the cloud 
chamber. 





Balloon flight 





Date: May 12, 1960 


Minneapolis, Minnesota 
(geomagnetic latitude 55°N) 


Location: 


12 hours 


4 to 6.5 g cem~. 
Average: 4.5 g cm~* 


Time at ceiling: 
Pressure altitude: 





Multiplate cloud chamber 





Number of lead plates: 5 


0.6 cm-7.5 g em~* 
(1.1 radiation lengths) 


Sensitive time per picture: (0.19 +0.01) sec 


Thickness of plates: 


Geometric factor for 
region bounded by illu- 
minated areas of top 


and bottom: (33.5 +1.5) cm? sr 





Shower events in which there was no trace of a 
track in the top section even though the axis was 
well illuminated were assumed to be initiated by 
gamma rays. There is no reason to believe that 
the upper section was insensitive at any time dur- 
ing the flight. The procedure used to determine 
the energies of the electrons and gamma rays is 
identical to that used by Roe and Ozaki? and is 
based on the Monte-Carlo calculations of Wilson.* 
Figure 2 and Table II summarize the data taken 
while the balloon was at ceiling altitude. Figure 2 
shows the differential energy spectra for the elec- 
trons and gamma rays. While it is certain that 
some electrons were incident upon the chamber, 
it is not possible to estimate the number of pri- 
mary electrons without first knowing the number 
of secondary electrons produced in the atmos- 
phere above the chamber. An upper limit for the 
number of secondaries can be computed from the 
observed number of gamma rays under the as- 
sumption that all the gamma rays arise from nu- 
clear interactions above the chamber. The dashed 
histogram in Fig. 2 shows the energy spectrum 
of secondary electrons computed from the ob- 
served gamma-ray spectrum under the following 
assumptions: (a) The gamma rays arise only as 
decay products of 7° mesons; (b) two charged 7 
mesons are produced for every 7° meson; (c) the 
charged 7 mesons and the 7° mesons are produced 
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FIG. 1. A cloud-chamber picture of a shower pro- 
duced by a high-energy electron. The incident electron 
is visible in the top section of the cloud chamber. 


with the same energy spectrum; and (d) the Michel 
parameter which describes the 1. -meson decay 
spectrum has the value pm =0.75.* Under these 
assumptions, most of the secondary electrons 
come from the decay of » mesons created when 
the charged 7 mesons decay (a small contribution 
from pair production by the gamma rays has also 
been included). In the simplest approximation, 
the number of gamma rays and secondary elec- 
trons should be equal because there are two gam- 
ma rays per 7° meson and two charged 7 mesons 
(each of which should ultimately give an electron) 
per 7° meson. A more detailed analysis shows 
that less than half of the » mesons decay before 
they reach the level of the balloon and that many 
of the electrons from these have energies so low 
that they could not be primary cosmic rays ar- 
riving at the latitude of Minneapolis. The net re- 
sult is that, in the energy region accessible to 
primary cosmic rays (above 0.5 Bev), only 1.5 

+ 0.4 secondary electrons are expected, while the 
total number of electrons observed was 11. There- 
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FIG. 2. Differential energy spectra for electrons 
and gamma rays. The dashed spectrum represents 
the secondary electrons that would be expected under 
the assumption that all the gamma rays are secondaries. 


fore, most of the observed electrons must be 
primary cosmic rays. The pronounced peak in 
the electron energy spectrum between 0.5 Bev 
and 1.0 Bev is further evidence in favor of this 
conclusion because such a peak would appear if 
primary cosmic rays were excluded below the 
geomagnetic cutoff energy. (The cutoff energy for 
electrons at Minneapolis is about 0.7 Bev.)® While 
there is a low-energy (about 0.25 Bev) instrumen- 
tal cutoff arising from the fact that very small 
showers are not easily recognized, more elec- 
trons than are observed would appear in the inter- 
val from 0 to 0.5 Bev if the spectrum of the ob- 
served electrons had the monotonic increase with 
decreasing energy which seems to be a charac- 
teristic of the spectrum of secondary electrons. 
A reasonable upper limit on the geomagnetic 
albedo (upward moving particles produced in the 
atmosphere which are returned to the earth by 
the geomagnetic field) is the upward flux at the 
point of observation. Since no upward moving 
electrons with energies greater than 0.5 Bev were 
recorded during the flight, it is very unlikely that 
the observed electrons arise from this source. 
(Most of the » mesons whose decay electrons con- 
stitute the upward albedo would be expected to 
decay before they reach the level of the balloon.) 
The value computed from the data in Tables I 
and II for the flux of electrons incident on the 
cloud chamber with more than 0.5-Bev energy is 


| 


L, =(32+10) particles m~™ sec™ sr 


No more than (13+ 10)% of this flux is due to 


Table Il. Tabulation of particles whose paths passed 
through illuminated areas of both top and bottom of 
cloud chamber. 





Total number of electrons 11 
Total number of gamma rays 6 
Number of minimum-ionizing particles which 
penetrated all 5 plates 284 
Number of penetrating particles corrected 
for nuclear interactions 380 
Total number of pictures at ceiling 541 





secondaries if the assumptions given earlier are 
correct. Since the correction for secondaries is 
small and since it is partially offset by absorp- 
tion of the primary electrons, this flux is a close 
approximation to the true flux of primary elec- 
trons at the top of the atmosphere. The flux of 
protons determined from the number of minimum- 
ionizing penetrating particles selected by the same 
geometric criterion that was applied to shower 
axes is 


1,=(1100¢ 100) particles m~? sec™ sr™. 


A correction based on a mean free path for nu- 
clear interactions of protons in lead of 140g cm™~ 
has been applied in computing this flux. The ra- 
tio of the electron flux to the proton flux is 


1/1, =(3 #1) %. 


In 1949 Critchfield, Ney, and Oleksa set an upper 
limit on this ratio of 0.6%. The reasons that the 
above-measured value is larger than this upper 
limit are the following: (1) The minimum energy 
of the accepted showers was larger in the earlier 
experiment (1.0 Bev vs 0.5 Bev); (2) the proton 
flux was higher in 1949 than in 1960 (2200 p m™ 
sec! vs 1100 p m™ sec™); and (3) the thickness 
of the air above the balloon was less in the pre- 
sent experiment (4.5g cm~ vs 15g cm™?). The 
conclusion that primary electrons exist was made 
possible in the present experiment by (1) a clear 
separation of electrons and gamma rays which 
enabled a rigorous calculation of background ef- 
fects to be made, and (2) the high altitude obtained 
which ensured that the flux of secondary electrons 
would be small. 

Calculations by Ginzburg® indicate that synchro- 
tron radiation from high-energy electrons would 
account for the observed intensity of galactic 
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radio emission if the density of electrons with 
energy greater than 1 Bev is on the order of 3 

x 1075 cm. The density computed from the 
measured electron flux is 1.3107” cm™, which 
appears to be more than enough to account for the 
observations of galactic radio noise. 

An upper limit for the flux of electrons arising 
from nuclear interactions of cosmic rays with 
interstellar hydrogen gas can be obtained from 
the data which were used earlier to determine 
the number of secondary electrons produced in 
the atmosphere. Although the average thickness 
of matter traversed by a cosmic ray within the 
galaxy is about 1g cm™ (this thickness is sug- 
gested as an upper limit by considerations based 
on the relative abundances of protons and heavy 
nuclei in the cosmic-ray beam), the number of 
electrons produced in this thickness of interstel- 
lar hydrogen is about the same as the number 
produced in the 4g cm~ of air above the balloon. 
This occurs because the interaction probability 
per g cm™ for hydrogen is about twice that for 
air and because all of the » mesons produced in 
space decay while only half of those produced in 
the atmosphere above the balloon do so. (The 
flux of gamma rays from interstellar nuclear 
collisions is very small compared to the electron 
flux because the time that the electrons are trap- 
ped by galactic magnetic fields is large compared 
to the time required for the gamma rays to leave 


the galaxy on straight paths.) If the mean distances 
traversed in the galaxy by protons and electrons 
are equal, the electron flux arising from nuclear 
interactions in interstellar space is no larger than $ 
the flux arising from nuclear interactions in the 
atmosphere above the balloon. Since this flux was 
estimated to be only 13 % of the total electron flux, 
it appears unlikely that all of the observed elec- 
trons arise from nuclear interactions of cosmic 
rays with interstellar hydrogen. 

The author would like to thank Professor E. P. 
Ney for his interest in this investigation and Pro- 
fessor J. R. Winckler for his comments on the 
manuscript. 
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MEASUREMENT OF THE ANOMALOUS MAGNETIC MOMENT OF THE MUON 


G. Charpak, F. J. M. Farley, R. L. Garwin,* T. Muller, J. C. Sens, V. L. Telegdi, t and A. Zichichi 
CERN, Geneva, Switzerland * : 
(Received January 16, 1961) ; 


By storing polarized » mesons in a magnetic 
field for as long as 1000 cyclotron periods it has 
been possible to measure directly the anomalous 
magnetic moment. We find the anomaly in agree- 
ment to within 2% (that is, 2x107° accuracy on 
the total magnetic moment) with that expected 
from the quantum electrodynamics of a Dirac 
particle. 

At present the muon appears to be a heavy elec- 
tron with no interactions except the electromag- 
netic and the weak. This concept gives no ex- 
planation for the muon-electron mass difference, 
but allows the muon magnetic moment to be cal- 
culated from the Dirac equation and quantum 
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electrodynamics as'~*® 
Lh =g(e/2Mc) x(h/2), (1) 


with g=2(1+a), the anomalous part of the moment, 

a, being 

G,,, = (g- 2)/2 =(a/2n) + 0.75(a*/n*) ++ + -= 0.001165, 
(2) 

with a” =hc/e? = 137.04, the fine-structure con- 

stant of atomic physics. The first term in Eq. (2), 

a/2n, arises from the emission and reabsorp- 

tion of single photons, which alter the Dirac 

moment in two distinct ways, by the muon recoil, 

and by the muon spin-flip during the period of 
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existence of the virtual photon. 

A “breakdown of quantum electrodynamics,” 
for instance a cutoff on the photon propagator at 
energy Am,c*, will modify Eq. (2) to (approxi- 
mately)’ 


a=(a/2m)(1 - $A) +++>, (3) 


while the existence of a “fundamental length” L, 
which might be expected to introduce cutoffs on 
the muon form factor and muon propagator, as 
well as on the photon propagator, will also intro- 
duce® a correction of form Eq. (3), but with 
A~h(4McL)"?. 

Thus an experiment to a precision of 10~° on g, 
or 1% of the a to be expected from Eq. (2), is sen- 
sitive to the behavior of the electromagnetic field 
alone to a distance of ~3 x107* cm, and to the 
existence of a fundamental length as small as 
7.5x107'5 cm. Since there is nothing special about 
a photon, other fields with which the muon inter- 
acts would contribute in just the same way to the 
anomalous moment, scalar fields giving only the 
recoil effect while vector fields yield also the 
spin-flip term. Thus a 1% experiment is sensi- 
tive to the existence of a coupling constant of order 
of magnitude 10~ to a field of mass equal to that 
of the muon, or 10~° for a field of mass equal to 
that of the nucleon. Of course, these coupling con- 
stants depend on the type of coupling assumed. 

The magnetic moment of the muon has been 
measured’ to high accuracy in absolute units, 
as has the muon mass’*»"* (needed for computing 
the muon magneton in order to calculate g). But 
to determine a to 1% of the expected value would 
require the mass to <107° accuracy, which seems 


FIG. 1. General 
plan of the 6-meter 
magnet. M: bending ' 
magnet; Q: pair of 
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impossible even with the advanced techniques 
employed so far. 

Fortunately there exists a direct way to com- 
pare the magnetic moment with the muon mag- 
neton using the principle employed for the elec- 
tron by Louisell, Pidd, and Crane” and already 
applied very roughly to the muon by Garwin et al.’* 
If muons are made to circle in a magnetic field, 
the spin turns (1+ya) times as fast as the mo- 
mentum vector, with the result that after time ¢ 
in a magnetic field B the relative angle between 
spin and momentum is changed by the amount 
6 =aBw,t, where w, is the cyclotron frequency 
for zero-energy muons in unit magnetic field.**»*® 
We wish to report here the first precision meas- 
urement of the muon anomalous moment by this 
method.*® 

The muon mean life, 2.2 microseconds for a 
muon at rest, sets a limit to the storage time ¢ 
which may be used. In fact, if the only loss is 
due to decay, the optimum storage time is two 
mean lives. In such an experiment the problems 
are therefore to know the direction of the polar- 
ization of the muons (relative to the momentum) 
before and after the precession time ¢, to know 
the mean field B, and to know the time ¢ spent in 
the field, all to an accuracy considerably better 
than 1%. In fact, we record the polarization 
angle @ as a function of storage time ¢, since in 
our method muons have a range of storage times 
from 2.0 usec to about 6.5 psec. 

In our scheme shown in Fig. 1, longitudinally 
polarized muons, formed by forward decay of 
pions in flight inside the cyclotron, are focused 
by a bending magnet and quadrupole pair into the 
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entrance channel of the “6-meter magnet” and 
cross the center line (y = 0) at roughly 90°. At 
this point the muons are moderated by a beryllium 
block to a mean energy corresponding to a radius 
of curvature R of 19 cm in the magnetic field B, 
of 15.8 kgauss. The field is not uniform but has 
at any x in the horizontal median plane the form 


Bly) =B,(1 +ay + by? +cy*), (4) 


in which the term ay produces a walking towards 
+x with a “step-size” per turn, S=2R%a, where 
R is the orbit radius. The term by? with <0 
adds vertical focusing, and the term cy*® changes 
the step-size to 


S =nR?(a +0.75cR?), (5) 


allowing one to choose c to give constant S for all 
R in the neighborhood of the mean radius. Such 
a field with c=0 and with a and b constant in x 
has been used at this laboratory with a smaller 
magnet to store muons for 25 turns and to meas- 
ure their electric dipole moment.” 

In the injection region, S~2 cm/turn to allow 
considerable tolerance in alignment. A slow 
transition is then made to a weaker gradient giv- 
ing a storage region with S=4 mm/turn, while 
in the ejection region the step-size is 10 cm/turn 
to allow the particle to cross the exit face of the 
magnet nonadiabatically and so to be ejected from 
the field. In practice the desired field shapes 
were obtained by shimming a magnet with flat 
pole faces with many layers of 0.5-mm steel 
shims. About 300 kg of shim is inside the mag- 
net vacuum chamber. 

Hall plate measurements showed that the nec- 
essary field form had been obtained, and the de- 
sired motion of the orbits was verified by study- 
ing the flux through a search coil of 18-cm radius, 
the flux through the walking orbit being an adia- 
batic invariant. The field is calibrated by proton 
resonance. A change of field of 0.1% from one 
end of the storage region to the other would pro- 
duce a 2.5-cm sidewise displacement of the orbit 
and would lead to a considerable loss in intensity. 
By correcting local inhomogeneities by shims as 
thin as 0.03 mm, we achieved a field in which the 
orbits do not wander by more than +1 cm later- 
ally. 

The timing to ~ 0.1% accuracy is done with a 
“digitron” developed from that of Swanson and 
Lundy.” Additional dead-time circuits are in- 
corporated to remedy some minor defects in the 
original instrument and the timing is now referred 
to a 10-Mc/sec crystal oscillator. It is of great 
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importance not only to time every muon to 1%, 
but also to reduce below 0.1% the probability of 
assigning the wrong injection count to an emerging 
uu, etc. The electronics which accomplishes this 
is too complex to discuss in full here, but it has a 
timing accuracy of 10~® sec and a flat spectrum 
between random events to «14%; and it rejects 
completely all cases of doubtful parentage or 
progeny of no matter how high order. 

The muons emerging from the magnet are 
stopped in a field-free region (<0.1 gauss) ina 
nondepolarizing, nonconductive target (methylene 
iodide). The electronics stores events associated 
with “backward decay electrons” (6-6’ coinci- 
dences with appropriate gating) in channels 1-50 
of a CDC “pulse-height analyzer,” as a function 
of the storage time ¢ of the parent muon in the 
magnet (2.0-6.5 usec); similarly events associ- 
ated with forward decay electrons in telescope 
7-7’ are classified according to storage time in 
channels 51-100. The muon spin direction is 
flipped in successive runs through + 90° in the 
1.0 usec following its arrival by means of a 
pulsed vertical magnetic field produced by an 
aluminum -tape coil wound on the target. From 
the counts cy, and cy_ in the mth channel for + 90° 
flipping, one can compute independently for each 
channel an asymmetry 


A ale -¢, e+e, ) =A sin(aBwt ), (6) | 


where A is proportional to the polarization of the 
incident beam and ¢, is the storage time cor- 
responding to the mth and (m + 50)th channel. 

Figure 2 shows the observed data as well as 
the computed curve of the form (6) which fits it 
best. In actual practice the muons available for 
injection are not exactly longitudinally polarized. 
There is a transverse component which varies 
with range (~5°/g cm~*), crossing zero at the 
peak of the range curve. This has been measured 
and included in the calculation, as has the angu- 
lar distribution of decaying muons with respect 
to the normal to the polarization analyzer. The 
parameter (g-2)/2 for this best fit is 
Soxp =a, (0.983 + 0.019) = 0.001145 + 0.000022, (7) 
and the x’ for this fit is 13.5, in comparison with 
9.7 expected.*® 

Of the 1.9% error quoted in Eq. (7), 1.7% is 
statistics, 0.6% comes from the uncertainty in 
the measured direction of the ejected beam, 0.6% 
from uncertainty in the part of the range spec- 
trum (and hence transverse polarization) of the 
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FIG. 2. Curve A (left-hand 0,|¢ sec 
scale for ordinate): Storage- 
time distribution of muons that 
stop in the polarization analy- 4r 
zer and give rise to decay 
electrons. Curve B (right- 3h 

hand scale for ordinate): The 

sinusoidal curve [see formula 2 
(6)] represents the best fit to 
the measured variation of the 
asymmetry with storage time. 
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the point near zero time, 
which is a very important 
datum for the fitting. 


injected beam, and 0.3% from the statistical un- 
certainty in the initial polarization measurement. 

We conclude that, using the available theoreti- 
cal estimates,'~* and barring accidental cancella- 
tion of the effects on the muon moment, the fol- 
lowing assertions can be made with 95% confi- 
dence: 

(1) Conventional quantum electrodynamics is 
applicable to distances as small as 7x107** cm. 

(2) The “radius” of the muon is less than 
4.5x107™ cm. 

(3) No fundamental length exists with magnitude 
larger than 2x107** cm. 

(4) The coupling constant G*/47 of the muon to 
an unknown field of nucleonic mass is less than 
~3x10°%, this limit varying approximately as 
(Ina)™, with A=M/my. 

Combining our result with the muon precession 
measurement,® we obtain for the mass of the 
muon 206.77+ 0.01 electron masses. 

We wish to thank Professor G. Bernardini in 
particular for several decisions which were 
essential to the successful conclusion of this ex- 
periment, and Professor C. M. York for his 
close collaboration. Mr. B. Nicolai was of great 
technical assistance, as were R. Bouclier and 
L. Magnani. Professor W. K. H. Panofsky con- 
tributed very significantly during the fall of 1959. 
We are grateful to several physicists in this la- 
boratory, in particular Dr. A. Petermann, for 
clarification of theoretical questions. 

A more complete report, with details of the 
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experimental techniques, of orbit theorems, etc., 
will be published. 
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Following the CERN Conference in 1958' there 
was considerable interest shown in the possibility 
of a peak at ~30-Mev K -meson laboratory en- 
ergy in the elastic scattering cross section of K™ 
mesons on protons. Matthews and Salam? pro- 
posed an interpretation of the peak in terms of a 
J=1/2 resonance of the K'-p system. However, 
Jackson and Wyld® pointed out that the peak could 
also be satisfactorily understood on the basis of 
S-wave zero-range K -p scattering theory and a 
repulsive K -p nuclear interaction. The reso- 
nance proposal has also been discounted by 
Dalitz and Tuan‘ who pointed out that the effective 
range of the K -p interaction would have to be 
unreasonably large in order to exhibit the ob- 
served elastic scattering cross sections at the 
low-energy values in the region of the so-called 
resonance. 

At the Kiev Conference in 1959, bubble cham- 
ber results® were presented which no longer ex- 
hibited a peak in the low-energy elastic K -p 
scattering cross section. However, further 
emulsion results® by the combined Bologna- 
Munich- Paris- Parma laboratories still tended 
to show the elastic cross-section flattening off 
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at the lowest energies. 

Emulsion data from the European K™ Collabo- 
ration groups’ have now been carefully scruti- 
nized for low-energy K -proton elastic scatter- 
ing events. All two-pronged events which could 
possibly be attributed to K -p elastic scattering 
have been analyzed. Low-energy elastic scatter- 
ing events on free protons can be readily recog- 
nized because these only will show coplanarity 
with complete balance of momentum and energy. 
Moreover, since the fraction of inelastic scat- 
tering events in which the K' meson re-emerges 
from emulsion nuclei is known® to be only a few 
percent for energies in the range 30-100 Mev, 
there is a negligibly small chance of an apparent 
K’ -p scattering occurring on an internal proton 
of an emulsion nucleus, especially for K” meson 
energies below 10 Mev. 

The results of our investigations, based on a 
study of 10850 K” mesons coming to rest and 
2060 K” mesons decaying or interacting in flight, 
are shown in Table I. Our results, which are 
shown with others of the bubble chamber® and 
emulsion® groups in Fig. 1, appear to indicate 
that the elastic K -p scattering cross section is 
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Table I. Summary of low-energy K -p elastic scattering events. 








p (Mev/c) 53 93 133 173 213 253 
Number of events 4 6 14 15 17 
Track lengths (cm) 659 1769 3977 6804 10497 
0 ,,(K -p) (mb) 187 105 108 68 50 





increasing in the region of lowest K -meson mo- 
mentum. Although the statistics are still quite 
poor, on the basis of our observations, there is 
only an approximate 10% probability that the 
K’-p elastic scattering cross-section curve flat- 
tens off at momenta lower than 120 Mev/c. 
There are possibly two explanations of the ob- 
served flattening in the earlier emulsion experi- 
ments. One could possibly be that because the 
track lengths in the lowest energy intervals were 
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FIG. 1. Total elastic K -p scattering cross sections 
versus K~-meson momentum. The curves, which 
were taken from the theory of Jackson and Wyld,’ show 
the Fol variation to be expected from constructive in- 
terference a+ and destructive interference a- between 
the Coulomb and nuclear parts of the K~-p scattering. 


so small, different experimental groups probably 
made energy cutoffs at different values. In our 
present investigation we have consistently worked 
down to a K-meson momentum of 53 Mev/c, for 
which the residual range is 100 microns. 

The other possible reason that many low-energy 
scattering events may not have been recorded is 
that they were classified as “double stars,” i.e., 
they were regarded as K-meson capture stars 
and secondary stars resulting from interactions 
or decays of hyperons or hyperfragments. Photo- 
graphs of two very low energy K -p scattering 
events are shown in Fig. 2. 

There is, of course, a third possibility: that 
the differences are largely statistical. 

At the very low K -meson energies considered 
in these experiments there may be, in point of 
fact, a significant contribution to the elastic scat- 
tering from the Coulomb field of the proton. The 
Coulomb effect on the total elastic scattering 
cross section has already been calculated, for 
the conditions of our experiment, by Jackson and 





FIG. 2. Two examples of low-energy K -p elastic 
scattering events. The energies of the incident K~ 
mesons are 5.5 and 5.3 Mev. Scale divisions (between 


photographs) are 10 microns apart. 
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Wyld.* They considered two cases, a+ anda-, 
of the nuclear K -p scattering phase shifts, and 
these two curves are shown in Fig. 1. It is seen 
that our results, as also those of the Berkeley 
bubble chamber group, suggest the attractive a+ 
set of K -p nuclear phase shifts. 





*Fulbright and Guggenheim scholar on leave from 
University of Illinois, Urbana, Illinois. 

Ton leave from University of Mysore, Mysore, 
India. 

ton leave from University of Warsaw, Warsaw, 
Poland. 

Miss V. Mayes and Mr. M. G. Bowler assisted in 
the early part of the work. 
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1958), p. 171. Both the Berkeley Bubble Chamber 
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= PRODUCTION BY K” MESONS 


William B. Fowler,* Robert W. Birge,t Philippe Eberhard,} Robert Ely, Myron L. Good, ! 
Wilson M. Powell, and Harold K. Ticho** 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received January 6, 1961) 


The Berkeley 30-inch propane chamber was 
operated in a 1.17-Bev/c K beam! at the Beva- 
tron to produce cascade particles. Out of the 
100 000 pictures taken, there have been 18 posi- 
tively identified negative cascade decays with the 
decay mode = +A+z. Of these, six were pro- 
duced with a K* meson and three were accom- 
panied by a visible 6, decay. Three of the events 
with a K* satisfy the kinematics for the hydrogen 
interaction 


K +p+= +K", 
and the other three are assumed to have been 
produced from protons in carbon. Two of the K* 
decayed in the chamber and the other four were 
identified by comparing momentum and ioniza- 


tion. Since only one-third of the K° produced in 
the interaction 


K™ 4+n-= +K°® 


decay into a visible mode, it is quite consistent 
to attribute the nine events in which no positive- 
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strangeness K particle was observed to the neu- 
tral mode of the 6, decay, to the 6, decay, or to 
a K* that could not be distinguished from a pro- 
ton. 

The 18 events were selected from a sample of 
47 in which a V particle was thought to come 
from a prong of aK interaction which appeared 
to decay into az. The selection consisted of a 
series of two single-origin constraints on an 
IBM-704 calculator: (a) the Q of the V particle 
was required to be consistent with that of a A, 
37 Mev, and its line of flight consistent with the 
assumed = decay origin; (b) the A and the 7” 
were required to be coplanar with the assumed 
= and to balance transverse momentum. Re- 
vised estimates of the momentum and angles of 
the A anda were obtained from the constraint 
program, and the average mass of the =" for the 
18 surviving events, weighted by the estimated 
probable error of each event, was found to be 

M,._- =1317.9+ 1.9 Mev. 
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FIG. 1. Bartlett’s 
S(t) function with a A 

lifetime correction for 
(a) the 18 observed =~ 








-~ 4 4 4 | 

from this experiment, = . 1 3 ? 4 
(b) this experiment cor- 1/r(* 10sec) 
rected for scanning bias, ofh~ 
and (c) the corrected - 
data of this experiment i 
combined with six events -2r- 
from Trilling and Neuge- 

: bauer. ® 

é -3r- 

The mean life of the = has been determined corrections for chamber size. The dashed line c 
by means of a modification of the statistical in Fig. 1 shows the S function for their six events 
method of Bartlett? which accounts for the neces- combined with our 25 (including corrections for 
sity of seeing the A decay before the = leaves scanning bias), and it is seen that the combined 
the chamber. The possible time for decay was mean life falls within our error. 


determined by calculating the time of flight to 

the boundary of the fiducial region, which was 100 
chosen to allow a minimum of 4 cm for the meas- | PT | | | 
urement of the A decay. In all cases the correc- 
tion due to finite chamber size is small, and the 
total correction amounts to less than 10% of the 
uncorrected mean of the actual flight times. The 
flight times include dE/dx corrections and the 
possible times were modified for those cases in 
which the = would have come to rest before 
leaving the chamber. Figure 1 shows a plot of 
Bartlett’s S function as a function of 1/7. This 
function is zero at the most probable lifetime 

and + 1 at the best estimate of the standard er- 
ror. Figure 2 is a cumulative histogram of the 

] calculated times of flight, ¢;, and indicates a 
considerable scanning bias against flight times 
less than 0.5x10-*° sec. An estimate of the total 
number of events was obtained by fitting those 
events with a time of flight greater than 0.5x10-*° 
sec to the exponential Ne~'/ T, The revised mean 
life consistent with this procedure required the 
addition of seven events in the interval from zero 
to 0.510" sec, and was found to be 





a 
aan f 


| 
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i) 


Number of events 




















T= (1.2879-42) x10"!° sec. ' . a 
The S function corrected for scanning bias is 0 05 lO 15 20 25 30 35 





shown as the solid line } in Fig. 1. This mean Time of flight,t; (10sec) 

life is considerably shorter than that found by FIG. 2. A cumulative histogram of the times of 
Trilling and Neugebauer.* However, their events flight of the =. The solid line is the best fit of those 
at high cascade momentum required much greater events with ¢;>0.5x 10-"* sec to an exponential. 
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A violation of parity conservation in the cas- 
cade decay will manifest itself as a fore-aft 
asymmetry in the A decay, the = serving to 
polarize the A with the strength a= which is then 
analyzed by the A decay. Denoting a unit vector 
in the direction of the A in the = rest frame as 
4 and a unit vector in the direction of the proton 
in the A center-of-mass system as /, the distri- 
bution of u-p is* 


P(a-p) d(u- b) = 3[1+ (cx,.a2 ,)(é- b) |d(u- 5). 


We have found, on the basis of the above 18 
events, 


(a-a,) = -0.65+ 0.35, 


using (aa, ) = (3/N)2.,(4-p); as the best method 
of estimation. The negative sign shows the 
helicity of the A from the = to be opposite to 
that of the proton from A decay. Birge and Fow- 
ler® find the helicity of the proton to be positive, 
implying that the helicity of the A is negative. 
Using the best estimate of the magnitude of |a A! 
as 0.95+ 0.07,° and the sign of a, as determined 
by Birge and Fowler, we find 


a,_.=0.69+ 0.36. 


a 
_ 
— 


The possibility of an up-down asymmetry of 
the decay 1” from = ~A+zn produced in the re- 
action 

K +N-= +K 


has been examined, and (a=P) is -0.28+ 0.40, 
where P is the polarization of the = in produc- 


tion and a— the decay asymmetry parameter. 
There seems to be no indication that the = pro- 
duced in this manner from carbon are polarized. 

In the 100000 pictures we had a total of (1.5+ 0.2) 
x10’ cm of K track length. After adding seven 
events to correct for scanning bias and then 50% 
to correct for the unobserved neutral decay mode 
of the A, we estimate the total number of = to 
be 37.5. On this basis the = production cross 
section in the reaction K +N-= +K is 


o,.- = 18+ 5 ub/nucleon, 


assuming an A”* correction for shielding in the 
carbon nucleus. 
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TEST FOR PION-PION RESONANCE 


G. Pinski 
Department of Physics and Astronomy, University of Rochester, Rochester, New York 
(Received January 9, 1961) 


It has been suggested’ that the two-pion effec- 
tive mass squared distribution in antinucleon 
annihilation could serve as an efficient test for 
the proposed pion-pion resonance. This current- 
ly fashionable resonance has been investigated 
on several fronts.*~™* In this note, we would like 





to suggest a particular decay mode for study, in 
which the effects of a pion-pion resonance would 
be most easily detected. 

The relative transition probability for annihila- 
tion into m mesons is given, according to the 
covariant statistical theory,® by 


n n 
w,~ J I letp, 9(0,00(? DY Py, Py dM Py + Py Ep). 
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In this formalism the invariant function f is taken 
to be independent of the nucleon and pion four- 
momenta. A pion-pion resonance may be incor- 
porated into the statistical theory by choosing f 
to be a resonant function in the effective mass of 


pairs of particles, peaked at the “resonant mass”: 





n 
nPsa(t= [IT lato, 0,290(0*-1)] fL+x° Dol p 


i>j 


in units of the pion mass |p, | =[(p;°)’ - |p;|7]*” =1. 
This quantity is as useful as the angular cor- 

relation, but being covariant, is incomparably 

easier to evaluate. Since the effective mass 

(p, +p)" is equal to 2+2p,-p,, it can be thought 

of as a measure of “covariant angle”: 


6 = Pi* be 
cov |p, 1 |pz2! 


This distribution has been calculated in the two 
extreme cases: 

(1) a=0. a” is the distribution for no resonant 
term. 

(2)A=. ,P’ is the distribution for resonant 
term only. 

For case 2 there are three contributions: 

(a) nP ry’ if the correlation is between the two 
particles which are in resonance. 

(b) nPyf' if the correlation is between one reso- 
nant and one nonresonant particle. 

(c) nPf ‘ if the correlation is between nonreso- 
nant particles. 

If one considers a pion-pion resonance in the 
T =1 state, there will be a decided difference be- 
tween the A=0 andA=«terms. In the T=1 state, 
like pions cannot be in resonance and the like- 


=p,* b2- 


pion correlation will not contain the term nPrr'(u?) 


which equals ,P°(u?)5(u?-M*). This term con- 
tains the effect of the resonance explicitly and is 
sharply peaked while the other terms are rela- 
tively broad. One therefore expects a marked 
difference between (++) or (--) and (+-) correla- 
tions. 

For unlike pions, the A= term is of the form 
ayPryy'+ bynPyf'+C¢ nPrf, while for like pions it 
is b’,»Pys'+C'nPr', the coefficients depending 
on the particular charge complexion of the anni- 
hilation. Only the ratios b/a, c/a, and c’/b’ are 
important here. Asn increases, }b and c increase 
relative to a, since there remains only one way 


f =1+072 F[(p.+p.)? -M’], 
“Uw 
where ( pj +pj)°=(p;° +p; - (Bj +,’ and 2°, M? 
are two constants with the dimensions of the 
square of a mass. The effective mass squared 
distribution, »P,,(u”), where J,,P,2(u7)d(u?)=W,, 
is given by 


n 
i+) -M*)}O[(p, +o) 7 yu] o“[P yt Py “2.041, 





for the pair being correlated to be the resonant 
pair, while the total number of ways of choosing 
the resonant pair increases. Therefore the effect 
of a pion-pion resonance is more prominent for 

n small. The annihilation B+n~-1~+2~ +7" 
should be ideal. Here the A= term for (--) is 
3Pyf’ and for (+-) is sPyy'+3Pyf’. 

In Fig. 1 these distributions are plotted along 
with the A =0 distribution, normalized to the same 
area against u. Most of the area under the (+-) 
curve is concentrated around the resonant mass 
so that the difference between (--) and (+-) curves 
is unmistakable. 

These distributions have been calculated for 
annihilation at rest, E =1.874 Bev or 13.4 pion 
masses, using M =4 as resonant mass. For reso- 
nant function the smeared-out 6 function, 


B(u? -M*)= 35 oe, 


was used, with half-width € taken as 1/2 pion 
mass. 

This result assumes that there are no 7°’s 
present in the annihilation. In approximately 20% 
of the observed® p+n events, the charged pions 
are (--+) but there is an estimated average of 
1.4 7°’s. If the experiment were performed in a 
deuterium bubble chamber, it would be easier to 
verify by momentum balance those events in 
which no 7° is present. 

The (+-) correlation is quite sensitive to the 
position of the resonance, and so should serve to 
fix the resonant mass. If sufficiently detailed ex- 
perimental information is available, it may be 
possible to determine the details of the resonance 
function. There has been speculation’ on the non- 
occurrence of the two-pion annihilation. If the 
annihilation proceeds through the resonant (A = ©) 
channel, a sharp resonance will reduce the phase 
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FIG. 1. Pion mass squared 
distributions for p+n— 17 
+n-+n"*, 
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space sufficiently to account for the suppression 
of this mode. 

Even if the number of such annihilations is 
small compared to the number of events with one 
and two neutral pions, the poorer statistics should 
be compensated for by the pronounced peaking of 
the effective mass, which would tend to be sup- 
pressed for higher multiplicities. 


The author would like to thank Professor E. C. G. 
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in this work, Dr. T. Yamanouchi and Dr. S. Iwao 

for useful discussions, and Mr. K. T. Mahantappa 
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TWO- PION EXCHANGE INTERACTION* 


G. Breit, K. E. Lassila, H. M. Ruppel, and M. H. Hull, Jr. 
Yale University, New Haven, Connecticut 
(Received January 12, 1961) 


There is evidence for the one-pion exchange 
(OPE) interaction in nucleon-nucleon scattering.~* 
The agreement of the coupling constant g with 
that from pion physics is satisfactory, and tests 
of the form of the OPE potential (OPEP) per- 
formed by adding a central field term and by 
varying proportions of spin-spin and tensor 
force terms have not indicated* a significant de- 
viation from expectation for angular momentum 
L>4. The inclusion of phase shift K, indicated* 
a statistically significant coefficient of a central 
force addition to OPEP and it has been stated* 
that this probably indicates the setting in of two- 
pion exchange (TPE) effects. The statement was 
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based on comparisons of Gupta’s TPE potential® 
(TPEP) with the OPEP at the classical turning 
point for L=4. Results of further comparisons 
and of employing dispersion relations® are re- 
ported below. 

Gupta’s potential has a direct relationship to 
Dyson’s scattering matrix S in terms of free- 
particle wave functions. In this respect the em- 
ployment of field theory is similar to that of dis- 
persion relations. The location of singularities 
does not enter the application of field theory, but 
unknown contributions’ of higher powers of g” 
interfere with applicability of theory except at 
high L. The numerical values® of the TPEP use 
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nonrelativistic approximation and do not include 
some additional effects.® 

Since the TPE as well as the sum of OPE and 
TPE effects on phase shifts in states with L>1 
is in most cases <0.2, accurate distinctions be- 
tween a phase shift 5, and sin6d or tan6 do not 
enter and effects of second and fourth order po- 
tentials V, and V, were added employing Taylor’s 
formula for each.* Effects of the order of a per- 
cent corresponding to 5=0.2 are of no interest 
here. The effect of V, was obtained with the aid 
of numerical values in Table I of reference 5. 
The integrand of the phase shift formula has a 
maximum close to the classical turning point, 
and for V, the relative importance of small 7 is 
higher than for V,. The effect of V, was calcu- 
lated as in reference 9. On account of the mean- 
ing of V,, questions of wave function distortion 
do not arise and the integral converges. 

According to Grashin,’° in triplet states with 
total angular momentum J =L+1 the energy de- 
pendence of the OPE phase shift is anomalously 
rapid at small energies, E. This may be seen 
as a cancellation of dominant terms employing 
the formulas in reference 9. These states are 
thus especially suitable for TPE detection. The 
‘Pp, phase parameter @, is higher than its OPE 
value for all recently published fits*’ which fall 
in a band narrow compared with the distance 
from OPE. For °F, only the rather improbable 
fit YRB3 crosses the OPE curve and falls below 





it above 180 Mev. All other fits give values of 
oF greater than those for OPE. Similarly for 
the n-p phase parameter” 9D, the OPE value is 
below all empirical fits and has the wrong sign. 
An exception to these relationships is e&, for 
which most fits give values below those of OPE. 
But one of the empirically probable fits falls 
above OPE at low E, crosses it at 180 Mev, and 
differs from the sum of OPE and TPE at 274 
Mev on the limit of error. 

In Table I the contributions of V, and V, are 
shown for T=1. For oF, V, definitely overcor- 
rects the difference between OPE and the fits, 
the standard deviation being small and all 
searches giving nearly the same values at the 
energies considered. Evidence below shows that 
this is probably due to short-distance effects. 
The disagreement of calculation with YLAM 
values of 6°, and oP is more marked and simi- 
larly in the case of K,. For K, the agreement of 
calculated and YLAM values is good at the two 
lower energies and is much better than between 
OPE and YLAM. At 274 Mev, V, overcorrects 
OPE, presumably because of short-distance ef- 
fects. For 0f, 2 the accuracy of YLAM and the 
agreement among searches are poor."' At the 
two higher energies calculation and YLAM dis- 
agree markedly, but on account of the poor ac- 
curacy, a conclusion is difficult. An additional 
uncertainty may enter on account of coupling of 
P, to F,. The disagreement is therefore not 


Table I. Comparison of combined effect of V,+V, (in radians) with fit YLAM (T=1). 








Phase parameter E (Mev) V2 V4 V,+V,-YLAM Standard deviation 
F 
0 > 68.5 0. 026 0.23 0.12 0.011 
137 0.05 0.48 0.30 0.009 
K, 68.5 0.025 0.018 0. 000 0.006 
137 0.035 0.067 0. 006 0.003 
274 0.044 0.176 0.055 0.007 
oF, 68.5 0. 0096 0.0014 0. 0040 0.009 
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Table I. 


Comparison of combined effect of V,+V, (in radians) with fit YLAN3M (T= 0). 








Phase parameter E (Mev) V2 V, V,+V,- YLAN3M Standard deviation 
D 
0; 68.5 -0.132 0.025 0. 066 0.011 
137 -0. 265 0.094 0.084 0.010 
274 -0.467 0.252 0.164 0.03 
D 
6 > 68.5 0.183 0.025 -0.07 0.018 
137 0.33 0.09 0.04 0.013 
D 
0 3 68.5 -0.021 0.025 -0.010 0.010 
137 -0. 050 0.094 0.010 0.006 
274 -0.10 0.25 0.07 0.022 
Ky 68.5 -0.031 0.002 0.002 
137 -0. 056 0.014 0.015 
274 -0. 073 0.059 0.045 0.030 





surprising. The F, search YLAM used the OPE 
value below 160 Mev. The disagreement at 137 
Mev is, therefore, not significant and bf was 
probably forced towards small values by the en- 
forcement of OPE. This view is supported by the 
almost exact agreement of the calculated value 
with fit YLA at 137 Mev and the location of YLA 
nearly halfway between YLAM and YRB1. Simi- 
larly for @¥,, the 137-Mev value falls on YRB2 
and overshoots YRB1 only slightly. Possibly in 
both cases YLAM has underestimated the values 
at 137 Mev. 

In Table II similar comparisons are made with 
fit YLAN3M to T=0 states.” Effects of V, are 
in the wrong direction to give agreement for *D,. 
Since °D, is coupled to °S, to which the applica- 
tion of V,+V, is inadequate, this disagreement is 
in the same class as that for oF. For 6D, the 
addition of the V, effects is helpful at 68.5 Mev 
and overshoots YLAN3M at 137 Mev. A number 
of other fits fall in between the calculated value 
and YLAN3M, however. There appears to be 
agreement on the whole in this case. For eD, 
the OPE value has the wrong sign but TPE brings 
about agreement at 68.5 and 137 Mev on the limit 
of error. At 274 Mev the TPE overshoots again. 
For K, the YLAN3M used OPE values in 0<E<175 
Mev. The calculated value exceeds YLAN3M by 
half the standard deviation but YLAN3M would 
be higher had OPE not been enforced below 175 
Mev. There is agreement on the whole and de- 
viations from V,+ V, follow similar patterns for 
T=0 and T=1. 
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Tables III and IV show comparisons of phase 
parameters for T=1 and T=0, respectively, ac- 
cording to the dispersion relations (DR) treat- 
ment of Galanin et al.® and of Grashin and Kob- 
sarev® employing their tabulated ratios of values 
of TPE to OPE. The TPE gives improvement 
for *P, and *P,, but not for *P,. The critical 
quantity (L/k)(m_c/h)=Lyu/k is only 1.5 at 100 
Mev and is hardly large enough to expect agree- 
ment. For K,, i.e., ‘D,, the potential gives by 
far the better agreement; for F, the DR are 
somewhat the better; for F, and F, the DR give 
the better fit to YLAM but the potential is con- 
ceivably the better taking YRB1, YLA into ac- 
count. For T=0 the DR are much the better for 
K,, i.e., \P,, and also for *D,, but the potential 
is the better below 200 Mev for °D, and much the 
better up to 300 Mev for *D,. For K, the DR 
agree best with YLAN3 but considering the en- 
forcement of OPE below 175 Mev for this fit the 
potential may be actually the better. For G,, 

G,, G, the two ways of treating TPE do not appear 
significantly different considering the uncertain- 
ties of the fits to data. 

In the calculations of the effect of V,, the maxi- 
mum of the integrand falls at x,, >1 (x = rm c/n) 
only for L 22 at the three energies used. At 274 
Mev, x,,>1 only for L>2.2, and cor L=2 this 
criterion is barely satisfied at 137 Mev. The 
disagreement of the potential with the empirical 
fits for P terms is thus understandable and so 





are the poorer values for L =3 at 274 Mev. The 
reasonableness of this situation combined with 
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Table III. Comparison of combined effect YLAM+A 
of one-pion and two-pion exchange effects (in radians) 
according to Galanin et al. and Grashin and Kobsarev,* 
with fit YLAM for T= 1. 








Table IV. Comparison of combined effect YLAN3M 
+4 of one-pion and two-pion exchange effects (in radi- 
ans) according to Galanin et al. and Grashin and Kob- 
sarev, * with fit YLAN3M for T= 0. 











State E (Mev) OP TP A 
SP, 10 -0.052 0.004 0.002 
40 -0.176 0.018 -0.016 
100 -0.341 0.051 -0.056 
3p, 10 0.002 0.002 -0.009 
40 0.014 0.028 -0.045 
100 0.040 0.080 -0.071 
'D, 40 0.016 0.0002 -0.007 
100 0.032 0.008 -0.033 
200 0.041 0.020 -0.071 
KF, 40 0.0042 0.0002 0.0013 
100 0.0160 0.0010 0.16 
5, 40 -0.094 0.002 0.002 
100 -0.031 0.0012 0.0012 
200 -0.059 0.005 0.0026 
5F, 40 0.0008 0.0002 0.0002 
100 0.0041 0.0016 0.0016 
200 0.0104 0.0073 0.0049 
300 0.0158 0.0158 -0.0008 





4See reference 6. 


the striking improvements for K,, D,, D, caused 
by V, indicate the reality of the TPE. This con- 
clusion is strengthened by the DR calculations. 
The disagreements of the potential with data at 
low L and the higher E indicate the presence of 
additional effects in the general region x <1 and 
closer examination shows the presence of spin- 
orbit like effects with an “inverted” order. 
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Energy Commission and by the Office of Ordnance Re- 
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INFLUENCE OF THE An RESONANCE ON CORRELATIONS IN THE REACTION K™ +d=A+7°+pt 


Orin I. Dahl, Nahmin Horwitz,* Donald H. Miller, Joseph J. Murray, and Paul G. White 
Lawrence Radiation Laboratory, University of California, Berkeley, California 
(Received January 9, 1961) 


When K” mesons are absorbed in deuterium, 
two possible reactions are 


K~ +d-A+m° +p, (1a) 
K~ +d-D°+n° +p 
~~ 
Arty. (1b) 


We have examined in detail a group of events of 
types (la) and (1b) obtained in an exposure of the 
15-in. deuterium bubble chamber to a 450-Mev/c 
separated K” beam.’ Application of an energy- 
momentum balance at the interaction vertex leads 
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FIG. 1. Distribution of proton and pion kinetic en- 
ergies in the reaction K”+d—~A+1~+p (282 events). 
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FIG. 2. Proton energy spectrum for the reaction 
K~+d-~A+n~+p. A: direct production. B: =-A 
conversion. C: resonance events (half-width +20 
Mev). 
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to an unambiguous separation of the A and £° 
events except for a small fraction of those cases 
in which the recoil proton is too short for meas- 
urement.” The fitting constraints also ensure a 
high degree of accuracy (~1%) in the determina- 
tion of the momenta of the reaction products. 
Figure 1 shows the distribution of the kinematic 
variables T, and T,- (proton and pion kinetic en- 
ergy, respectively) for 282 events which (a) were 
produced by absorption of stopped K’s as deter- 
mined by the fitting procedure and (b) did not re- 
sult from the decay of a £° [reaction (1b)|. In 
addition, the pion and proton spectra for these 
events are plotted separately in Figs. 2 and 3. 
Two prominent features of this absorption re- 
action are reflected in the pion peaks at T_- =92 
and 147 Mev. We consider the higher momentum 
peak first. Because of the loose structure of the 
deuteron, it is expected that a group of A’s will 
appear as a result of K' absorptions on single 
nucleons. This group of events may be described 
by an impulse model which includes the effects of 
interactions in the final state. In reaction (1a) 
the situation is particularly favorable: since the 
pion-nucleon system is produced in the weakly 
interacting isotopic spin J=1/2 state, only the 
Ap interaction need be taken into account. The 
appropriate modifications of the impulse model 
have been developed by several authors.*»* They 
point out that for stopped K~’s captured from 
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FIG. 3. Pion energy spectrum for the reaction 
K~+d—-A+n~+p. 
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atomic S orbitals,® the YN (hyperon- nucleon) sys- 
tem is produced predominantly with relative 
angular momentum L =0, so that interactions in 
other than S waves of the YN system may be 
neglected. Using the formulas given by Kotani 
and Ross,‘ we find that while the angular correla- 
tions in the final stage depend sensitively on the 
strength of the Ap interaction, the total produc- 
tion rate for given T, is only slightly affected. 

It is apparent from . 2 that this model quali- 
tatively accounts for absorptions with T,<10 
Mev; however, it predicts that less than 20% of 
the directly produced events will have Tp >10 
Mev. 

The lower energy peak occurs at pion energies 
expected for production of £’s in the K” absorp- 
tion reaction. Since a A rather than a > emerges 
from the absorption vertex, a two-step process 
is suggested: the K” interacts with one nucleon 
to form a = and 7, and then the = interacts with 
the other nucleon and converts to a A via the re- 
action 5+N-A+p.°**® For K capture from 
atomic S orbitals, it may be shown that the ob- 
served conversion rate is readily accounted for 
by reasonable choice of parameters in the zero- 
effective-range theory for the S-wave =N inter- 
action. A typical prediction of the theory for the 
conversion- pion spectrum is shown in Fig. 3.’ 
Since conversion occurs predominantly in the S 
wave of the ZN system, the angular distribution 
of the Ap relative momentum will be isotropic 
with respect to the direction of the recoil 7°. 

For such a distribution, the frequency of events 
along lines of constant T_- in Fig. 1 will be in- 
dependent of Tp- The calculated proton spectrum 
is shown in Fig. 2. Again, although the theory 
accounts qualitatively for the main features of 
the conversion process, there remains a group 
of events with T, >10 Mev and T_->100 Mev 
which are not readily associated with either di- 
rect absorption or internal conversion. 

Alston et al. have recently presented strong 
evidence for the existence of a resonance in the 
An system (hereafter called Y*) at a total mass 
of My«= 1380+ 5 Mev and with a half-width less 
than + 32 Mev.® More recent data support the 
existence of this resonance and suggest that the 
width may be as little as + 15 Mev.® Production 
of this state in the K°d absorption reaction would 
result in a peak in the proton spectrum in the re- 
gion near T, =30.5 Mev. Examination of Figs. 1 
and 2 indicates that the anomalous events may 
readily be attributed to production of the reso- 
nant Az system.*® Under this assumption, we 





have attempted to estimate the relative import- 
ance of the processes contributing to this reac- 
tion by folding together the three distributions in 
Fig. 2 appropriately normalized to reproduce the 
observed proton spectrum. Possible interference 
effects have been neglected. The resonance 
curve used corresponds to My*=1382 Mev, with 
a half-width of +20 Mev. From this we estimate 
that of the 282 events examined, roughly 87 are 
direct A’s, 102 are internal conversions, and 93 
are associated with the resonant channel." 

The spin of the resonant Az state has not yet 
been established definitely. In the absence of 
final-state interactions, the resonant channel 
will result in the An being produced in a pure 
angular-momentum state. For J=1/2, the distri- 
bution of the c.m. Az relative momentum will be 
spatially isotropic. For J=3/2, the distribution 
is not uniquely determined by conservation of 
angular momentum alone, but must still be sym- 
metric about a plane perpendicular to the direc- 
tion of the recoil proton.” However, in the re- 
gion (large Tp) where the effect of the Ap inter- 
action is likely to be small, the angular distribu- 
tion is distorted by overlap with the internal- 
conversion events. 

A region of particular interest is 5< Tp < 20 Mev, 
where both the Y* and direct production are im- 
portant, while the contribution from internal 
conversion is negligible. The observed events 
show a marked tendency for the Ap system to be 
produced with low relative momentum (large 
T,,- on Fig. 1). In particular, for the 48 events 
in this group we obtain 


W, N,)/ W, 
where Ny (Np) is the number of events for which 
the c.m. angle of the Az system is less than 
(greater than) 90 deg with respect to the recoil 
proton direction. If the quantum numbers of the 
Y* are the same as those of the zero- energy 

K~p system (S” for odd KA parity, P” for even), 
so that this is the resonance suggested by Dalitz 
and Tuan,’* the distribution of events for Tp <20 
Mev is readily understood in terms of an attrac- 
tive S-wave Ap interaction. In using the final- 
state interaction theory, the Am resonance was 
taken into account through an enhancement of the 
single-nucleon transition operator, (K"n|T\|Az"), 
as the energy of the An system was decreased 
below the Kn threshold. For a Ap potential of 
Gaussian form and range corresponding to two- 
pion exchange, the data suggest that the volume 


+N,)=0.44+0.18, 
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integral of the average potential is 280+ 90 Mev- 
fermi’. A strong S-wave Ap interaction could 
also be present if the resonant state were P*? 
and the KA parity even. IfJ is one-half, but the 
parity of the Y* differs from that of the zero- 
energy K~p system, or if the Y* state is P** and 
the K°A parity is odd, the recoil proton must be 
produced in an odd angular-momentum state with 
respect to the Y*. Since the effect of S-wave Ap 
scattering will then be small, the angular distri- 
bution can probably be accounted for by an inter- 
ference between the two production modes. In 
this case, the A will in general be polarized with 
respect to the production plane; however, the 
data are statistically inadequate to determine 
whether a significant effect is present. 

It is of interest to compare the rate for Y* 
production with the rate for the nonmesonic pro- 
cesses 


K” +d-D" +p; D°+n; A+n, 


which occur with a combined frequency of ~1.2% 
when stopped K~’s are absorbed."® Since reac- 
tion (la) constitutes 21.5% of all zero-energy 
K~ absorption events,’* we estimate that ~6% of 
the absorptions proceed through the resonant An™ 
channel. By charge independence, an additional 
3% proceeds through the An° channel. The sup- 
pression of nonmesonic absorption reflects the 
much smaller volume of the deuteron in which 
these events may occur because of the large mo- 
mentum (> 500 Mev/c) present in the final state. 
It is a pleasure to thank Professor Luis Alva- 
rez for his advice and encouragement throughout 
this work, and the assistance of the bubble-cham- 
ber and scanning staffs of the Lawrence Radia- 
tion Laboratory. 
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K-1 RESONANCE AND LIFETIME OF THE K’ PARTICLE* 


Mirza A. Baqi Bég and Paul C. De Celles 
Brookhaven National Laboratory, Upton, New York 
(Received January 12, 1961) 


Recent experimental data’ seem to indicate a 
resonance in the 7 +K system at an energy of 878 
Mev and full width at half maximum of 234+ 20% 
Mev. The spin and isobaric spin of the resonance 
have not yet been determined’; however, it is 
certainly attractive to consider the possibility 


that it might be due to the conjectured K’ particle,® 


a scalar isobaric doublet. The theoretical rea- 
sons for postulating this particle have been dis- 
cussed recently by Gell-Mann,° Tiomno,* and 
Bernstein and Weinberg.*® On the other hand, it 
has been suggested® that the small width of the 
experimental 7K resonance might be an indica- 
tion that it is a P-wave resonance’ (possibly due 
to a vector particle, say, Ky). 

In this Letter we wish to point out that the cross 
section for K’ (or Ky) production in K-p scatter- 
ing can be simply related to its lifetime and spin. 
The relations we derive may then be used to 
distinguish between the possibilities K’ and Ky, 
and provide a method for measuring the lifetime 
of the particle involved.® 

For the scalar case we have, on using Lorentz 
invariance and parity conservation, 


(4k ky (Kd 1K’) =A(9), (1) 


where J, is the pion current, q=k-k’ is the dif- 
ference of K, K’ 4-momenta, and As(q*) is a scalar 
function of ¢g’. The lifetime is simply® 
2 
- 87 Mi 
"s°¥(m") a’ 
1 





(2) 


where, with an cbvious notation for particle 
masses, 
_ 1 2 2 2\2 2 21,2 
°F OK? ~My -m_ ) 4m my | . (3) 
The cross section for K+p~-K’+p, when calcu- 
lated in the one-pion approximation in the bary- 
centric frame, is 
do\ fp?) ikl . 
oa) => —— » (4) 


64n°E" |k!| | -m | 


where E is the total energy of the system, and 


if (q?)N'y,N =(N’ \dq IN) (€€’/m?)™, @ =(k'-k)’, 





with k’ and k the 4-momenta of K’ and K, respec- 
tively, and with € and ¢€’ the initial and final pro- 
ton energies, respectively. Providing g* is such 
that A(q*) =a(m,”) and f(q*) =f(m,”)=G, the re- 
normalized pion-nucleon coupling constant 

(G?/4n =14.5), then 


1 2 ike oti 
my, Ik’ | g 


do G? 
(za), - (ae) are Tk Wa?=m 7. 
" (5) 


Using the experimental values of mx’ =878 Mev 
and £ corresponding to the data of reference 1, 
one finds for the total cross section 


o(mb) =1.9/(7, x10), (7, in sec). (6) 


For the vector case, the analog of Eq. (1) is 
v2 = . 
(Aege y) (Kld_ IK) =AL (a? E-k), (7) 


where ¢ is the “polarization” vector of Ky and 
we have used the usual subsidiary condition for 
a Proca field, viz., ky-&=0. Ay(q’) is the ap- 
propriate form factor. The lifetime of the Ky is 
then 
8n 3m KP 
vo z > 
V Ay (m_ ) @ 


(8) 





where we assume that mx,,=mx: =878 Mev. 
Corresponding to Eq. (5), we now have for the 
differential cross section 


(=), 


The total cross section, corresponding to 
Eq. (6), is now 


o(mb) =24/(7 1x10), (r, in sec). (10) 





Gy 3 1. Ik -¢ 
-(5) 270° EE? ik Ig-m 70 
x al4m ,?m,,? -(¢ Me, -m,°)'}. 
(9) 


The angular distributions of the recoil protons in 
the c.m. frame is plotted in Fig. 1. A measure- 
ment of the backward to forward ratio possibly 
could decide whether the resonance is to be 
identified with K’ or Ky. Although the angular 
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FIG. 1. Angular distribution of recoil protons, in 
barycentric frame, for the process K+p—K’'(Ky)+p. 
Both curves are normalized to the same total cross 
section of 1 mb. 


distributions plotted in Fig. 1 are clearly dis- 
tinguishable, these theoretical cross sections 
are due to single-pion-exchange graphs only. It 
is improbable that such terms dominate the cross 
section over the entire range of cos@. A pro- 
nounced backward peaking of the protons, how- 
ever, would be strong evidence for the K’, but 
lack of it would not rule out the same possibility. 
Alternatively, a precise knowledge of o and the 
resonance width, I, could distinguish between 
the two alternatives. With the available data, 
o =(1+0.42) mb and thus 


1.34< T X10 <3.28, 
16.9<1_x10%<41.4, (11) 


whereas 2.38 <(1/T'exp) x 107°<3.58. Therefore, 
present experimental evidence supports the K’ 
alternative. However, it has been said that the 
present experimental width Texp could be entire- 
ly due to experimental uncertainties. Hence, 

one cannot yet exclude the possibility that the 





resonance is due to a vector particle. 

It should be remarked that the experimental 
values of E and ¢ used above are not optimal for 
this kind of analysis. The use of the one-meson 
approximation in Eqs. (5) and (9) is best when the 
momentum transfer to the proton is smallest. 
The ¢ in the data of reference 1 is -0.12 (Bev/c)? 
for protons scattered in the forward direction, 
which is to be compared with m,” = +0.0196 (Bev)’. 
One should therefore consider the inequalities 
(11) as rather crude estimates until more data 
corresponding to smaller ¢* and larger E are 
taken. 

We wish to thank J. Bernstein and R. Marr for 
useful discussions of this work and I. Cole for 
assistance with the computations. 









*Work performed under the auspices of the U. S. 
Atomic Energy Commission. : 
1M. Alston, L. W. Alvarez, P. Eberhard, M. L. j 
Good, W. Graziano, H. K. Ticho, and S. G. Wojcicki, , 
Phys. Rev. Letters 5, 11, 520 (1960); also M. L. 


Good [private communication to one of the authors ‘ 
(P.C.D.)]. : 
?From parity considerations one rules out a pseudo- i 


scalar or pseudovector particle as the possible cause 
of the resonance. E 

3M. Gell-Mann named such a particle K’ |Proceed- B. 
ings of the Tenth Annual International Rochester Con- 
ference on High-Energy Physics, 1960 (Interscience 
Publishers, New York, 1960), p. 509]. 

4{J. Tiomno, reference 3, pp. 466, 467, and 513. 

5J. Bernstein and S. Weinberg, Phys. Rev. Letters 
5, 10, 481 (1960). 

®M. Gell-Mann (private communication to J. Bern- 
stein). 

™Previous theoretical calculations using dispersion 
relation techniques have indicated a possible P-wave 
resonance in the 7K scattering. B. W. Lee, Phys. 
Rev. 120, 325 (1960); and M. Gourdin, Y. Noirot, 
and Ph. Salin, Nuovo cimento 18, 651 (1960). These 
analyses cannot rule out “kinematical” S-wave reso- 
nance [a phrase used by G. F. Chew and S. Mandelstam, 
Phys. Rev. 119, 467 (1960)]. For our present pur- 
pose, it makes no difference whether or not we asso- 
ciate a particle with the observed resonance as long 
as the resonance exists in a pure state of spin and iso- 
baric spin. 

SEssentially similar techniques were employed by 
H. Primakoff [Phys. Rev. 81, 899 (1951)] in his dis- 
cussion of the r° lifetime. 

*This is similar to what one obtains in lowest order 
perturbation theory from a three-field interaction 
Lint =A"KK’, but questions of renormalization and the 
inapplicability of perturbation theory do not arise in 
our approach. 
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NONC LASSICAL STRUCTURE OF THE ENERGY- MOMENTUM TENSOR OF A POINT 


MASS SOURCE FOR THE SCHWARZSCHILD FIELD* 


Frank R. Tangherlini 
Institute of Field Physics, University of North Carolina, Chapel Hill, North Carolina 
(Received November 28, 1960; revised manuscript received December 22, 1960) 


Although the Schwarzschild solution to the field 
equations has been known for many years, be- 
cause of the complicated nonlinear character of 
the equations, there is still no generally accepted 
form for the structure of the energy-momentum 
tensor T,," when the latter is that of a point mass 
source. From a classical standpoint, one would 
surmise that 7," in a static coordinate system 
reduces to T,°, with 7,° proportional to a delta 
function, but from the structure of the field 
equations it is by no means obvious that other 
components are not also present. 

Indeed, as we shall now show, and briefly dis- 
cuss, this tensor does not have only one non- 
vanishing component, i.e., 7T,° (as one might ex- 
pect classically for a point particle at rest), but 
four: T,°=T,', T,2=T,°=-2T7,°, with T,°=m5(r)/ 
2nr*; the notation is described below. 

In order to arrive at this result, we shall first 
show that, for systems with static spherical 
symmetry, when a certain condition is satisfied 
by the g, ,, (as is the case for the Schwarzschild 
solution}, the field equations reduce to a pair of 
dependent linear equations, most simply ex- 
pressed in terms of the scalar potential U (de- 
fined by g,,=1+2U). Using these equations, the 
above structure for 7,” may be immediately in- 
ferred. 

Upon introducing the coordinate system asso- 
ciated with the following standard line element 
for systems with static spherical symmetry, 


ds? = g,,(r) dt® + g,,(r) dr? - r7(de?+sin?6 dg”), (1) 


one finds for the nonvanishing components of the 
field equations’: 


G,= -r-*[1+ ve") , = “KT ®, 


6,’ s 6,"* me r a8" 800, ? ‘ anal 
Go=-2™e EV" g5.)-[E+ ar(ingy,) ,](G,'- G,”) 


= KT,’, (2) 
and by spherical symmetry, G,*=G,?. We note 
the following about the equations: (a) The first 
equation is already linear in g™; (b) the other 





two equations would be linear if gg"" were a con- 
stant,® since this is equivalent to g" o ST te 
=0; by the asymptotic condition‘ this constant ” 
may be set equal to -1; (c) if g°g"*=-1, the trans- 
verse equation reduces to a Newtonian one, re- 
lating V*g,, to the transverse stresses; (d) in 
order for gg" to be constant, it is necessary 
and sufficient that 7,°=7,'. When the latter is 
the case, the requirement that T,’. m =0 reduces 
to 

wey 2rT,*=0, (3) 
and it follows that there is only one independent 
component for the Fe which most conveniently 
may be replaced by the trace, T. One then has 
in terms of T for the components 


T,°=T,=r-* [Tr dr+C/r, 


T,2=T,°=3T - T,°, . (4) 


where C is a constant of integration.°® 

Let us now consider the Schwarzschild field, 
i.€., Zy9=1-2Gm/r, g=-g5. This field is de- 
rived under the assumption that outside the 
source T,,"=0; hence in this region, the require- 
ments on Ag for linearity are satisfied. We 
shall continue the solution down to the source 
point at r=0, noting that g,,(=- g'") is finite and 
differentiable even when passing through the 
Schwarzschild singularity, and the regularity of 
the field equations is not altered by the fact that 
&o0= 9 (or g,,=%) at r=2Gm, since this is not a 
singularity in the differential equations them- 
selves. ® 

Upon introducing the scalar potential, -g''=g,, 
=1+2U, the field equations reduce to the follow- 
ing pair of dependent linear equations’: 


« -2 = ) 1 
2r (rv) , KT) ( KT, ), 
-vzy- 2(- 3 
v?U KT, ( aT, ), (5) 
where, for the case of a point particle, the source 
terms are singular distributions which vanish for 


y>0O. For the Schwarzschild solution, i.e., 
U=-Gm/r, one readily obtains the following 
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values for the he which may be conveniently 
summarized as 
T + diag (1, 1, -4, -4) m(r)/2ar" (6) 
v ag >") we & r mY , 
where 6(r) is the radial delta function i d(r) dr 
=4. Since r6’(r)=- 5(r), these distributions 
satisfy the requirement (3). 

Although the m appearing in the above equations 
is up to this point of the discussion to be thought 
of only as a gravitational mass, it follows di- 
rectly from the action principle that it is also a 
dynamical quantity, the rest mass. We have for 
the gravitational action, A, the expression 
J(-g)“*«x""R d*« which by the field equations may 
be written {(-g)“*Td*x. Hence, since T=T7,° 
=m6(r)/2nr?, the action integrated for a finite 
time interval over any spatial region w contain- 
ing the particle “ 


A(t; w) = [ af moe) dr =m(t-t). (7) 


The surfaces of constant action are therefore 
parallel to the surfaces of constant time, and 
for this case, and in general for static coordi- 
nate systems,® the normal derivative to the ac- 
tion surface is the rest mass, i.e., dA/dt=m 
= {(-g)T d*x.° 

The result is therefore quite analogous to spe- 
cial relativity because in the above case T,'+ T,? 
+T,°=0, so that 7,°=T as in special relativity 
for a particle at rest. However, despite this 
analogy, we encounter a difficulty if we try to 
describe Pe classically. We would expect 

= p(dx/ds)(dx,/ds), where p is the trace or 

proper density m6(r)/2nr?._ However, the pre- 
viously found values for T would imply, for- 
mally, in the limit r=0, 


Boo(dt/ds)* = 1, &,,(dr/ds)* =1, rv? (d6/ds )? = 3, 
y® sin?6(d¢/ds)* =4, (8) 


instead of as in special relativity (dt/ds)*? =1, 
with the other terms vanishing. Moreover, again 
formally, we have (since |g,,.|==, g,,=0, at 
v=0) dt/ds=0, |\dr/ds|=«, in contrast with the 
classical result dt/ds =1, dr/ds =0. The two 
angular terms may be formally equated to yield 
d6? = sin*6 dg”, which has the particular integral, 
sin@ exp| @| =constant, corresponding to spirals 
on the unit sphere.*° 

While these results follow because we are no 
longer working in the classical region of eynee- 
time, and because of the structure of a , it 
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would appear that we have reached a fundamental 
limitation here in understanding the point par- 
ticle in general relativity as a purely classical 
object (i.e., a structureless mass point, with the 
spatial stresses vanishing).** We leave as an 
open question whether this nonclassical structure 
has some relation to the nonclassical behavior 
of particles as exhibited in the quantum domain. 
In conclusion, to clarify the physical principle 
underlying the condition on the T,,", we note that 
one may arrive at this condition in the following 
way: Let a particle be freely falling radially in 
a gravitational field described by the line ele- 
ment (1); what is the condition on g,, and g,, that 
will make the radial acceleration independent of 
the energy integral E =g,,dt/ds ? One readily 
finds that the condition is g,,g,,=constant, and 
upon examining the field equations, one finds 
T,°=T,!. In other words, we are dealing here 
with a consequence of the principle of equivalence 
as it manifests itself in the field equations.” It 
is certainly remarkable that this condition should 
also make the equations, for this case, linear 
and, upon transforming to rectangular coordi- 
nates, make the coordinate system “proper” (to 
within a normalization constant). We hope to 
discuss these results further in a subsequent 
publication. 





*This research was supported by the United States 
Navy under a contract monitored by the Office of Naval 
Research. 

‘Because of the final linearity in go) and g"', it will 
not be convenient to make the customary substitution 
£9 = expu(7), &1,= -exprv(y). However, if one adopts 
static cylindrical coordinates, one obtains an equation 
linear in Ingo, as first shown by H. Weyl, Ann. 
Physik 54, 117 (1917); 59, 185 (1919). See also T. 
Levi-Civita, Rend. Accad. nazl. Lincei, (1917-1919); 
R. Bach and H. Weyl, Math. Zeit. 13, 134 (1922). We 
should like to emphasize the restricted character of 
our results, i.e., to systems with static spherical 
symmetry, and to call attention to the more general 
problem of finding rigorous linearizing conditions for 
the field equations, and determining, by experiment, 
whether or not they do occur. 

*We set c=1, x=8nG, where G is the Newtonian gra- 
vitational constant, also f vr =af/ar, Ve=r- {rf ar 
'This condition is not, however, necessary for lin~ 
oe €.g-, &o9= constant yields equations linear in 

, although it, of course, requires a different struc- 
ie for the T,, # than the Schwarzschild oe 

‘That is, pa ee @, Buy =diag (1, -1, -1, -1), 
to the extent that we may ap pect pre 

‘The term C/r*, with C =e?/2, where e is the charge 
of the electron in Heaviside units, occurs in the 
Reissner-Nordstrom solution, for which the T,# also 
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satisfy the above linearity condition, as does also the 
cosmological term Aé,". 

SMoreover, we note that the Schwarzschild solution 
for go is linear in m; hence we could alternatively 
reflect m and obtain a solution for which go) > 1 for 
0<r<@, One should therefore draw a clear distinction 
between singular points or regions in the field equa- 
tions, and singularities or “unphysical” behavior in 
ds?=g),,dx"dx” due to the values of the g,,,,. 

"Note that upon substituting the eugrenolens for T,4 
in (5) into the requirement (3), it is satisfied identi- 
cally, so that a reduced version of the Bianchi identi- 
ties has been preserved. 

5If one adopts the isotropic coordinates by setting 
r=(1+ Gm/27)*7, one finds the interesting result that 
the particle cannot be located at 7=0 in the isotropic 
system, and that rather 7=- Gm/2 corresponds to the 
particle’s location. Moreover, the isotropic system 
for the range 0<7<~ covers the space outside the 
Schwarzschild radius twice, so that 7=0 actually cor- 
responds to spatial infinity, as may be seen from the 
above transformation or by observing that the isotropic 
line element is form-invariant under the substitution 


q= G’m?/47¥. Compare the results, different from ours, 


obtained by R. Arnowitt, S. Deser, and C. W. Misner, 
Phys. Rev. Letters 4, 375 (1960), who assume the 
particle to be located at 7=0. 

SWe have also calculated the surfaces of constant 





action for Lemaitre’s nonstatic coordinate system 

[G. E. Lemaitre, Ann. soc. sci. Bruxelles, Ser. A. 
53, 51(1933)], and find A=m(t’+ r’)/2. The original 
static observers satisfy t’ -7’ = constant. Differentia- 
ting A in the direction dt’ =dr’, we have again m. 

An alternative classical representation for the 7,4 
is given by 7,4 = (p +p) (dx /ds) (dx,/ds) -p5,", where 
pis ascalar pressure.” One then has p=p/2 and 
again, formally, go(dt/ds)*= g,,(dr/ds)*=1, with the 
angular “motion” vanishing. 

‘One is perhaps reminded here of Einstein’s feeling 
that one should not assume a priori a classical repre- 
sentation for the singularities of the gravitational field 
in connection with deriving the equations of motion 
from the field equations. See the interesting accounts 
in L. Infeld, Helv. Phys. Acta. Suppl. IV, 206 (1956); 
Revs. Modern Phys. 29, 398 (1957). The above re- 
sults bring out the urgency of finding exact solutions 
to the two-body problem. 

12The question of the relationship of this condition to 
the principle of equivalence has been raised recently 
by H. Bondi and S. Kilmister, Am. J. Phys. 28, 508 
(1960), in their Letter commenting on the recent ap- 
plication of the principle by L. I. Schiff, Am. J. Phys. 
28, 340 (1960). Note, however, that Schwarzschild’s 
interior solutions do not in general satisfy this condi- 
tion, in contrast with the interior solutions that would 
be obtained from (5). 
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ABSTRACTS 





In this section are printed the abstracts of Articles 
that have been forwarded to The American Institute of 
Physics for publication in THE PHYSICAL REVIEW. 
In quoting information obtained from this section be- 
fore the appearance of the corresponding Article, ref- 
erence should be made to “Physical Review (to be 
published)” rather than to this Journal. 


DYNAMICS OF ATOMIC MOTIONS IN LIQUIDS 
AND COLD NEUTRON SCATTERING. A Rahman, 
K. S. Singwi, and A. SjSlander, Argonne National 
Laboratory, Argonne, Illinois (Received No- 
vember 16, 1960). 


The self-correlation function for the motion of 
an atom in a liquid has been constructed ina 


manner which simulates the behavior of this func- 


tion between the two extreme cases: viz., that in 
a solid and that in a simple diffusive model. Us- 
ing this function, the differential scattering cross 
section for cold neutrons has been calculated for 
liquid lead and compared with available experi- 
mental results. It appears that even in a simple 
liquid like lead the solid-like behavior of atomic 
motions persists for a considerable time. 


SOLUTION FOR THE TWO-ELECTRON CORRE- 
LATION FUNCTION IN A PLASMA. Burton D. 
Fried, Space Technology Laboratories, Incor- 
porated, Los Angeles, California, and H. W. 
Wyld, Jr., University of Illinois, Urbana, Illi- 
nois (Received November 23, 1960). 


The two-electron correlation function, g, re- 
sponsible for “collisional” corrections to the 
correlationless (or Vlasov) description of a 
plasma, is investigated. It is shown that an ex- 
act solution of the integral equation for g can be 
found for a fairly wide class of spatially homo- 
geneous, one-electron distribution functions, f 
(the ion dynamics being neglected). This is 
carried out in detail for the simplest member of 
the class (the resonance shape), and the Landau 
damping of g to its asymptotic (t~) form is ex- 
hibited explicitly. It is shown that correlations 
between particles separated by more than the 
Debye length are damped in a time which exceeds 
the period of plasma oscillations, w,~', and that 
these make an appreciable contribution to the 
“collisional” rate of change of f. It is concluded 
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that for rapidly varying f (as in problems in- 
volving plasma oscillations) conventional treat- 
ments of the “collision” term should be replaced 
by a self-consistent solution of the coupled equa- 
tions for f and g. 


FINE STRUCTURE AND MAGNETO-OPTIC EF- 
FECTS IN THE EXCITON SPECTRUM OF CAD- 
MIUM SULFIDE. J. J. Hopfield,* Bell Telephone 
Laboratories, Murray Hill, New Jersey and 
Laboratoire de Physique, Ecole Normale Su- 
périeure, Paris, France, and D. G. Thomas, 
Bell Telephone Laboratories, Murray Hill, New 
Jersey (Received November 15, 1960). 


The valence band of cadmium sulfide is split 
by spin-orbit and crystal field effects into three 
nearly degenerate bands at k=0. The magneto- 
optic absorption spectrum of direct excitons 
formed from the top valence band and the con- 
duction band has been studied in detail. Most of 
the experiments reported have been performed 
in light polarized parallel to the hexagonal axis. 
In this geometry, the exciton series consists of 
weak lines amenable to magneto-optic experi- 
ments. When the magnetic field and the wave 
vector of the light are perpendicular to each 
other and to the hexagonal axis, the reversal of 
the magnetic field produces large changes in the 
absorption spectrum. This effect can be quanti- 


' tatively understood as an interference effect be- 


tween allowed exciton transitions (optical matrix 
elements independent of the wave vector of the 
light) and forbidden exciton transitions (optical 
matrix elements proportional to the wave vector 
of the light). It is shown that in CdS the for- 
bidden processes having a principal quantum 
number 2 are somewhat stronger than allowed 
processes of the same quantum number. By 
using group theory and the effective-mass ap- 
proximation, the electron and hole anisotropic 

g values and masses are determined from an 
analysis of the exciton spectrum. The electron 
mass, 0.20.m (almost isotropic), determined in 
this analysis is compatible with the assumption 
that the k=0 conduction band valley is the lowest 
conduction band valley. The hole masses for the 
top valence band are my» , = 0.7m and Mp\ = Sm. 
An experimental upper limit on the slope of the 
conduction band at K=0 is obtained. 


*Now at the Ecole Normale Supérieure, Laboratoire 
de Physique, Paris, France. 
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LORENTZ FORCE ON SCREW DISLOCATIONS 
AND RELATED PROBLEMS. Jens Lothe,* Me- 
tals Research Laboratory, Carnegie Institute of 
Technology, Pittsburgh, Pennsylvania (Received 
October 14, 1960; revised manuscript received 
December 15, 1960). 


The concept of a Lorentz force on screw dis- 
locations, first introduced by Nabarro, is anal- 
yzed. It is concluded that the Nabarro Lorentz 
force on screws and the Lorentz force of electro- 
magnetism are not analogous, and that the term 
“a Lorentz force on screws” should be dropped 
in order to avoid confusion. Only when the screw 
is constrained to move on one slip plane is the 
analogy with electromagnetism complete. Total 
quasi-momentum is not generally conserved 
when screw dislocations interact with elastic 
waves. 


*On leave from Universitetets Fysiske Institutt, 
Blindern, Oslo, Norway. 


PHOTOCHEMICALLY PRODUCED COLOR CEN- 
TERS IN KCl AND KBr. John A. Cape,* Depart- 
ment of Physics, University of Illinois, Urbana, 
Illinois (Received September 26, 1960). 


KCl and KBr crystals were exposed to unfiltered 
mercury arc radiation at 15°K. If an “OH” band 
was present in the crystals before irradiation, 
the ultraviolet irradiation produced an optical 
absorption spectrum similar to that produced by 
x rays at 15°K. The optical absorption spectrum 
and the changes in the absorption produced by 
annealing at higher temperatures were measured 
and compared with the spectra observed during 
similar annealing of x-irradiated crystals. In 
KCl the 335-my band, formed by the uv irradia- 
tion, bleaches thermally at 56°K as does the 
335-my H band of x-irradiated KCl. In KBr the 
381-my band bleaches in steps at 35°, 46°, 56°, 
and 80°K as compared with the 381-my H band 
in x-irradiated KBr which bleached at 30°, 46°, 
56°, and 80°K. In both KCl and KBr the V, band 
appears with the disappearance of the H band. 
ilumination in the V, band causes regeneration 
of the H band, as occurs in x-irradiated KCl and 
KBr. In both KCl and KBr the photoproduced H 
band grows by about 10% at approximately 25°K. 
It is concluded that the H centers bleach ther- 
mally by diffusing to and combining with other 
color centers. Recombination with F and a cen- 
ters annihilates the H centers and leads to the 
formation of Vx centers, while recombination 


with a third center (possibly a positive-ion va- 
cancy) results in the formation of V, centers. 

U, U,, U,, O°, and @ centers are produced by 
the uv irradiation as well as F and H centers. 
The photoproduced F band may be bleached opti- 
cally with negligible effect on the H band. 


*Present address: Atomics International, Division 
of North American Aviation, Inc., Canoga Park, Cali- 
fornia. 


ENERGY RELEASE IN REACTOR-IRRADIATED 
COPPER. II. 600 TO 700°K RELEASE. T. H. 
Blewitt and S. T. Sekula, Solid State Division, 
Oak Ridge National Laboratory, Oak Ridge, 
Tennessee and J. Diehl, Max-Planck- Institut 
fur Metallforschung, Stuttgart, Germany (Re- 
ceived October 17, 1960). 


The energy release associated with the re- 
covery peak occurring between 600 and 700°K in 
neutron-irradiated copper was measured uti- 
lizing a new technique, that of nuclear heating. 
Following a bombardment at 40°C of 1.7x10”° 
fast neutrons of a 1/E distribution which raised 
the critical shear stress to 12.8 kg/mm? at 
4.2°K (5.2 kg/mm? at 300°K), a release of 7.7 
cal/mole was measured. Using this measured 
value of the energy release, it is possible to 
estimate the number of defects annihilated if it 
is assumed that the annealing is the result of 
the migration and subsequent annihilation of a 
single defect. In this way the number of inter- 
stitials, vacancies, interstitial-vacancy pairs, 
and dislocation lines required to account for the 
measured energy release were estimated. The 
values were, respectively, 5x10'* per mole, 
2x10”° per mole, 4x10'* per mole, and 1x10" 
cm per mole. 


FERROMAGNETIC RESONANCE LINEWIDTH IN 
COBALT-SUBSTITUTED FERRITES. C. Warren 
Haas and Herbert B. Callen, Department of 
Physics, University of Pennsylvania, Philadel- 
phia, Pennsylvania (Received November 14, 1960). 


The effect of small concentrations of substitu- 
tional ions on the ferromagnetic resonance line- 
width of a host ferrite is calculated. The reso- 
nance linewidth arises from magnon scattering 
induced by the variation from ion-to-ion of the 
spin-orbit interaction. This interaction is unique- 
ly large for cobalt ions because of the orbital 
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degeneracy of the ground state of the ion ina 
trigonally-symmetric crystalline field. The re- 
sultant contribution to the linewidth is found to 
be isotropic and of the order of 20-30 oersteds 
for each percent of cobalt in normal (i.e., non- 
inverse) ferrites or in ferrous ferrite. In other 
ferrites the effect is diminished by the lifting of 
the ground-state orbital degeneracy at some of 
the cobalt sites; this diminution is calculated as 
a function of the degree of inversion. 


DIRECT CURRENT ELECTROLUMINESCENCE 
AT LOW VOLTAGES. W. A. Thornton, Lamp 
Division, Westinghouse Electric Corporation, 
Bloomfield, New Jersey (Received November 
29, 1960). 


Electroluminescence due to dc excitation oc- 
curs in activated ZnS films at 2.0 volts. Since 
in the de case no ambiguity is introduced by 
possible transient potential differences within 
the phosphor layer, these experiments show that 
electroluminescence can occur at applied volt- 
ages corresponding to about half the band gap 
(3.8 electron volts) of the ZnS phosphor. The 
acceleration-collision theory of electrolumines- 
cence is thus ruled out at low voltage, and since 
no appreciable difference is found to be charac- 
teristic of the electroluminescence at low volt- 
ages, the acceleration-collision mechanism may 
not be important in any voltage range. 


THEORY OF AUGER NEUTRALIZATION OF 
IONS AT THE SURFACE OF A DIAMOND-TYPE 
SEMICONDUCTOR. Homer D. Hagstrum, Bell 
Telephone Laboratories, Murray Hill, New 
Jersey (Received November 21, 1960). 


The two-electron, Auger-type transitions which 
occur when an ion of sufficiently large ionization 
energy is neutralized at the atomically clean sur- 
face of a diamond-type semiconductor are dis- 
cussed. Consideration ofthe basic elements of 
the problem leads to a computing program which 
enables one to calculate the total electron yield 
and kinetic energy distribution of ejected elec- 
trons in terms of a number of parameters. It is 
possible to account for the experimental results 
for singly-charged noble gas ions incident on the 
(111) faces of Si and Ge and the (100) face of Si. 
The fit of theory to experiment is unique in its 
principal features yielding numerical results con- 
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cerning (1) the state density function for the va- 
lence bands of Si and Ge, (2) the energy depend- 
ence of the matrix element as it is determined 
by symmetry of the valence-band wave functions, 
(3) the effective ionization energy near the solid 
surface, (4) energy broadening, and (5) electron 
escape over the surface barrier. Over-all width 
of the valence band is found to be 14-16 ev for both 
Si and Ge. Width of the degenerate p bands is 
5.1 ev in Si, 4.3 ev in Ge. The matrix element 
for p-type valence electrons is 0.3 times that for 
s-type valence electrons. Effective ionization 
energy is 2.2 ev less than the free-space value 
for 10-ev He* ions and decreases linearly with 
ion velocity. Energy broadening is small for 
10-ev ions and increases approximately linearly 
with ion velocity. Probability of electron escape 
is several times that predicted for an isotropic 
distribution of excited electrons incident on a 
plane surface barrier. A general theory of Auger 
neutralization is given in which the conclusions 
of the fit to experiment are interpreted. Investi- 
gation of the matrix element as a Coulomb inter- 
action integral involving wave functions whose 
general characteristics are known but which are 
not explicitly evaluated leads to an understanding 
of its principal dependences on energy and angle. 


VIBRATIONAL SPECTRUM OF GOLD. S. K. 
Joshi and M. P. Hemkar, Physics Department, 
The University, Allahabad, India (Received No- 
vember 28, 1960). 


A theoretical lattice vibration spectrum is 
calculated for gold on the basis of a three-force- 
constant model as elaborated by Bhatia. The 
lattice specific heat is calculated and compared 
with the experimental data. 


CYCLOTRON RESONANCE IN INDIUM ANTI- 
MONIDE AT HIGH MAGNETIC FIELDS. B. Lax, 
J. G. Mavroides, H. J. Zeiger, and R. J. Keyes, 
Lincoln Laboratory, Massachusetts Institute of 
Technology, Lexington, Massachusetts (Re- 
ceived November 28, 1960). 


The room-temperature pulsed magnetic field 
infrared cyclotron resonance data of Keyes and 
co-workers in indium antimonide are interpreted 
using the k-p perturbation technique of Kane, 
which is extended to include the effects of a dc 
magnetic field. The energy levels are calculated 
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for a carrier in a band, in the presence of a mag- 
netic field, which is interacting with neighboring 
bands. From the energy expression, the varia- 
tion of the apparent effective mass with magnetic 
field is calculated. The theory is consistent with 
the room-temperature measurements of Keyes 

et al. as well as with other cyclotron resonant 
results obtained at low temperatures. 


ANNIHILATION OF POSITRONS IN LiH. A. T. 
Stewart* and R. H. March,! Dalhousie Univer- 

sity, Halifax, Nova Scotia (Received November 
23, 1960). 


The angular correlation of photons from posi- 
tron annihilation in LiH and NaH has been meas- 
ured. The data yield a wave-function product 
density distribution much wider than the outer- 
shell electron density around the negative ion. 
This result is in contrast with the observations 
for other alkali halides for which these two dis- 
tributions are much alike. 


*Now at the University of North Carolina, Chapel 
Hill, North Carolina. 

tNow at Clarendon Laboratory, Oxford University, 
Oxford, England. 


CURRENT NOISE AND DISTRIBUTED TRAPS 

IN CADMIUM SULFIDE. James J. Brophy, 
Physics Division, Armour Research Foundation, 
Chicago, Illinois (Received November 30, 1960). 


The high-frequency portion of the current 
noise spectrum observed in lightly doped CdS 
single crystals under uniform 5200A illumina- 
tion is characteristic of electron trapping transi- 
tions to shallow levels. In many crystals the 
noise spectra have a characteristic 1/f trend 
when the electron quasi-Fermi level is not 
located near a discrete trap. The 1/f trapping 
noise observed in one crystal at temperatures 
of 10, 27, and 52°C and for 20 different positions 
of the electron quasi-Fermi level between 0.5 and 
0.3 ev below the conduction band can be repre- 
sented by a single expression of the form 
(1/w)tan“'wr, where 7 is determined by the low- 
frequency turnover of the 1/f trapping noise. 
From these experimental values of 7, trap depths 
are calculated which are in good agreement with 
the positions of discrete trapping levels deter- 
mined from other measurements. Since the low- 
frequency turnovers of the 1/f spectra are thus 


related to the discrete traps, rather than to the 
position of the electron quasi-Fermi level 
directly, it appears that the 1/f noise may not 
be associated with a postulated continuous dis- 
tribution of traps in energy, but rather with a 
dispersion of capture and release times into the 
discrete traps. 


PARAMAGNETIC BEHAVIOR OF METALLIC 
CERIUM AND EUROPIUM. R. V. Colvin, Sigurds 
Arajs, and J. M. Peck, Edgar C. Bain Laboratory 
for Fundamental Research, United States Steel 
Corporation Research Center, Monroeville, 
Pennsylvania (Received November 28, 1960). 


The magnetic behavior of metallic polycrys- 
talline cerium and europium has been studied 
above room temperatures. The measured para- 
magnetic behavior of cerium can be explained 
using the interacting Ce*** ion model (6 = - 50°K) 
and the Van Vleck theory of paramagnetism with 
an additional temperature-independent paramag- 
netic term x, =1.00x10"* g~! cm® resulting from 
the conduction electrons. Europium metal has 
an unusual magnetic behavior. Its paramagnetic 
properties in the solid state definitely cannot be 
explained on the basis of the Eu**+ model. Near 
the melting point metallic europium behaves as 
a collection of weakly interacting Eu* ions. The 
Bohr magneton number of liquid europium is 
very close to that of Eu** ions. 


ABSOLUTE MEASUREMENT OF THE ATOMIC 
SCATTERING FACTORS OF IRON, COPPER, 
AND ALUMINUM. Boris W. Batterman, Bell 
Telephone Laboratories, Murray Hill, New Jer- 
sey, and David R. Chipman and John J. DeMarco, 
Ordnance Materials Research Office, Watertown, 
Massachusetts (Received November 16, 1960). 


The x-ray atomic scattering factors of iron, 
copper, and aluminum were carefully remeas- 
ured to obtain more reliable information on the 
outer electron charge densities in these elements. 
The scattering factors were obtained from meas- 
urements of the integrated Bragg intensities of 
powder samples using monochromatic Mo Ka 
radiation. The intensities were put on an abso- 
lute scale by direct measurements of the power 
in the primary x-ray beam. Extinction, surface 
roughness, and preferred orientation effects 
were shown to be negligible in the samples used. 
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The ratios of the measured scattering factors of 
the three elements agree with those calculated 
from Hartree-Fock theory to within 1%. This 
substantiates the findings of Batterman, and, in 
contrast with the previous results of Weiss and 
DeMarco, indicates that there is no large dis- 
crepancy between the electronic structures of 
copper and iron. The absolute values of the 
measured scattering factors, however, lie about 
4% below theory in the region of low sin@/,. It 
is pointed out that the high theoretical values for 
iron and copper could result from known differ- 
ences in electronic structure between a free 
atom and one in the solid, but that present theory 
probably cannot account for the discrepancy in 
the case of aluminum. 


CLUSTER SIZE IN RANDOM MIXTURES AND 
PERCOLATION PROCESSES. C. Domb and 

M. F. Sykes, Wheatstone Laboratory, King’s 
College, London, England (Received November 
30, 1960). 


It is pointed out that the topological problem 
which arises in the theory of critical concentra- 
tion in ferromagnetic crystals, and discussed 
recently by Elliott et al., is identical with a 
problem arising in the statistical theory of ran- 
dom mixtures. Additional terms of series ex- 
pansions are given, and these enable improved 
estimates of the critical concentration to be 
made. A distinction is drawn between “site mix- 
tures” and “bond mixtures,” the latter having 
been previously considered under the heading of 
“percolation processes.” The critical concen- 
trations for the two problems are usually dif- 
ferent. Some rigorous results regarding the 
critical concentrations are quoted. 


MICROWAVE ZEEMAN EFFECT OF FREE 
HYDROXYL RADICALS. H. E. Radford, National 
Bureau of Standards, Washington, D. C. (Received 
November 7, 1960). 


Paramagnetic resonance absorption at 3-cm 
wavelength has been observed in the products of 
an electric discharge in low-pressure H,O and 
D,O vapor. The spectra are of the electric dipole 
type, and arise from A-type doubling transitions 
in low-lying rotational levels of the free OH 
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and O**D radicals. The theory of the Zeeman 
effect in 711 levels of light diatomic radicals is 
extended to the general intermediate-coupling 
case, and is used for a detailed analysis of the 
observed spectra. Numerical results of this 
analysis include molecular g factors precise to 
within 3 parts in 10°, and magnitudes of the A-type 
doubling intervals in several low rotational levels. 
The measured g factors are compared with theory, 
including small corrections for molecular rota- 
tion, the anomalous spin magnetic moment of the 
electron, and estimated relativistic effects. This 
comparison yields the value 0.66+ 0.01 for the 
molecular matrix element (1 ILy |=), and also 
brings to light serious discrepancies between the 
present experimental results and earlier meas- 
urements of the A-type doubling in OH and OD. 
The paramagnetic resonance spectra also ex- 
hibit hyperfine structure, from which are derived 
the following molecular constants: 


(1/r*), = (0.616 + 0.006) x10" cm=, 


((3 cos?x - 1)/”*) y = -(0.21+ 0.04) x10 cm“, 


y?(0) = (0.184 + 0.004) x10™ cm™, 


which describe the distribution of unpaired elec- 
trons about the H or D nucleus. 


ANALYTICAL METHODS IN HARTREE-FOCK 
SELF-CONSISTENT FIELD THEORY. S. Huzi- 
naga,* Laboratory of Molecular Structure and 
Spectra, Department of Physics, University of 
Chicago, Chicago, Illinois (Received November 
7, 1960). 


Two alternative schemes are proposed for the 
determination of electronic self-consistent field 
(SCF) orbitals in atoms and molecules. They 
are designed to be applied principally to elec- 
tronic configurations consisting of two open 
shells. Both schemes are based upon the idea 
that the SCF orbitals are expanded in terms of 
adequate basis functions but they are different 
in the way of solving the SCF problem. An at- 
tempt is made to rate relative merits of both 
schemes, though they have not yet received any 
actual application. 


*Permanent address: General Education Depart- 
ment, Kyushu University, Fukuoka, Japan. 
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DETERMINATION OF ELECTRON AND POSI- 
TRON HELICITY WITH M@LLER AND BHABHA 
SCATTERING. J. D. Ullman,* H. Frauenfelder, 
H. J. Lipkin, f and A. Rossi, University of 
Illinois, Urbana, Illinois (Received November 16, 
1960). 


The determination of the helicity of electrons 
and positrons from beta decay by means of elec- 
tron-electron (Mgller) and positron-electron 
(Bhabha) scattering is discussed. The theoreti- 
cal background, the apparatus, and the experi- 





Radioisotope ps? 


mental procedure are treated in detail. The 
apparatus included, in addition to the conven- 
tional parts, a beta monochromator with a mo- 
mentum resolution of 16% for the investigation 
of the energy dependence of the helicity. Ex- 
periments were performed with P*, Au’®, RaE, 
and Ga™. In all cases, the helicities P were 
found to be proportional to v/c. The measured 
helicities, in units of v/c and averaged over the 
observed range of energies, are summarized 

as follows. (The errors given are statistical 
only; systematic errors are estimated to be less 
than + 3%.) 

Ga*®® 


Au!% RaE 





Particles e 
660 - 990 
-1.00 +0. 02 


Energy interval in kev 
Helicity P/(v/c) 


*Present address: Centre de Recherche Nucléaire, 
Department de Physique Nucléaire, Université de 
Strasbourg, Strasbourg, France. 

tPresent address: Weizmann Institute of Science, 
Rehovoth, Israel. 

tPresent address: Centro Ricerche Nucleari, Sorin, 
Saluggia, Vercelli, Italy. 


SOME CONSIDERATIONS CONCERNING THE 
NUCLEAR MATRIX ELEMENT 27 1B;;|? IN 

BETA DECAY. Don S. Harmer, Engineering 
Experiment Station, Georgia Institute of Tech- 
nology, Atlanta, Georgia, and Morris L. Perlman, 
Chemistry Department, Brookhaven National 
Laboratory, Upton, New York (Received Novem- 
ber 29, 1960). 


Evidence that the nuclear matrix element 
2) |B;j|? contributes to first-forbidden transi- 
tions may be deduced from the electron capture- 
positron emission ratios observed in these tran- 
sitions and from an analysis of the energy and Z 
dependence of the coefficients of all the nuclear 
matrix elements involved. Tne analysis of these 
coefficients for the capture and for the emission 
processes shows that it is reasonable to expect 
that, except for >) |B,j\’, these transitions would 
be characterized by essentially the same capture- 
positron ratios as those characteristic of allowed 
transitions. Consequently the deviations of the 
observed ratios from the allowed values would 


460 - 810 
-0.98 +0. 03 


+ 


e e e 
1030 - 1300 
+0.99 +0.09 


520 - 950 
-0.75 +0. 03 





be due to contributions from )) IB; -|?; and the 
magnitude of the deviation, in a particular case, 
together with lifetime information may be used 
to compute values for this matrix element and 
for the sum of the others. Values are thus com- 
puted from available experimental data for three 
2--~ 2+ transitions, and these are compared with 
values derived directly from the theory of Rose 
and Osborne. In each case the value of 2B; j |? 
computed from experiment is nearly as large 

as the theoretical )/ |B;j|* value calculated for 
nucleon spin change two and is considerably 
larger than those calculated for spin change zero 
or one; furthermore, the “observed” value for 
the sum of the other matrix elements is con- 
siderably smaller than the largest single matrix 
element calculated for spin change zero or one. 
These results are consistent with the conclusion 
that the nuclear states involved are of such a 
nature that the transitions proceed unhindered 
only via 21; ; |?, 4j=2. One arrives at the same 
conclusion on the basis of the shell-model de- 
scription of these nuclear states. For 2-~0+ 
transitions, which can be effected only by 

2 |B;;|?, the “observed” and theoretical values 
of the matrix element are in good agreement; 
and it is of interest to note that both the theo- 
retical and “observed” values of >) 1B,;\" in 
each of these 2- + 0+ transitions are about the 
same or smaller than the values for the 2- + 2+ 
transition in the same nucleus. 
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SPIN-ORBIT SPLITTING IN NUCLEI DUE TO 
TENSOR INTERACTION. Paul Goldhammer, 
Physics Department, University of Nebraska, 
Lincoln, Nebraska (Received October 31, 1960). 


The effect of the tensor force in nuclei with 
closed shells plus one nucleon has been investi- 
gated using second-order perturbation theory. 

It is found that one can explicitly exhibit the 
spin-orbit splitting due to the tensor force using 
some simple identities. The spin-orbit splitting 
in He® is computed, and found to be 3.4 Mev 
compared with an experimental value of 2.6 Mev. 


SMALL-ORDER SHAPE FACTORS IN In™, P*, 
AND Y”. R. T. Nichols,* R. E. McAdams, and 
E. N. Jensen, Institute for Atomic Research and 
Department of Physics, Iowa State University, 
Ames, Iowa (Received November 7, 1960; re- 
vised manuscript received January 4, 1961). 


The beta spectra of In™*, P**, and Y® have 
been studied closely in an intermediate- image 
beta-ray spectrometer and compared to theoret- 
ical predictions in terms of a linear shape fac- 
tor of the form (1+aW). The values obtained 
for a were (+0.0036+ 0.0021)/mc? for In’, 
(-0.0133 + 0.0011)/mc? for P®, and (-0.0047 
+ 0.0008)/mc? for Y®, all for electron kinetic 
energies from about 200 kev up to near the maxi- 
mum beta energies. Tests were made to give 
indications for spectrometer fidelity. Because 
of the linearity of the shape-factor plots and the 
similarity in energy range, the comparative re- 
sults from In’, P*?, and Y® are taken as a de- 
finite indication that for at least two of these 
activities the shape factors have nonzero slopes, 


irrespective of questions of instrumental fidelity. 


*Present address: Gustavus Adolphus College, St. 
Peter, Minnesota. 


ANALYSIS OF SOME DEUTERON-INDUCED 
REACTIONS IN OXYGEN-18. J. C. Armstrong* 
and K. S. Quisenberry, Radiation Laboratory, 
University of Pittsburgh, Pittsburgh, Pennsyl- 
vania (Received November 21, 1960). 


The reactions O"*(d, t)O'’, O'*(d,d’)O"**, and 
0"*(d, p)O”* are studied using 15-Mev deuterons 
and magnetic analysis of reaction particles. 
Absolute cross sections are determined for all 
reactions studied and the Butler-Born approxi- 
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mation is used to extract reduced widths when 
possible. Angular distributions of triton groups 
corresponding to the ground, 0.871-, 3.846-, 
4.555-, 5.083-, and 5.378-Mev states of O'” are 
obtained. An estimate of the configuration ad- 
mixtures in the O” ground state is made from 
analysis of the reduced widths and indicates the 
presence of a sizable (about 6%) (1f,,.”), com- 
ponent. The experimentally determined admix- 
tures are compared with several theoretical 
estimates. All O” levels observed in the inelas- 
tic deuteron scattering have been previously 
reported—the known 5.01-Mev state is not ob- 
served. The angular distribution of inelastic 
deuterons corresponding to the 1.982-Mev state 
of O” is obtained and comparison of the absolute 
cross section with theory provides an estimate 
of the O”* deformation. Proton groups from 
O"*(d, p)O"* reactions are observed correspond- 
ing to O'* excitations of 0, 1.469, 3.164, 3.948, 
(4.123), (4.586), (4.706), (5.165), 5.45, 5.707, 
and 6.279 Mev , where assignment of the levels 
in parentheses to O** is uncertain. The known 
0.096-Mev state is not observed and the proton 
group corresponding to 5.45-Mev excitation con- 
tains contributions from at least two states. 
Angular distributions leading to the O"* ground, 
1.469-, 3.164-, 3.948-, 5.707-, and 6.279-Mev 
states are obtained and reduced widths extracted. 
The J, values for these angular distributions are 
ambiguous except for the ground-state reaction 
(l, =2) and the 1.469-Mev state reaction (J, =0). 
Analysis of the data suggests that J” (ground 
state) = 5/2* and J"(0.096-Mev state) =3/2*. 
Using parameter values estimated from the O'* 
energy level spectrum or obtained from neigh- 
boring nuclei, a description of this nucleus in 
terms of the strong-coupling unified model agrees 
with the data. 


*Now at Physics Department, University of Mary- 
land, College Park, Maryland. 


14.4-Mev (n,2n) CROSS SECTIONS. L. A. Ray- 
burn, * Argonne National Laboratory, Argonne, 
Illinois and Department of Physics and Astron- 
omy, University of Georgia, Athens, Georgia 
(Received November 14, 1960). 


Cross sections for the (m,2n) reaction have 
been measured at an incident neutron energy of 
14.4+0.3 Mev for 27 nuclides. These measure- 
ments were made relative to the cross section 
for the Cu®*(n, 2n)Cu™ reaction. The relative 
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cross sections were then converted to absolute 
cross sections by using the weighted mean of 
several Cu®(n, 2n)Cu™ reaction- cross-section 
measurements made by other investigators. 


*Present address: Department of Physics and 
Astronomy, University of Georgia, Athens, Georgia. 


AVERAGE RADIATIVE CAPTURE CROSS SEC- 
TIONS FOR 7- TO 170-kev NEUTRONS. J. H. 
Gibbons, R. L. Macklin, P. D. Miller, and J. H. 
Neiler, Oak Ridge National Laboratory, Oak 
Ridge, Tennessee (Received August 17, 1960). 


Measurements of neutron radiative capture 
cross sections in the kev region have been made 
using fast (millimicrosecond) time-of-flight tech- 
niques and a large liquid scintillator tank. Two 
series of measurements have been completed on 
a number of nuclides. These are determinations 
of (1) cross sections relative to that of indium at 
30 kev and at 65 kev for 49 nuclides, and (2) cross 
sections as a function of neutron energy for the 
following nuclei: Br, Nb, Pd, Ag, Cd, In, Sb, I, 
Pr, Sm, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Ta, 
W, Pt, and Au. Curve fits, using the statistical 
model, have been obtained for Br, Nb, Ag, In, 
Sb, I, Pr, Tb, Ho, Tm, Lu, Ta, and Au. The 
results demonstrate the presence of the 2p giant 
resonance near A = 100 predicted by the optical 
model. The average nuclear parameters obtained 
are in good agreement with recent low-energy 
total cross-section results, but are in poor agree- 
ment with earlier results. Possible reasons for 
these disagreements are discussed. 


DECAY OF Na”*. E. L. Robinson, B. T. Lucas, 
and O. E. Johnson, Physics Department, Purdue 
University, Lafayette, Indiana (Received Novem- 
ber 18, 1960). 


The bombardment of high-purity natural mag- 
nesium targets with fast neutrons from the 
Li’(d,n)Be® reaction (E,, 424 Mev) yielded, in 
addition to activities associated with well-known 
isotopes, an activity having an experimental 
half-life of 1.03+ 0.06 sec. Scintillation spec- 
trometer measurements showed a (1.82 + 0.03)- 
Mev gamma ray and a (6.7+ 0.3)- Mev beta group 
decaying with the 1.03-sec half-life. The 6.7- 
Mev beta group was found to be in coincidence 
with the 1.82-Mev gamma ray. The assignment 
of the (1.03 + 0.06)-sec half-life to Na®* produced 


in the reaction Mg**(n, p)Na”* and the proposed 
decay scheme are supported by internally con- 
sistent arguments based on the known character- 
istics of the reaction products and their decay 
modes, the half-life studies using both beta and 
gamma radiation, the features of experimental 
beta and gamma spectra, and nuclear system- 
atics. The decay of the ground state of Na”® 
takes place in part by a (6.7+ 0.3)- Mev beta 
transition to the (2+) first excited state of Mg”®. 
The intensities of the beta transitions to the (2+) 
second excited state and the (0+) ground state 
are less than 0.1 of the intensity of the (6.7+ 0.3)- 
Mev transition. The possible ground-state spin 
and parity assignment for Na” is either 1+, 2+, 
or 3+. A weak theoretical argument against 
spin 1 is presented. The present experimental 
measurements alone will not reduce the ambi- 
guity in the spin assignment. A Na?°— Mg? mass 
difference of 8.5+ 0.3 Mev is derived from the 
known level structure of Mg”® and the beta tran- 
sition energy measured in the present investiga- 
tion. 


MECHANISM OF THE REACTION O"*+p—p+4a 
AT 29 Mev. Oscar C. Kolar, Lawrence Radia- 
tion Laboratory, University of California, 
Berkeley, California (Received October 26, 
1960). 


An expansion cloud chamber containing oxygen 
gas at 1/3 atmosphere pressure was used to 
study the reaction O'*+p—p+4a at a bombard- 
ing energy of 28.9 Mev. Two hundred and twelve 
events were obtained that satisfied the criteria 
of energy and momentum balance. Ninety-one 
of these had all five outgoing prongs visible, 
while the remaining 121 had but four prongs 
visible, the fifth being obscured by the beam. 
Slightly more than half of all the events showed 
the presence of the ground state of Be®. Of 
these, five events showed the presence of two 
Be® nuclei in the ground state. The events ex- 
hibiting the presence of a single ground-state 
Be® were interpreted according to the mecha- 
nism O'*+p—p+2a+Be®; Be®~2a. The possi- 
bility of a compound state was considered. If 
such an intermediate state did occur, it was 
such that it did not obey strictly the predictions 
of the compound nucleus theory. The remaining 
half of the events did not show evidence for any 
intermediate nuclei (with the possible exception 
of the appearance of the 1.4- Mev state of B® in 
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5%, or fewer, of the cases) and could be inter- 
preted only on the basis of the direct quadripar- 
tition of the oxygen nucleus. 


CLUSTER NATURE OF Li’ AND Be’. T. A. 
Tombrello and G. C. Phillips, Rice University, 
Houston, Texas (Received November 17, 1960). 


Measurements of the capture y-radiation pro- 
cesses, mass 3+ a— mass 7+y and nucleon +Li®— 
mass 7++y, give information about the cluster 
structure of the mirror nuclei Li’ and Be’. The 
cluster model predicts that the ground state and 
low excited states of these nuclei should have 
large reduced widths @,,,” for the configuration 
mass 3+a particle and small reduced widths 
6,,,_° for the configuration nucleon + Li®. Scatter- 
ing experiments provide accurate initial, cap- 
turing, wave functions, and an assumption of the 
cluster nature of the final, bound states allows 
the electromagnetic capture cross sections to be 
calculated and compared to experiment. The re- 
duced widths deduced show that @,,,” is large, 
6,4. is small, and that the ground states and 
first excited states of Li’ and Be’ are primarily 
of the two-body cluster form, mass 3 + a@ par- 
ticle. 


ANGULAR YIELD OF NEUTRONS FROM THE 
T(d,n)He* REACTION FOR 6- TO 11.5-Mev 
DEUTERONS. Murrey D. Goldberg and James M. 
LeBlanc, Lawrence Radiation Laboratory, Uni- 
versity of California, Livermore, California 
(Received November 21, 1960). 


The angular yield of monoenergetic neutrons 
from the T(d,)He* reaction has been measured 
with 6.2-, 7.9-, 9.1-, 10.2-, and 11.4-Mev 
deuterons. The neutrons were detected with a 
proton recoil telescope which provided discri- 
mination against neutrons from the T(d,np)T 
breakup reaction. The yield curves are all 
peaked forward, with a second maximum at about 
65° which becomes more pronounced with in- 
creasing energy, and a back-angle rise. A quasi- 
absolute determination of the 0° cross section for 
this reaction is described. The measured yield 
curves are compared to those for the companion 
He*(d,p)He* reaction and a strong similarity in 
shape and magnitude is noted. A simple strip- 
ping model is inadequate to describe these data. 
A distorted-wave calculation is required. 
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TWO-NUCLEON INTERACTION FROM DOUBLET 
SPLITTINGS. F. C. Barker, Australian National 
University, Canberra, Australia (Received No- 
vember 8, 1960). 


The splitting of closely spaced doublet levels 
of spin J+ 4 andJ - 5 is investigated in cases 
which can be approximated by a core of spin J 
and ans nucleon. Allowance is made for the 
Thomas shifts of the levels. Using j-j coupling 
wave functions, the exchange parameters and 
strength of the effective interaction between nu- 
cleons in the nucleus are determined by a least- 
squares fit to the doublet splittings, and some 
improvement is made over results obtained with 
previously accepted parameters. It is suggested 
that measurement of the sign of some doublet 
splittings would clarify the interaction further. 


K~ -MESON SCATTERING IN NUCLEAR EMUL- 
SION. R. D. Hill, J. H. Hetherington, and D. G. 
Ravenhall, Department of Physics, University 
of Illinois, Urbana, Illinois (Received October 
5, 1960). 


Elastic single scatterings of K” mesons by 
nuclei of a dilute G5 photographic emulsion have 
been experimentally measured for energies be- 
tween 30 and 80 Mev. Accurate numerical cal- 
culations of the differential scattering cross 
sections have been carried out by solving a 
Klein-Gordon wave equation for an optical- model 
potential of the Woods-Saxon form. Experimen- 
tal and theoretical results are in good agreement 
both on an absolute and on a relative- variation 
basis if the real part of the K -nuclear potential 
is attractive. The connection with low-energy 
K~-nucleon scattering is discussed. 


a SCATTERING FROM COMPLEX NUCLEI. 
R. M. Edelstein,* W. F. Baker,f and J. Rain- 
water, Columbia University, New York, New 
York (Received November 29, 1960). 


Differential cross sections were measured for 
nm -carbon scattering at 69.5 and 87.5 Mev and 
a” -oxygen scattering at 87.5 Mev from 20° to 
125° extending the technique of Baker, Rainwater, 
and Williams. The energy resolution was suffi- 
cient to measure pure elastic as well as 5- and 
10- Mev inelastic cross sections. The modified 
Kisslinger optical- model equation was used to 
fit the elastic cross-section data. A ,? analysis 





VoLuME 6, NUMBER 3 


PHYSICAL REVIEW LETTERS 


Fesruary 1, 1961 





for the 69.5-Mev carbon data gave a nuclear 
radius parameter 7,=1.05+ 0.02 fermis and a 
fall-off parameter ¢=1.16+0.07 fermis. These 
parameters give good fits to the other data as 
well. An energy dependence in the strength 
parameters for carbon is observed in qualitative 
agreement with prediction.- 


*Present address: Carnegie Institute of Technology, 
Pittsburgh, Pennsylvania. 

tPresent address: Brookhaven National Laboratory, 
Upton, New York. 


STARS FRODUCED BY x CAPTURE INA 
HYDROGEN BUBBLE CHAMBER CONTAINING 
DISSOLVED HELIUM. Michel Schiff,* The Enrico 
Fermi Institute for Nuclear Studies and Depart- 
ment of Physics, University of Chicago, Chicago, 
Illinois, Roger H. Hildebrand, Argonne National 
Laboratory, Argonne, Illinois and The Enrico 
Fermi Institute for Nuclear Studies, University 
of Chicago, Chicago, Illinois, and Clayton Giese, 
Department of Physics, University of Chicago, 
Chicago, Illinois (Received November 28, 1960). 


One-pronged stars produced by 7 capture in 
a hydrogen bubble chamber containing dissolved 
helium have been investigated. The distribution 
of prong lengths in the interval 0.029 to 0.64 g/cm’ 
is presented. About one third of the prongs in this 
interval are found to have a unique range corre- 
sponding to tritons from the reaction 1~ + He* 
H*+n. Some prongs lying beyond the triton peak 
are identified as protons from the reaction 7~ +He* 
~+H'+3n. The fraction of pions producing stars 
is found to be approximately equal to the helium 
concentration. 


*Present address: 7 bis Rue Alexandre Parodi, 
Paris 10, France. 


TOTAL CROSS SECTIONS FOR NEGATIVE 
PIONS ON PROTONS AT 230, 290, 370, 427, 
AND 460 Mev. John C. Caris,* Lester K. Good- 
win, ' Robert W. Kenney, Victor Perez-Mendez, 
and Walton A. Perkins, II,+ Lawrence Radiation 
Laboratory, University of California, Berkeley, 
California (Received November 21, 1960). 


Total cross sections for negative pions on pro- 
tons were measured at laboratory energies of 
230, 290, 370, 427, and 460 Mev. The measure- 
ments were made in the same pion beams as, and 


at energies identical with those of our 1~-p dif- 
ferential scattering experiments. Comparisons 
of the total and differential scattering can be 
made with the dispersion theory at a given energy 
without introducing the systematic errors that 
would normally enter due to uncertainties in the 
parameters of more than one pion beam. The 
measured total cross sections are found to agree 
within statistics with other measured values, and 
with the sums of elastic, inelastic, and charge- 
exchange cross sections measured at this labora- 
tory. The results are: 


E (Mev) Stota] (mb) E (Mev) Stota) (mb) 





27 +2 
28 +2. 


58 +9 
33 +2 
27 +2 


427 +10 
460 +20 


230 +6 
290 +7 
370 +9 


*Present address: Aeronutronic, Newport Beach, 
California. 

t Present address: Physics Department, Hebrew 
University, Jerusalem, Israel. 

t Present address: Lawrence Radiation Laboratory, 
University of California, Livermore, California. 


RELATIVISTIC PARTICLE SYSTEMS WITH IN- 

TERACTION. Leslie L. Foldy, Case Institute of 
Technology, Cleveland, Ohio (Received Novem- 

ber 17, 1960). 


The possibility of covariantly describing a 
system of a fixed number of particles interacting 
directly is explored by attempting a direct “inte- 
gration” of the cummutation relations for the 
inhomogeneous Lorentz group under restrictions 
appropriate to the term “system of a fixed num- 
ber of particles.” By direct interaction is meant 
the fact that interaction between the particles is 
expressed directly in terms of coordinates, 
momenta, and spins for the particles rather than 
through the agency of a mediating field. The in- 
tegration is carried out in considerable generality 
with the assumption that the infinitesimal gene- 
rators of the group have expansions in inverse 
powers of the square of the velocity of light. The 
result coincides with that obtained earlier by 
Bakamjian and Thomas, but the method employed 
yields greater insight into the generality of the 
result, as well as into how further conditions 
beyond covariance, such as the property which 
is here called “separability of the interaction,” 
can be incorporated in the result. The relation- 
ship of the result to the complete reducibility of 
a representation of the inhomogeneous Lorentz 
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group is pointed out. Possible generalizations 
and applications of the procedures here employed 
are discussed. 


TIME-ORDERED GREEN’S FUNCTIONS AND 

PERTURBATION THEORY. M. Muraskin and 
K. Nishijima, Department of Physics, Univer- 
sity of Illinois, Urbana, Illinois (Received No- 
vember 16, 1960). 


A formulation of field theories based on the 
generalized unitarity condition and parametric 
dispersion relations is presented. In the per- 
turbation theory we discuss the connection be- 
tween the present scheme and the Lagrangian 
theory and derive the renormalizability condition 
in our formulation. Finally we show for typical 
processes in the first, second, third, and fourth 
orders that our theory can reproduce the renor- 
malized Feynman perturbation theory. 


TIME-ORDERED GREEN’S FUNCTIONS AND 
ELECTROMAGNETIC INTERACTIONS. K. Nishi- 
jima, Department of Physics, University of Illi- 
nois, Urbana, Illinois (Received November 16, 
1960). 


In the present paper we study various aspects 
of the Ward- Takahashi equations. In perturba- 
tion theory the equivalence between this set of 
equations and the requirement of gauge invariance 
is shown. It is then shown that these equations 
are valid for composite particles as well as for 
elementary particles. Based on our new formu- 
lation the definition of composite particles is 
given, and then we show with the aid of the Ward- 
Takahashi equations that the photon is an ele- 
mentary particle. 


THEORETICAL ASPECTS OF NONLEPTONIC 
HYPERON DECAYS. A. Pais,* Lawrence Radia- 
tion Laboratory, University of California, Ber- 
keley, California (Received November 14, 1960). 


Recent experimental results on nonleptonic 
hyperon decays are studied on the basis of a 
doublet approximation for strong and weak in- 
teractions, with the implied suggestion that this 
higher symmetry may be more easily discernible 
in such reactions in which K particles do not 
occur explicitly. The doublet approximation is 
characterized by a doublet spin J which is 3, 1, 

0 for baryons, 7, K, respectively, and by a K 
spin. It is not necessary to assume that the 
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strong K interactions are weak compared to the 
strong 7 interactions. For the mentioned re- 
actions it is necessary to assume that the strong 
interactions which do not conserve J play a minor 
role compared to those which conserve J. 

The following refinement of the nonleptonic 
AT =1/2 rule is proposed. (T=isotopic spin.) 
The weak nonleptonic interactions consist of two 
parts H®, H™ with AJ=0, 1, respectively. In 
the doublet approximation H® and H™ separately 
conserve parity in the presence of all strong 7 
and K interactions. H® and H™ together do not 
conserve parity, however. In addition to AJ=1, 
H™ should in general satisfy a further constraint, 
but there are classes of graphs for which AJ=1 
is sufficient. 

Current x current structures for H and H™ 
are examined. Results of a foregoing paper can 
be viewed as a special case of the AJ=0, 1 rule. 
The same is true for results obtained by Feldman, 
Matthews, and Salam and by Wolfenstein. The 
considerations of these authors can be extended 
to wider classes of graphs. 

Odd relative helicity and the relation between 
rates for A+p+m~, £+~+p+7° are consequences 
of the AJ=0, 1 rule only. So is the prediction 
that = decay is strongly P-nonconserving. 

The parity properties of H,H™ are sufficient 
conditions. It is a delicate question whether they 
are necessary. For a subset of graphs they are 
not necessary, but this set seems arbitrary. 

If it is assumed that the parity conditions are 
necessary, the schizon scheme is ruled out. 

It is noted that the nonleptonic weak interactions 
may be generated by the strong interactions in 
terms of the following prescription. H™ is gen- 
erated by assuming that the 7 (K) fields have 
small K (7) components. An H™ is generated by 
assuming that the doublets N, (N,) have small 
N, (N,) components; likewise for N, and N,. This 
procedure also generates a nonelectromagnetic 
AT =3/2 interaction. This last coupling is small 
in the sense that it only contributes to Ki2" to 
the extent that the doublet approximation is not 
valid. 


*Permanent address: Institute for Advanced Study, 
Princeton, New Jersey. 


MANDELSTAM REPRESENTATION WITH ANOM- 
ALOUS THRESHOLDS. R. J. Eden, P. V. Land- 
shoff, and J. C. Polkinghorne, Department of 
Applied Mathematics and Theoretical Physics, 
University of Cambridge, Cambridge, England, 
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and J. C. Taylor, * Department of Physics, Im- 
perial College, London, England (Received 
November 28, 1960). 


It is proved that fourth order diagrams provide 
necessary and sufficient conditions for the 
Mandelstam representation to be valid for every 
finite order in perturbation theory. 


*Now at the Department of Applied Mathematics and 
Theoretical Physics, University of Cambridge, Cam- 
bridge, England. 


EFFECTS OF K-MESON INTERACTIONS ON THE 
NUCLEON ANOMALOUS MAGNETIC MOMENT. 
K. Nishimura, Department of Physics, Rutgers 
University, New Brunswick, New Jersey (Re- 
ceived November 23, 1960). 


A previous calculation of the physical nucleon 
wave function in static-source pion theory is ex- 
tended by including the pseudoscalar K-meson 
interactions, the motivation being that this should 
increase the scalar part of the nucleon anomalous 
magnetic moment, thus improving the result of 
previous calculations. In constructing a trial 
function for the physical nucleon, the method of 
moments was used and terms containing up to 
three mesons were included. Calculation shows 
that a too strong K- meson coupling is detrimental 
to the vector part, and that the scalar part can be 
increased approximately by 10% if the K-meson 
interaction is made moderately low. 


INTERMEDIATE VECTOR BOSON AND RADIA- 
TIVE LEPTON DECAY OF THE K MESON. 
Akira Kanazawa* and Masao Sugawara, Physics 
Department, Purdue University, Lafayette, 
Indiana and Katsumi Tanaka,! Physics Division, 
Argonne National Laboratory, Argonne, Illinois 
(Received October 20, 1960). 


K-meson decay into electron, neutrino, and 
photon is analyzed in the lowest order perturba- 
tion with respect to weak and electromagnetic 
interactions, but without making any approxima- 
tion regarding the strong interaction. The weak 
interaction is assumed to be transmitted by a 
single charged intermediate vector meson which 
interacts with the weak current in the conventional 
way. It is pointed out here that a certain angular 
and momentum distribution of decay particles 


could reveal almost unequivocally whether the 
intermediate vector meson does exist. It is 
shown also that other lepton decays of the K 
meson which include 1 mesons and 7 mesons 
cannot be used for the same purpose. 


*On leave of absence from Physics Department, 
Hokkaido University, Sapporo, Japan. 

tTemporary address for 1960-1961: Istituto di 
Fisica Teorica dell’ Universita di Napoli, Mostra 
dell’ Oltremare, Napoli, Italy. 


NONLINEAR ELECTRODYNAMICS IN GENERAL 
RELATIVITY. Asher Peres, Department of 
Physics, Israel Institute of Technology, Haifa, 
Israel (Received November 30, 1960). 


General relativistic field equations are derived 
from a gauge-invariant electromagnetic Lagran- 
gian, which does not involve derivatives of the 
field, nor any charge density, but otherwise is 
completely arbitrary. These equations are ex- 
plicitly solved in the static spherically symmetric 
case, and it is shown that there are solutions 
which are everywhere regular and behave, at 
large distances, like the gravitational and elec- 
tromagnetic fields of a point charge. Some wave- 
like solutions are also derived. 


INTERNAL SYMMETRIES OF STRONG BARYON- 
MESON INTERACTIONS. Yoshio Shimamoto, 
Brookhaven National Laboratory, Upton, New 
York (Received November 3, 1960). 


A simple model considered previously by Pais 
in which the = and A hyperons are regarded as 
a mass-degenerate supermultiplet in the strong 
pion interaction is reconsidered. It is shown 
that recognition of the symmetry exhibited by 
these hyperons in their pionic coupling leads to 
certain prescriptions which may be used to break 
the symmetry via the strong K-meson interac- 
tions. The symmetry reduction schemes de- 
scribed make possible the construction ofa strong 
baryon-meson interaction Hamiltonian which re- 
quires no more than four coupling constants 
(rather than the customary eight) and which in 
no way imposes severe restrictions on the strong 
interactions. Finally, production and scattering 
amplitudes based on the 4-symmetry are dis- 
cussed. 





